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15.  SUBJECT  TERMS 

Respiration  and  cervical  spinal  cord  injury,  intraspinal  microstimulation,  epidural  stimulation,  Neurochip,  phrenic. 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

Unclassified 

88 

USAMRMC 


19b.  TELEPHONE  NUMBER  (include  area 
code) 


2 


Table  of  Contents 


Paqe(s) 


1.  Introduction .  4 

2.  Keywords .  4 

3.  Accomplishments .  5-7 

3a.  Opportunities  for  Training  and  Professional  Development....  5 

3b.  Dissemination  of  Results .  5-7 

3c.  Plans  for  Accomplishment  of  Goals .  7 

4.  Impact .  7 

4a.  Impact  on  Principal  Disciplines .  7 

4b.  Impact  on  Other  Disciplines .  7 

4c.  Impact  on  Technology  Transfer .  7 

4d.  Societal  Impact .  7 

5.  Changes/Problems . 

5a.  Changes  in  Approach .  7-8 

5b.  Anticipated  Problems  and  Proposed  Corrective  Action . 8 

5c.  Changes  Affecting  Expenditures . 8 

5d.  Changes  Affecting  Use  and  Care  of  Animals . 8 

6.  Products . 8 


7.  Participants  &  Other  Collaborating  Organizations 


7a.  Participating  Investigators . 8-9 

7b.  Changes  in  Active  Support/Personnel . 9 

8.  Special  Reporting  Requirements . 

8a.  Ouad  Chart . 9  (Attached  PDF) 

9.  Appendices . 9  (Attached  PDFs) 


3 


1.  INTRODUCTION: 


Spinal  cord  injuries  (SCIs)  within  mid-  to  upper-neck  regions  (i.e.,  at  or  above  the  fourth  cervical  level) 
often  result  in  severely,  if  not  completely,  compromised  breathing  due  to  disrupted  communication  between 
respiratory  centers  in  the  brainstem  and  circuits  in  the  spinal  cord  that  control  ventilatory  muscles.  The 
diaphragm,  which  is  often  considered  the  primary  muscle  of  inhalation,  receives  its  motor  innervation  from 
the  phrenic  nucleus  located  at  spinal  C3-C5.  Interruption  of  respiratory  signals  from  the  brainstem  to  the 
phrenic  nucleus  results  in  uni-  or  bilateral  paralysis  of  the  diaphragm.  Severe  impairment  of  diaphragm 
function  requires  individuals  to  receive  some  form  of  assisted  ventilatory  support  (e.g.,  mechanical,  phrenic 
nerve,  or  diaphragm  pacing).  Even  if  some  individuals  are  weaned,  patterns  of  breathing  rarely  approximate 
normal,  and  pulmonary  infection-related  hospitalizations  and  deaths  remain  common.  A  need  thus  exists  for 
interventions  that  can  lead  to  a  greater  and  more  normal  restoration  of  diaphragm  and  other  respiratory 
muscle  functions.  Studies  being  conducted  under  this  Department  of  Defense  award  are  based  upon 
increasing  evidence  showing  that  electrical  stimulation  of  the  injured  spinal  cord  can  elicit  some  degree  of 
motor  function  (e.g.,  standing,  limb  movements)  in  animals  and  humans  even  with  motor  complete  SCIs.  The 
most  common  strategy  involves  epidural  stimulation  in  which  an  electrical  current  is  delivered  to  an  electrode 
array  placed  on  the  surface  of  the  spinal  cord.  Other  studies  have  shown,  however,  that  electrical  stimulation 
via  tiny  microwires  inserted  directly  into  spinal  cord  tissue  (intraspinal  microstimulation,  ISMS)  may  produce 
more  natural  motor  output.  Epidural,  high  frequency  spinal  stimulation  has  been  shown  to  offer  potential  for 
improving  aspects  of  respiratory  function,  but  it  is  still  dependent  upon  delivery  of  programmed  signals  (i.e., 
open-loop).  We  are  now  proposing  to  explore  the  capacity  of  ISMS  to  improve  diaphragm  function  when 
delivered  to  phrenic  circuitry  in  a  more  closed-loop-like  fashion.  An  innovative  feature  of  our  approach  is  that 
a  custom-designed  microcircuit  (i.e.,  the  Neurochip)  will  be  used  to  capture  electromyographic  (EMG)  signals 
from  unaffected  respiratory  muscles.  Those  EMGs  will  then  be  used  to  trigger  stimulation  of  the  phrenic 
circuit  in  a  physiologically-appropriate  fashion  under  normal  and  more  demanding  ventilatory  conditions  after 
SCI.  The  intriguing  possibility  also  exists  this  strategy  may  lead  to  persistent  recovery  by  directing  functional 
and  anatomical  circuit  remodeling  such  that  functional  improvements  will  no  longer  require  extrinsic 
stimulation.  We  will  also  compare  outcomes  via  epidural  vs.  ISMS.  Overall,  we  are  proposing  the  first 
investigation  of  a  promising  approach  for  stimulating  plasticity  in  respiratory  circuits.  Success  of  this  project 
will  lead  to  a  significant  shift  in  current  approaches  for  managing  respiratory  dysfunction  following  cervical 
SCIs.  Knowledge  obtained  from  this  study  also  may  have  broader  implications  related  to  the  use  of  ISMS  for 
treating  other  aspects  of  SCI  alone  or  in  combination  with  other  interventions. 

2.  KEYWORDS: 

respiration 

cervical  spinal  cord  injury 

intraspinal  microstimulation 

epidural  stimulation 

Neurochip 

phrenic  nucleus 

phrenic  motoneurons 

phrenic  neurograms 

hypoglossal  nerve 

hypercapnia 

hypoxia 

neuroplasticity 

C2  hemisection 
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3.  ACCOMPLISHMENTS: 


3.a.  Opportunities  for  Training  and  Professional  Development 

This  project  has  provided  training  opportunities  for  two  postdoctoral  fellows  (Gonzalez-Rothi,  Streeter) 
and  two  graduate  students  (Mercier,  Sunshine).  In  each  case,  the  individuals  acquired  advanced  technical 
training  in  neurophysiology  and  respiratory  biology.  Individuals  were  also  directly  involved  in  data  analysis 
and  manuscript  preparation. 

3.b.  Dissemination  of  Results 

Two  papers  have  been  published  and  one  is  in  revision  as  cited  under  Item  6  below.  PDF  copies  of  the 
published  manuscripts  are  appended  to  this  document.  An  overview  of  the  paper  being  revised  is  presented 
below  (Item  3.b.2).  The  paper  from  the  Moritz  (Univ.  Washington)  laboratory  is  still  in  preparation  and  was 
reviewed  extensively  in  a  previous  Quarterly  Progress  Report.  Relevant  details  of  the  two  published  and  one 
under  revision  studies  performed  at  the  University  of  Florida  are  as  follows: 

3.b.1.  Details  for  Mercier  et  al.,  2017  publication: 

As  noted  in  the  last  Progress  Report,  our  efforts  have  been  directed  primarily  at  showing  initial  proof  of 
concept  efficacy  of  naturally-triggered  (i.e.,  closed-loop  like)  cISMS  in  the  C2Plx  model.  We  have  recently 
published  our  findings  in  that  regard  in  the  Journal  of  Neurophysiology  (Mercier  et  al.,2016).  A  corroborating 
and  complementary  study  has  also  been  performed  at  the  University  of  Washington  by  our  subcontract 
project  leader,  Dr.  Chet  Moritz,  and  a  paper  describing  those  results  is  now  being  prepared  for  journal 
submission  (Sunshine  et  al.). 

Thirteen  female  Sprague-Dawley  rats 
(255+16  g;  Plarlan,  Indianapolis,  IN)  were 
distributed  between  two  experimental  groups: 

1)  ISMS  after  acute  SCI,  tested  before  and 
immediately  after  lateral  hemisection  of  the 
spinal  cord  at  C2  (C2Hx)  ( n  =  8),  and  2)  ISMS 
after  subacute  SCI,  tested  5-21  days  after 
C2Plx  (n=  5).  All  data  were  collected  in 
terminal  procedures  using  respiratory 
electrophysiology  protocols  approved  by  the 
Institutional  Animal  Care  and  Use  Committee 
at  the  University  of  Florida  and  the  US  Army 
Medical  Research  and  Materiel  Command 
Animal  Care  and  Use  Review  Office.  Technical 
details  are  in  the  appended  manuscript 
attached  to  this  document.  Emphasis  below  is 
on  experimental  design  and  results. 

The  focus  of  these  experiments  was  on 
deter-mining  whether  cISMS-related 
diaphragm  EMG  activity  could  be  triggered  by 
an  endogenous,  spared  respiratory  signal.  For 
the  immediate  proof  of  concept  purposes,  we 
used  EMG  signals  from  the  tongue  to  trigger 
cISMS.  It  is  well  established  that  tongue 
muscle  activity  is  coordinated  with  respiration. 
We  also  investigated  whether  cISMS  would 
result  in  off-target,  non-respiratory  muscle  co¬ 
activation.  One  muscle  of  interest  because  of 
proximity  of  motoneuron  pools  was  the  external  carpi  radialis  (ECR)  muscle  of  the  forelimb. 


jTongue 


Tongue 

CO 


1 1 M  'Ml" 


Right 

Diaphragm 

(RD) 


Diaphragm 

(ID) 


Figure  1.  Cervical  ISMS  activates  the  ipsilateral  diaphragm  A: 
schematic  diagram  illustrates  the  EMG  recording  sites  relative  to  the 
placement  of  the  stimulating  electrode  (S).  The  stimulating  electrode 
was  placed  in  the  immediate  vicinity  of  the  phrenic  motor  nucleus,  and 
ISMS  was  initiated  via  a  trigger  signal  based  on  inspiratory  tongue 
(genioglossus)  EMG  activity.  LI,  left  inter-costal.  B:  pre-C2Hx 
representative  EMG  recordings  from  the  genioglossus,  extensor  carpi 
radialis  (ECR),  and  both  sides  of  the  diaphragm  showing  baseline 
(prestimulation)  activity  and  a  period  of  genioglossus-triggered  cervical 
ISMS  (100  Hz,  200  |iA;  0.3-ms  pulse  duration,  250-ms  train  duration) 
represented  by  a  gray  box  in  the  tongue  trace.  C:  expanded  trace  of 
the  EMG  recording  from  the  left  diaphragm  during  the  period  of  ISMS 
shows  each  stimulus  artifact,  indicated  by  asterisks,  and  each 
subsequent  MUAP.  D:  overlay  of  each  elicited  MUAP,  aligned  by  the 
stimulus  artifact,  demonstrates  constant  latency  and  amplitude. 
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Figure  1A  shows  a  representative  example  of  triggered  cervical  ISMS  (directed  to  the  left  side  of  the 
spinal  cord)  during  the  inspiratory  phase  in  a  spinal-intact  rat.  Note  that  compound  MUAPs  in  the  left 
diaphragm  are  clearly  discernible  (Fig.  1C)  and  are  entrained  to  each  stimulation  pulse  (e.g.,  Fig.  ID).  In 
this  example,  the  average  latency  between  stimulus  artifact  and  MUAP  (motor  unit  action  potential)  peak 
is  2.14+0.05  ms,  which  is  consistent  with  a  synaptic  delay  before  motor  unit  activation.  In  the  absence  of 
spontaneous  inspiratory  EMG  activity,  it  can  be  appreciated  that  ISMS  targeting  the  left  phrenic  motor 
nucleus  evoked  a  marked  response  in  both  the  ECR  and  diaphragm  but  had  minimal  impact  on  the  left 
tongue  or  right  diaphragm. 


To  confirm  that  the  very  small-amplitude  EMG  signals  in  the  tongue,  diaphragm,  and  intercostals 
represented  activation  of  motor  units,  we  employed  a  previously  published  “detectability  index”  for  evoked 
potentials.  This  assessment  indicated  that  ISMS  evoked  a  small  but  statistically  significant  increase  in 
entrained  and  averaged  EMG  activity  in  the  tongue  (5  of  8  animals),  left  intercostal  (8  of  8  animals),  and 
right  diaphragm  (8  of  8  animals).  It  should  be  noted,  however,  that  the  amplitudes  of  these  off-target 
evoked  potentials  were  modest  compared  with  those  of  the  left  diaphragm  and  ECR  muscles. 


Figure  2.  ISMS  after  acute  C2Hx.  A:  illustration  of  the  electrode  placement 
relative  to  the  C2Hx  lesion.  The  solid  dark  bar  represents  the  injury  site;  the 
dashed  lines  represent  severed  pathways.  B:  examples  of  EMG  activity  in  the 
right  and  left  diaphragm  before  C2Hx.  C:  Examples  of  left  and  right  diaphragm 
EMG  activity  after  acute  C2Hx.  Note  that  activity  is  abolished  in  the  left 
diaphragm.  In  B  and  C,  3  spontaneous  breaths  are  shown  followed  by  3  breaths 
during  triggered  cervical  ISMS  (solid  bars;  100  Hz,  200  A).  D:  stimulus-triggered 
averages  of  right  and  left  diaphragm  EMG  activity  obtained  from  the  period  of 
ISMS  high-lighted  by  dashed  boxes  in  B  and  C.  Stimulus-triggered  averages  were 
scaled  to  the  same  values  and  represent  25  stimulus  triggers.  In  D,  the  stimulus 
artifact  is  highlighted  by  the  gray  boxes.  (Figure  3  in  Mercier  et  al.,  2017). 


Experiments  were  next  per-formed 
immediately  after  C2Flx  (acute 
injury).  The  general  experi-mental 
paradigm  is  illustrated  in  Fig.  2A 
Diaphragm  MUAPs  were  first 
evaluated  with  the  spinal  cord  intact 
(Fig.  28).  The  ISMS-evoked  po¬ 
tentials  in  the  spinal-intact  condition 
were  similar  to  those  described  in 
the  preceding  section. 

Acute  C2Hx  abolished  spon¬ 
taneous  inspiratory  EMG  activity  in 
the  ipsilateral  hemidiaphragm  in  all 
animals  (Fig.  2C).  However,  cervical 
ISMS  still  produced  clear  MUAPs  in 
the  ipsilateral  diaphragm  after  the 
acute  lesion  (Fig.  2D).  Collectively, 
the  data  are  consistent  with  a  slight 
reduction  in  cervical  ISMS-evoked 
MUAP  amplitude  in  most  animals 
after  acute  C2Hx.  The  off-target 
impact  of  ISMS  after  C2Hx  (i.e., 
evoked  activity  in  ECR,  tongue,  and 
right  diaphragm)  was  indistinguish¬ 
able  from  that  reported  above  for  the 


spinal-intact  condition  (data  not  shown). 

For  the  last  phase  of  this  study,  we  tested  cISMS  in  rats  that  were  5-21  days  post-C2Hx  (subacute 
injury).  These  temporally  advanced  lesion  experiments  were  performed  to  determine  whether 
degenerative  processes  triggered  by  C2Hx  (e.g.,  axonal  retraction)  prevented  or  mitigated  the  impact  of 
ISMS  on  diaphragm  activation.  Spontaneous  ipsilateral  diaphragm  EMG  activity  was  absent  in  rats  studied 
5-21  days  after  C2Hx  (Fig.  3 A),  and  subsequent  histological  evaluation  of  the  spinal  cord  indicated 
anatomically  complete  hemilesion  in  all  animals  (Fig.  38). 

The  ISMS  procedure  evoked  clearly  discernible  MUAPs  in  the  ipsilateral  (paralyzed)  hemidiaphragm 
(Fig.  3C).  Compared  with  the  acute  injury  group,  the  only  apparent  difference  in  the  ISMS  responses  in 
the  animals  with  the  subacute  lesions  was  a  trend  for  more  variable  and  longer  latencies.  On  average, 
however,  there  were  no  statistical  differences  between  the  groups  for  latency  (acute  C2Hx:  1.69+0.57  ms, 
subacute  C2Hx:  2.50+0.67  ms),  but  the  increase  in  latency  in  the  subacute  group  was  close  to  threshold 
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for  significance  (P_  0.065).  The  MUAP  amplitude  was  variable,  but  with  no  evidence  for  a  difference 
between  the  two  groups  (acute  C2Hx:  1 .62+1 .45  mV,  subacute  C2Hx:  1 .76+2.65  mV,  P_  0.833). 

ISMS  caused  the  expected  negative 
deflection  in  tracheal  pressure  in  the  subacute 
C2Hx  animals,  thus  confirming  that  there  was 
a  biomechanical  impact  of  the  stimulation  on 
breathing  (Fig.  3D).  The  off-target  muscle 
activation  was  similar  to  what  was  obtained  in 
the  spinal-intact  and  acute  C2Hx  animals. 
Thus  very  small,  but  statistically  significant, 
changes  in  EMG  activity  during  ISMS  were 
detected  in  the  tongue  (4/5  animals), 
intercostal  (3/4  animals),  and  right  diaphragm 
(4/5  animals).  The  left  ECR  showed  greater 
EMG  responses  in  five  of  five  animals. 


Stimulus 

Triggered 

Average 


Pre-StirmJafon 


Figure  3.  ISMS  after  subacute  C2Hx.  A:  representative  diaphragm 
EMG  activity  after  subacute  C2Hx.  In  the  example  traces,  genio- 
glossus-triggered  ISMS  was  delivered  during  the  breaths  marked  by 
the  black  bars.  B:  histological  section  of  the  C2  spinal  cord  stained 
with  cresyl  violet.  The  example  demonstrates  an  anatomically 
complete  hemilesion  extending  to  the  midline  of  cervical  cord.  C: 
stimulus-triggered  averages  from  the  ipsilesional  (solid  line)  and 
contralesional  diaphragm  (dashed  line);  data  were  obtained  from  the 
period  indicated  by  the  dashed  boxes  in  A.  These  traces  illustrate 
activation  of  the  diaphragm  ipsilateral  to  the  C2Hx  lesion.  D:  average 
change  in  tracheal  pressure  during  lung  inflation  at  prestimulation 
baseline  and  during  ISMS.  *P  <0.05.  (Figure  5  in  Mercier  et  al.,  201 7) 


While  our  primary  intent  was  to  determine 
whether  ISMS  evoked  diaphragm  activity  during 
the  period  of  stimulation,  we  noted  that 
spontaneous  diaphragm  EMG  activity  often  was 
present  after  the  stimulus  was  turned  off.  Thus,  in 
six  of  eight  animals  after  acute  C2Hx,  the 
hemidiaphragm  that  was  electrically  silent  before 
ISMS  showed  both  tonic  and  inspiratory-related 
activity  during  and  after  the  1-min  period  of 
stimulation  (Fig.  4/4).  In  these  experiments, 
“activity”  was  defined  as  one  or  more  clearly 
discernible  motor  unit  potentials  that  were 
discharging  phasically  during  the  inspiratory 
period.  The  duration  of  the  effect  was  variable 
and  did  not  persist  beyond  2  min  in  any 


experiment.  On  average,  ISMS-induced  spontaneous  diaphragm  motor  unit  activity  lasted  for  43  ±  38 
respiratory  cycles  after  ISMS  was  terminated.  The  first  breath  immediately  after  cessation  of  ISMS  typically 
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Figure  4.  Example  recordings  illustrate  short  term 
potentiation  of  ipsilateral  diaphragm  EMG  activity 
following  ISMS.  A.  tongue  and  ipsilateral  (left) 
diaphragm  EMG  activity  after  acute  C2Hx,  during  ISMS, 
and  immediately  after  cessation  of  ISMS.  Note  that 
clear  phasic  (inspiratory)  activity  can  be  seen  after 
ISMS,  whereas  the  baseline  showed  no  such  activity. 
B.  a  similar  response  can  be  observed  in  a  subacute 
C2Hx  animal.  In  these  examples,  the  ISMS  was 
triggered  by  the  inspiratory  genioglossus  EMG  signal 
and  was  delivered  at  100  Hz  for  1  min  of  respiratory 
efforts  (Figure  6  in  Mercier  et  al.,  2017). 
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demonstrated  the  greatest  amount  of  motor  unit  recruitment,  and  this  was  followed  by  a  gradual  decrease  in 
motor  unit  activity  over  subsequent  respiratory  cycles.  In  addition,  we  noted  an  increase  in  tonic  activity  (i.e., 
activity  persisting  across  the  entire  respiratory  cycle)  after  ISMS  that  followed  a  time  course  similar  to  the 
phasically  active  motor  units.  A  similar  response  to  ISMS  was  observed  in  three  of  five  rats  after  subacute 
C2Hx  injury  (Fig.  46). 

Conclusions 


This  study  has  demonstrated  the  feasibility  of  using  endogenous  medullary  output  to  trigger  ISMS  to 
the  midcervical  spinal  cord.  The  ISMS  paradigm  was  successful  in  recruiting  diaphragm  motor  units  in 
spinal-intact  animals  and  also  after  both  acute  and  subacute  interruption  of  bulbospinal  inspiratory  drive  to 
the  phrenic  motoneuron  pool  following  high  cervical  SCI.  Further  refinement  of  this  approach  may  enable 
development  of  a  “respiratory  neuroprosthesis”  that,  by  virtue  of  the  endogenous  trigger,  could  adapt  to 
temporally  changing  metabolic  demands.  An  interesting  and  potentially  important  observation  that  was 
not  part  of  our  original  hypothesis  was  that  ISMS-induced  inspiratory  bursting  in  the  hemiparetic 
diaphragm  persisted  for  several  respiratory  cycles  after  cessation  of  the  stimulation.  This  occurred  in  the 
majority  of  acute  and  subacute  C2Flx  lesions  and  is  similar  to  the  previously  described  phenomenon  of 
respiratory  short-term  potentiation  (STP)  (Lee  et  al.  2015).  Tongue  EMG  output  enabled  a  rigorous  test  of 
the  fundamental  hypothesis,  and  we  suggest  that  future  work  should  focus  on  additional  potential  sources 
of  respiratory  output  to  serve  as  endogenous  trigger  signals  or  alternative  detection  algorithms.  Future 
investigations  involving  awake,  spontaneously  breathing  animals  will  be  crucial  for  determining  whether 
chronic  ISMS  delivery  will  function  solely  as  a  neuroprosthetic  or  can  be  used  short  term  to  promote 
functional  and  anatomical  neuroplasticity  (Moritz  et  al.  2007)  alone  or  in  combination  with  other 
therapeutic  approaches  leading  to  long-lasting  improvement  in  respiratory  function. 

3.b.2.  Details  for  Gonzalez-Rothi  et  al.,  in  revision  publication: 

Per  the  recommendations  of  the  panel  that  reviewed  our  proposal,  we  have  also  explored  epidural 
stimulation,  and  that  has  led  to  a  second  manuscript  which  has  been  tentatively  accepted  by  the  Journal  of 
Neurophysiology  pending  minor  revisions  recommended  by  the  referees. 

Overview  of  Approach.  Adult  rats  were  anesthetized  and  mechanically  ventilated  and  HF-ES  (300  Hz) 
was  applied  to  the  ventrolateral  C4  or  T2  spinal  cord  ipsilateral  to  C2Hx  for  60-sec  at  currents  from  100-1000 
pA  via  a  tungsten  microwire.  Arterial  blood  pressure,  bilateral  phrenic  nerve  activity,  and  ipsilateral 
hypoglossal  nerve  activity  were  recorded  before  and  after  high  frequency  epidural  stimulation  (HF-ES). 

Phrenic  and  XII  neural  output  following  cervical  HF-ES  in  rats  with  C2Hx  lesions.  Representative  traces  of 
phrenic  and  XII  nerve  output  following  brief  periods  of  C4  epidural  stimulation  are  provided  in  Fig.  5.  These 
examples  first  illustrate  the  recording  configuration  precluded  assessment  of  respiratory  nerve  activity  during 
the  period  of  stimulation  due  to  a  large  stimulus  artifact  that  obviated  the  underlying  signal.  Second,  robust 
short-term  potentiation  of  the  ipsilateral  phrenic  signal  is  evident  following  HF-ES.  The  lowest  current  (100 
pA)  typically  evoked  minimal  if  any  potentiation  of  respiratory  nerve  output,  but  stimulation  at  500  and  1000 
pA  both  resulted  in  a  substantial  increase  in  the  overall  output  of  the  ipsilateral  phrenic  nerve.  Also,  there  was 
no  discernable  impact  on  XII  or  contralateral  phrenic  output,  even  at  the  higher  stimulus  currents. 
Superimposition  of  phrenic  waveform  averages  (not  shown),  obtained  immediately  following  HF-ES,  further 
illustrates  the  potentiation  of  both  phasic  (inspiratory)  and  tonic  phrenic  output  that  occurred  following  the 
higher  stimulus  currents. 

The  average  changes  in  phasic,  tonic  and  total  activity  in  the  phrenic  and  XII  nerve  recordings  following 
C4  HF-ES  are  presented  in  Fig.  6.  In  addition  to  the  mean  and  standard  error,  each  individual  data  point  is 
provided  in  keeping  with  recent  recommendations  for  data  reporting  in  preclinical  research  (Landis  et  al. 
2012).  Evaluation  of  phasic  inspiratory  activity  in  the  left  (ipsilateral)  phrenic  nerve  revealed  a  statistical 
interaction  (P=0.028)  between  the  stimulus  current  {i.e.,  100,  500,  and  1000  pA)  and  the  time  post-injury 
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Figure  5.  Following  a  left-sided  C2  spinal  cord  hemisection  injury,  HF-ES 
was  delivered  to  the  ventrolateral  epidural  surface  of  the  C4  spinal  cord, 
ipsilateral  to  the  side  of  injury.  A)  Representative  compressed  traces 
depict  integrated  phrenic  and  hypoglossal  nerve  activity  prior  to  and 
following  FIF-ES  at  100-,  500-,  and  IOOOuA.  B)  Expanded  traces  depict 
both  raw  and  integrated  phrenic  and  hypoglossal  nerve  activity  at  baseline 
and  following  each  bout  of  stimulation.  Scale  bars  represent  1  min  (A)  and 
1  s  (B). 
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(i.e.,  2-  vs.  12-wks  post-C2Hx).  Thus, 
higher  C4  stimulus  currents  evoked 
increases  in  phasic  ipsilateral  phrenic 
activity  at  12-  but  not  2-wks  following 
C2Hx  (Fig.  6A).  A  similar  response 
occurred  in  the  phrenic  nerve  contralateral 
to  the  lesion  (P<0.05,  Fig.  6B);  however, 
neither  the  time  post-injury  (P=0.051)  nor 
the  time-stimulus  current  interaction 
(P=0.079)  reached  the  threshold  for 
statistical  significance.  Although  there  was 
no  significant  impact  of  time  post-injury, 
contralateral  phrenic  output  did  increase 
with  stimulus  current  (P=0.014). 

The  C4  stimulation  paradigm  also 
influenced  the  medulla  as  confirmed  by 
small  increases  in  XII  motor  output 
following  FIF-ES.  The  relative  changes 
were  small  as  compared  to  the  ipsilateral 
phrenic  signal,  but  the  XII  inspiratory  burst 
(Fig.  6C)  increased  with  the  stimulus 
current  (P=0.027),  and  potentiation  of  XII 
bursting  tended  to  be  greater  at  12-  vs.  2- 
wks  post-C2Plx  (P=0.070). 
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Figure  6.  Average  changes  in  phrenic  and  hypoglossal  nerve  activity  following 
C4  HF-ES.  Changes  in  phasic  (A),  tonic  (B),  and  total  (C)  phrenic  and  hypoglossal 
nerve  activity  following  HF-ES  of  the  C4  spinal  cord,  delivered  at  2  and  12  weeks  post¬ 
injury.  Data  were  calculated  over  the  three  minutes  following  each  bout  of  stimulation 
and  are  represented  as  a  %  change  from  baseline  activity.  HF-ES  delivered  to  the  C4 
spinal  cord  at  12  weeks  post-C2Hx  resulted  in  an  increase  in  phasic  phrenic  activity 
(*p,<0.05  different  from  100pA,  #p<0.05  different  from  500|iA).  No  other  group 
differences  were  observed,  although  dose  dependent  increases  in  amplitude  (percent 
change  above  baseline)  resulted  following  HF-ES,  with  the  greatest  amplitude 
changes  occurring  following  the  highest  stimulation  intensities  (500  and  lOOOpA). 
Note  that  profound  increases  in  tonic  activity  comprised  much  of  the  increase  in  total 
activity  and  that  the  ipsilateral  phrenic  nerve  demonstrated  the  most  dramatic  changes 
in  nerve  amplitude  relative  to  baseline  (Note  the  scale  of  the  y-axes  on  the  right 
phrenic  and  left  XII  plots  is  lower  than  the  axes  for  the  L  phrenic  plots). 


phasic  and  tonic  components)  are  presented  in  Fig.  5G-I.  All 


The  average  changes  in 
tonic  neural  activity  following  C4 
stimulation  are  provided  in  Figs. 
6D-F.  Note  that  the  relative 
magnitude  (%baseline)  of  the 
increased  tonic  dis-charge  was 
considerably  greater  for  tonic  as 
corn-pared  to  phasic  output. 
Both  the  left  (Fig.  5D,  P=0.001) 
the  right  phrenic  nerves  (Fig. 
5E,  P=0.002)  showed  increases 
in  tonic  bursting  as  the  stimulus 
— current  increased.  The  most 
•  robust  potentiation  occurred  in 

,  I  I  the  left  phrenic  nerve  with  in- 
n  rfiltl  creases  of  approximately  500% 

'  of  the  baseline  output  at  the 
highest  stimulation  intensity. 
There  was  no  evidence  for  a 
time-dependent  effect  in  either 
the  left  or  right  phrenic 
neurogram  (i.e.,  tonic  output 
was  similar  at  2-  and  12-wks 
post-C2Plx).  Tonic  XII  discharge 
tended  to  be  greater  following 
C4  stimulation  at  12-  vs.  2-wks, 
but  this  was  not  statistically 
significant  (Fig.  5F,  P=0.083). 
Changes  in  total  phrenic  and  XII 
activity  (i.e.,  including  both  the 
three  recordings  showed  a  progressive 
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increase  in  total  activity  as  the  stimulus  current  was  increased.  Although  the  XII  data  suggest  an  impact  of 
HF-ES  on  medullary  neurons  (Fig.  6),  no  changes  in  inspiratory  burst  frequency  were  observed  following 
stimulation  of  the  cervical  spinal  cord  (Table  1,  P=0.402). 

Phrenic  and  XII  neural  output  following  thoracic  HF-ES.  The  impact  of  thoracic  (T2)  stimulation  on 
phasic,  tonic  and  total  activity  in  left  and  right  phrenic  nerves,  and  the  left  XII  nerve  is  summarized  in  Fig. 

7.  At  the  100  pA  stimulus  current, 
most  rats  showed  a  small 
increase  in  phasic  discharge  in  all 
three  motor  outputs.  The  only 
exception  was  in  the  left  phrenic 
recording  at  the  2-wk  time  point 
where  only  3  of  6  rats  showed  an 
increase  in  phasic  discharge  (Fig. 
7A).  The  magnitude  of  poten¬ 
tiation  increased  with  stimulus 
current  in  the  left  phrenic  (Fig.  7A, 
P=0.004),  right  phrenic  (Fig.  7B, 
P=0.011)  and  XII  recordings  (Fig. 
6C,  P=0.021).  Plowever,  time- 
depend-ent  effects  ( e.g .,  2-  vs. 
12-wks)  on  phasic  output  were  not 
detected  for  any  of  the  recordings 
(all  P>0.265). 

Tonic  output  was  potentiated 
in  all  three  nerves  following  T2 
stimulation  in  nearly  all  ex¬ 
periments  at  both  the  2-  and  12- 
wk  time  points  (Figs.  7D-F).  In  the 
left  phrenic  signal,  tonic  discharge 
increased  with  stimulus  intensity 
(P=0.004),  and  reached  approxi¬ 
mately  500%  of  baseline  at  the 
highest  stimulus  current.  In¬ 
creases  in  tonic  discharge  in  the 
right  phrenic  and  left  XII  record¬ 
ings  were  typically  in  the  range  of 
50-100%  (Fig.  6E-F).  Time  dependent  effects  on  the  tonic  response  to  T2  stimulation  were  not  statistically 
significant  for  any  of  the  neural  outcomes  assessed  (all  P>0.557).  The  magnitude  by  which  total  activity 
was  potentiated  increased  with  stimulus  intensity  in  all  recordings  (all  P<0.007),  and  there  were  no 
significant  differences  in  the  total  activity  between  the  two  time  points  (all  P>0.075). 

Similar  to  the  cervical  PIF-ES  data,  changes  in  inspiratory  burst  frequency  were  not  observed  following 
stimulation  of  the  thoracic  spinal  cord  (Table  1,  P=0.281). 

Cardiovascular  parameters  during  and  following  HF-ES.  The  heart  rate  (FIR)  remained  stable  during 
all  procedures,  with  the  exception  of  C4  stimulation  at  12-wks  following  C2Flx.  Under  those  conditions, 
increases  in  FIR  ranged  from  1-14%  of  baseline  (Fig.  7C),  and  the  mean  increase  was  statistically 
significant  (P<0.001).  The  mean  arterial  pressure  (MAP)  increased  during  C4  stimulation  at  500pA 
(P=0.033)  and  1000  pA  (P=0.001),  but  with  no  evidence  for  a  time-dependent  change  between  2-  and  12- 
wks  (P-0.259).  The  MAP  also  increased  during  T2  stimulation  at  1000  pA  (P=0.034);  no  time-dependent 
changes  were  detected.  At  both  2-  and  12-wks,  both  FIR  and  MAP  returned  to  baseline  values  following 
cessation  of  PIF-ES  regardless  of  the  level  of  stimulation  (all  P>.627). 

Arterial  blood  was  assessed  for  PaC02,  Pa02,  and  pH  values  prior  to  and  following  the  stimulation 
protocol.  No  differences  between  groups  were  detected  at  either  time  point  (all  P>0.38).  All  rats  remained 


jj  Iflo 


u: 


Is" 

£  iso 


rjirfi  Mrji  li]_ 


ij  S _ Ufa- 1.4- . 


IdOpiA  500  jA 


RIGHT  PHUEMC 


r!-rii  rTM  rll 


rji4n  rjii f i  rjirfi 


rfrii  rHril  rr 


IOOuA  50D|jA  1000mA 


P  40 0 


riiiii  rfnii  fil^l 


Jixfc  rfirii 


-jjillL 


rFrfi  i+i 


iA 


IOCKjA  SOOlA  lOOOu-A 


Figure  7.  Average  changes  in  phrenic  and  hypoglossal  nerve  activity  following 
T2  HF-ES.  Changes  in  phasic  (A),  tonic  (B),  and  total  (C)  phrenic  and  hypoglossal 
nerve  activity  following  HF-ES  of  the  T2  spinal  cord,  delivered  at  2  and  12  weeks 
post-injury.  Data  were  calculated  over  the  three  minutes  following  each  bout  of 
stimulation  and  are  represented  as  a  %  change  from  baseline  activity.  A  dose- 
dependent  increase  in  amplitude  (percent  change  above  baseline)  was  observed  in 
all  nerves  (*p<0.05  different  from  lOOpA,  #p<0.05  different  from  500  pA)  across  all 
outcomes  (phasic,  tonic,  and  total  nerve  activity),  with  the  greatest  amplitude 
changes  occurring  following  the  highest  stimulation  intensities.  Note  that  profound 
increases  in  tonic  activity  comprised  much  of  the  increase  in  total  activity  and  that  the 
ipsilateral  phrenic  nerve  demonstrated  the  most  dramatic  changes  in  nerve  amplitude 
relative  to  baseline  (Note  the  scale  of  the  y-axes  on  the  right  phrenic  and  left  XII  plots 
is  lower  than  the  axes  for  the  L  phrenic  plots).  No  significant  group  differences  were 
observed  in  any  nerve  in  any  of  the  outcomes  assessed. 
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well  oxygenated  throughout  the  experimental  protocol  with  Pa02  values  greater  than  98  mmHg  in  all 
cases. 

The  impact  of  the  serotonergic  receptor  agonist  quipazine  on  ipsilateral  phrenic  output  before  and  after 
HF-ES.  Serotonergic  neurons  innervate  the  spinal  phrenic  circuitry  and  serotonin  is  a  powerful  stimulant  of 
respiratory  motor  output  following  SCI.  In  a  subset  of  experiments,  we  examined  the  impact  of  i.v.  delivery 
of  the  5-HT  receptor  agonist  quipazine  on  ipsilateral  phrenic  motor  output  before  and  following  HF-ES. 
Consistent  with  prior  studies  of  the  5-HT  precursor  5-Hydroxytryptophan  (Ling  et  al.  1994;  Mitchell  et  at. 
1992;  Zhou  and  Goshgarian  2000),  i.v.  quipazine  triggered  a  robust  increase  in  both  the  tonic  and  total 
output  of  the  ipsilateral  phrenic  nerve  at  both  2-  and  12-wks  after  C2Hx  (all  P<0.028).  The  relative 
magnitude  of  the  acute  increases  following  quipazine  exceeded  the  short-term  potentiation  of  ipsilateral 
phrenic  output  that  was  induced  by  epidural  stimulation  prior  to  quipazine.  No  time-dependent  differences 
in  the  relative  increase  in  ipsilateral  phrenic  bursting  induced  by  quipazine  (%  baseline)  were  detected 
across  2-  vs.  12-wks  for  either  phasic,  tonic,  or  total  phrenic  output  (all  P>0.44). 

Following  i.v.  quipazine,  C4  and  T2  HF-ES  potentiated  tonic  phrenic  output  at  both  2-  and  12-wks 
post-injury  (Fig.  9),  and  the  magnitude  of  potentiation  was  dependent  upon  the  stimulus  intensity  (all 
P<0.014).  However,  quipazine  pretreatment  did  not  enhance  the  relative  magnitude  of  tonic  potentiation 
following  stimulation  when  compared  to  the  baseline  response  ( e.g .,  without  potentiation  without  quipazine 
treatment;  P>0.88  for  all  groups).  Conversely,  the  relative  increase  in  phasic  inspiratory  phrenic  bursting 
that  followed  epidural  stimulation  was  blunted  following  quipazine  (i.e.,  reduced  potentiation,  all  P>0.092 
vs.  pre-stimulus  baseline)  in  all  but  the  12-wk  T2  group.  Thus,  the  large  increase  in  overall  tonic  output  at 
this  quipazine  dose,  appears  to  have  largely  precluded  the  ability  to  detect  short-term  changes  in 
inspiratory  phrenic  bursting. 

Conclusions 


This  is  the  first  study  to  examine  the  impact  of  high  frequency  epidural  stimulation  (HF-ES)  on  phrenic 
motor  output  following  incomplete  cervical  spinal  cord  injury.  While  a  detailed  head-to-head  comparison  with 
ISMS  is  premature  since  both  approaches  have  yet  to  be  optimized,  we  have  shown  that  (High  Frequency- 
Electrical  Stimulation)  HF-ES  of  the  ventrolateral  C4  or  T2  spinal  cord  can  cause  an  increase  in  phrenic 
discharge  that  persists  after  stimulation.  The  short-term  potentiation  of  ipsilateral  phrenic  bursting  illustrates 
the  potential  for  spinal  cord  stimulation  to  trigger  respiratory  plasticity,  and  this  is  a  concept  that  could  be 
useful  in  designing  neurorehabilitation  paradigms  (Kasten  et  al.  2013;  McPherson  et  al.  2015).  However,  we 
also  observed  large  increases  in  tonic  phrenic  output  with  epidural  stimulation  which  indicates  that  HF-ES  is 
unlikely  to  have  application  as  a  “respiratory  neuroprosthesis”  for  activation  of  the  respiratory  muscles  after 
incomplete  SCI. 

3.b.3.  Details  for  Streeter  et  al.,  201 6  publication: 

Lastly,  as  part  of  Major  Task  1  in  which  we  are  electrophysiologically  mapping  changes  in  phrenic 
circuitry  following  cervical  SCI,  we  have  just  published  another  paper  in  the  Journal  of  Neurophysiology  that 
describes  a  method  for  identifying  multi-electrode  array  recording  sites.  Several  groups  have  advanced  the 
hypothesis  that  mid-cervical  spinal  interneurons  can  modulate  phrenic  motoneuron  excitability  and  thereby 
influence  the  neural  control  of  the  diaphragm.  A  significant  hurdle  in  testing  that  hypothesis,  however,  is  the 
difficulty  of  studying  the  “functional  connectivity”  in  diffuse  spinal  cord  networks.  A  major  challenge  is  to 
simultaneously  record  numerous  cells.  This  can  be  addressed  with  multi-electrode  array  (MEA)  technology. 
But  histological  identification  of  each  recording  location  (versus  the  electrode  track)  while  preserving  tissue 
integrity  poses  a  further  challenge. 

Basic  Experimental  Approach.  Thus,  the  initial  thrust  of  the  current  work  was  modification  and  validation 
of  a  silver  labeling  technique  to  enable  post-recording  deposition  of  a  small  amount  of  silver  (i.e.,  for 
histological  marking)  from  the  tip  of  each  electrode  in  an  MEA.  Additionally,  we  developed  an  electrical  circuit 
to  enable  the  use  of  small  currents  (100  nA)  for  silver  electroplating  and  deposition  to  prevent  tissue  and 
electrode  damage  associated  with  high  levels  of  current.  Using  the  electroplated  MEA,  we  recorded 
discharge  from  ensembles  of  neurons  in  the  mid-cervical  (C3-C5)  spinal  cord  in  adult  rats  and  demonstrated 
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a  practical  application  of  this  technique  by  “matching”  the  array  electrodes  to  the  corresponding  anatomical 
locations  marked  by  silver. 

Most  experiments  (n=12)  were  performed  with  untreated,  spinal  intact  rats.  A  subset  of  rats  (n=2) 
received  a  cervical  spinal  injury  (C3/C4  lateralized  contusion,  force:  205  KD,  displacement:  1225  pm;  Infinite 
Horizon  pneumatic  impactor,  Precision  Systems  &  Instrumentation)  using  published  methods  from  our  group. 
Spinal  injured  rats  were  allowed  to  recover  for  12  weeks  prior  to  electrophysiology.  Histological  results  from 
spinal  injured  animals  were  used  to  determine  effectiveness  of  silver  labeling  and  to  compare  micro-motor 
depths  relative  to  silver  labeling  depths.  Since  there  is  evidence  suggesting  neurophysiological  properties 
such  as  connectivity  may  be  altered  following  SCI,  neurophysiological  analysis  was  limited  to  spinal  intact 
animals  (n=4).  See  appended,  published  manuscript  for  methodological  details.  Images  below  are  from  the 
published  manuscript  although  figure  numbers  in  text  differ  from  those  in  the  publication. 

Optimization  of  Silver  Deposition  from  ME  As  and  Subsequent  Histological  Development  to  Enable 
Visualization  of  Recording  Locations.  We  first  assessed  the  impact  of  two  critical  variables  on  the  histological 

appearance  of  silver  in  the  spinal  cord: 

1)  the  duration  of  deposition  current,  and 

2)  the  duration  of  histological 
development.  Consistent  and  discrete 
labeling  at  a  resolution  between  -25-50 
pm  was  obtained  when  current  used  to 
deposit  silver  was  applied  for  50  sec  and 
tissue  was  developed  for  7-10  minutes 
(Figures  8A  and  B).  When  deposition 
current  of  less  than  50  sec  was  applied 
few  apparent  labels  were  detected, 
presumably  due  to  an  inadequate 
amount  of  silver  ions  necessary  for 
visualization.  In  contrast,  deposition  of 
80-100  sec  resulted  in  a  “halo  effect” 
which  increased  the  silver  labeled  area 
to  -100-150  pm  (Figures  8C  and  D), 
most  likely  due  to  a  broad  distribution  of 
deposited  silver.  Similarly,  histological 
development  of  less  than  7  minutes  was 
not  long  enough  to  “develop”  the 
deposited  silver  required  to  visualize  the 
label;  whereas  prolonged  exposure  to 
the  developing  solution  tended  to  darken 
large  areas  of  the  section  thus 
confounding  differentiation  due  to 
insufficient  contrast  (data  not  shown). 

We  next  verified  that  sites  of  silver  deposition  represented  the  location  of  the  tip  (i.e.,  recording  location) 
of  a  given  electrode.  In  a  subset  of  animals  (n=2  with  silver  processing,  n=1  unprocessed  positive  control), 
two  of  the  sixteen  silver  plated  electrodes  were  coated  with  the  fluorescent  lipophilic  dye,  Dil  to  visualize  the 
tracks  associated  with  the  electrode  insertion.  While  Dil  was  readily  detected  in  unprocessed  tissue  (data  not 
shown),  Dil  could  not  be  detected  after  the  histological  processing  required  for  silver  labeling.  Therefore,  we 
conclude  these  two  techniques  are  not  compatible.  However,  in  several  histological  sections  an  “electrode 
track”  could  be  identified  terminating  at  the  silver  labeled  site  (Figure  8E).  This  is  consistent  with  the 
interpretation  that  the  silver  labeled  site  is  located  at  the  tip  of  the  electrode. 

Another  objective  was  to  determine  if  the  silver  labeling  method  could  be  coupled  with  fluorescent 
immunohistochemistry  (IHC)  procedures.  Accordingly,  we  performed  IHC  with  two  markers  commonly  used 
in  the  central  nervous  system  (GFAP  and  NeuN)  either  prior  to,  or  immediately  after  histological  processing 
for  silver  development.  These  markers  were  chosen  due  to  their  ability  to  label  astrocytes  and  neurons. 


Duration  of  Histological  Development 


Figure  8.  Representative  images  of  silver  labeling.  A-D:  Representative 
histological  sections  from  the  cervical  spinal  cord  containing  silver  labeled  sites 
counterstained  with  cresyl  violet  and  high  resolution  images  of  the  boxed  areas. 
These  images  depict  the  variability  of  silver  labeling  at  different  durations  of 
silver  deposition  (50-100  sec)  and  histological  development  (7-10  minutes).  E: 
Images  of  a  cervical  spinal  section  depicting  an  electrode  track  (black  arrows) 
coursing  through  the  spinal  tissue  and  terminating  at  the  silver  labeled  site.  CC: 
central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  A-D:  0.5  mm  and  50 
pm  (callouts);  E:  0.5  mm,  100  pm,  and  50  pm,  respectively. 
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Fluorescent  staining  was  detected  in  both  cases,  but  was  more  robust  when  IHC  was  performed  after  the 
silver  staining  procedures  (data  not  shown),  indicating  IHC  (at  least  for  these  commonly  used  markers)  is 
compatible  with  the  histological  processing  steps  necessary  to  “develop”  the  silver.  Using  this  staining  order, 
we  determined  the  location  of  deposited  silver  relative  to  NeuN  positive  cells.  Representative 

photomicrographs  from  a  C4  spinal 


A. 


OH 
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Figure  9.  Silver  labeling  coupled  with  fluorescent  immunohistochemistry. 

A:  Representative  cervical  spinal  section  of  positive  silver  labeling  (brown/black) 
and  high  resolution  image  of  the  boxed  area.  B:  Fluorescent  labeling  of  neurons 
stained  with  NeuN  (green)  and  nuclei  stained  with  DAPI  (blue)  of  the  same 
section  presented  in  panel  A,  and  high  resolution  image  of  the  boxed  area 
illustrating  the  silver  labeled  site  was  in  close  proximity  of  a  NeuN  positive  cell. 

CC:  central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  0.5  mm  and  50 
pm  (callouts). 

locations  of  the 

recording  electrodes,  we  obtained  positive  labeling  in  35  of  39  attempts.  This 
indicates  that  the  silver  labeling  method  can  be  used  to  identify  the  anatomical 
locations  of  array  electrodes  with  a  success  rate  of  approximately  90%.  Positive 
silver  labeling  was  identified  in  the  spinal  gray  matter  between  C3-C5.  Silver  labeled 
sites  were  identified  between  lamina  IV-X,  with  the  greatest  number  found  in  lamina 
VII  (Figure  10A).  The  locations  of  each  silver  labeled  site  in  spinal  intact  and  SCI 
animals  was  plotted  (unilaterally  for  simplicity)  according  to  their  anatomical 
locations  (Figure  10B). 

The  success  of  the  discrete  silver  labeling  afforded  the  chance  to  compare  the 
intended  recording  location  (i.e.,  the  micro-motor  coordinates  used  during  the 
neurophysiology  experiment)  with  actual  the  location  of  the  labeling  (Figure  10B,  C). 
Each  silver  labeled  site  was  measured  from  the  dorsal  surface  of  the  cord.  Based 
on  estimates  of  prior  investigations  regarding  tissue  shrinkage  during 
paraformaldehyde  fixation  and  subsequent 
tissue  processing,  the  measured  histological 
depth  of  each  silver  labeled  site  was  adjusted 
by  10%.  This  analyses  indicated  that  the  micro¬ 
motor  depth  is  likely  to  overestimate  the  actual 
depth  of  the  electrode,  and  perhaps  even  more 
importantly  that  this  relationship  is  altered  by 
SCI.  Linear  regression  analysis  indicated  a 
significant  relationship  between  recording 
coordinates  and  histological  staining  in  spinal- 
intact  animals  (motor:  1.6  +  0.04  mm; 
histological:  1.2  ±  0.05  mm;  P  =  0.0003),  but 


section  depicting  positive  silver  labeling 
in  the  intermediate  gray  is  shown  in 
Figure  9A.  Fluorescent  labeling  of  the 
same  section  stained  with  NeuN  and 
DAPI  are  shown  in  Figure  9B.  High 
magnification  images  of  the  silver 
labeled  site  suggest  the  recording 
electrode  was  in  close  proximity  to  a 
NeuN  positive  neuron  (Figure  9B). 
These  qualitative  results  demonstrate 
that  the  silver  labeling  technique  can  be 
coupled  with  fluorescent  IHC. 

Using  Silver  Labeling  to  Identify 
Anatomical  Locations  of  ME  As.  In  six 
animals  (n=4  spinal  intact  and  n=2  SCI), 
we  used  the  optimized  deposition  and 
development  parameters  discussed 
previously  to 
label  select 
electrodes. 

Overall,  using 
this  approach 
to  identify  the 
anatomical 
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Figure  10.  Electrode  and  silver  labeled  depth.  A:  Representative 
C4  section  of  silver  labeled  sites  in  spinally  intact  (•)  and  spinally 
injured  (o)  animals  projected  onto  one  side  of  the  cord  for  simplicity. 
Each  lamina  was  shaded  according  to  the  number  of  silver  labeled 
sites  within  that  lamina.  B:  Scatter  plot  of  the  micro-motor  depths 
versus  the  depths  of  the  corresponding  silver  labeled  sites  in  spinally 
intact  (•)  and  spinally  injured  (o)  animals.  Linear  regression  analysis 
was  applied  to  determine  the  line  of  best-fit  and  linear  equation  for 
each  group.  The  line  of  identity  is  displayed  to  indicate  the  location 
where  micro-motor  and  silver  labeling  depths  are  equal.  C:  Normalized 
depth  calculated  as  the  difference  between  the  micro-motor  depth 
silver  labeled  site  for  spinally  intact  (•)  and  spinally  injured  (o)  animals. 
*  PcO.OOl  Unpaired  T-test. 
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not  after  chronic  SCI  (motor:  1 .4  +  0.07  mm;  histological:  0.6  ±  0.05  mm;  P  =  0.3014;  Figure  10B).  The  depth 
of  the  silver  labeled  site  differed  from  the  micro-motor  coordinate  in  both  groups  (Figure  10C),  but  was 
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significantly  increased  following  SCI  (0.8  ±  0.07) 
compared  to  spinal  intact  animals  (0.4  +  0.04). 
These  data  highlight  the  need  to  histologically 
identify  multi-electrode  recording  sites  especially 
following  experimental  conditions  such  as  SCI 
when  tissue  fibrosis  and  scarring  can  be 
expected  to  alter  electrode  movement  within  the 
CNS.  Future,  studies  can  apply  the  linear  fit  as  a 
proxy  to  more  accurately  calculate  the  “actual 
depth”  of  the  electrode  tips. 

Coupling  Spinal  Discharge  with  Anatomical 
Locations.  All  silver  electroplated  electrodes 
could  discriminate  single  units,  indicating  that  the 
pretreatment  had  a  minimal,  if  any,  impact  on  the 
ability  to  record  and  discriminate  extracellular 


Representative  G4  section  containing  a  silver  labeled  site 
counterstained  with  cresyl  violet  and  high  resolution  image 
containing  the  silver  labeled  site  located  in  the  medial 
aspect  of  the  ventral  horn  in  lamina  IX  (callout).  B: 
Corresponding  raw  spinal  discharge  and  integrated  phrenic 
output  during  baseline,  hypoxia,  and  post-hypoxia  depicting 
discharge  during  the  inspiratory  phase.  C:  Integrated 
phrenic  output,  raw  neuronal  discharge,  and  sorted  spikes 
(waveform)  during  hypoxia.  D:  Cycle-triggered  histogram 
during  50  consecutive  breaths  overlaid  with  the  average 
integrated  phrenic  waveform  during  hypoxia.  E:  Spike- 
triggered  average  of  the  raw  and  rectified  phrenic  nerve 
depicting  a  delay  of  0.76  ms.  CC:  central  canal;  DH:  dorsal 
horn;  VH:  ventral  horn.  Scale  bars:  A  0.5mm  and  50pm 
(callout);  B,  C:  1 .5  sec. 

signals.  Three  silver  labeled  sites  were  identified 
within  lamina  IX  of  the  ventral  horn  (Figure  10A 
and  11  A).  Spike-triggered  averaging  (STA)  of  the 
raw  and  rectified  phrenic  nerve  activity  in  relation 
to  neuronal  discharge  produced  a  distinct  peak 
with  an  average  lag  time  of  0.45  ±  0.16  ms  (Figure 
1 1 E),  and  therefore  indicated  that  the  recorded 
discharge  was  from  phrenic  motoneurons.  All  three 
of  these  cells  were  active  primarily  during  the 
inspiratory  phase  (Figure  1 1B  and  C),  which  is  the 
typical  firing  pattern  of  phrenic  motoneurons  in  this 
preparation.  The  silver  labeling  data  verified  that 
the  recordings  were  from  the  region  of  the  phrenic 
motor  nucleus.  This  close  matching  between 
neurophys-iological  data  and  anatomical  data 
illustrates  the  effectiveness  of  the  silver  labeling 
method. 

The  remaining  silver  labeled  sites  were  located 
between  laminae  IV-X.  Spike-triggered  averaging 
provided  no  evidence  of  discharge  synchrony  in 
relation  to  phrenic  motor  output;  thus  these  28 
recordings  were  considered  to  represent 
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Figure  12.  Mapping  anatomical  locations  of  mid-cervical 
spinal  neurons.  A:  Representative  photomicrograph  of  a  C4 
section  containing  a  silver  labeled  site  (from  electrode  nine) 
counterstained  with  cresyl  violet  and  high  resolution  image  of  the 
silver  labeled  site  (callout).  B:  Camera  Lucida-style  drawing  of  the 
cervical  spinal  cord  summarizing  the  anatomical  locations  of  eight 
silver  labeled  sites  obtained  in  one  animal.  C:  Integrated  phrenic 
motor  output  and  mid-cervical  spinal  discharge  on  the  left  and  right 
hemi-cord  corresponding  to  identified  silver  labeled  sites  in  panel  B 
during  baseline  and  hypoxia.  Electrode  one  recorded  phrenic 
motoneuron  discharge  and  the  remaining  electrodes  recorded 
interneuron  discharge.  D:  Representative  C4  section  summarizing 
the  total  number  of  silver  labeled  interneurons  within  each  lamina 
in  spinal  intact  animals.  E:  Average  discharge  frequency  (Hz)  of 
spinal  interneurons  within  each  lamina  represented  in  panel  D.  CC: 
central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  A,  B: 
0.5  mm  and  50  pm  (callout);  C:  1 .5  sec. 


interneurons  (see  Figure  7  in  appended 

publication).  At  baseline,  most  recorded  interneurons  (18/28  or  66%)  fired  tonically  throughout  the  respiratory 
cycle.  These  tonically  discharging  cells  were  not  restricted  to  a  particular  lamina,  but  were  instead  recorded 
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throughout  the  cervical  grey  matter.  A  smaller  proportion  of 
interneurons  (5/28  or  18%)  primarily  fired  during  the 
inspiratory  phase,  and  these  cells  were  found  in  lamina  VI 
(n=1),  VII  (n=3),  and  X  (n=1).  During  hypoxia,  two  of  these 
neurons  (in  lamina  VII  and  X)  switched  to  a  tonic  firing 
pattern,  but  then  resumed  an  inspiratory  pattern  post¬ 
hypoxia.  Due  to  phase  switching,  the  proportion  of  tonic 
firing  interneurons  increased  to  82%  (23/28)  during  hypoxia 
and  remained  elevated  post-hypoxia  (20/28  or  71%). 


Figure  13.  Cross-correlation  analysis  of  mid-cervical 
spinal  interneurons.  A:  Representative  cross¬ 
correlations  obtained  from  three  pairs  of  neurons 
depicting  a  significant  central  peak  (left),  1.8  ms  offset 
peak  (middle),  and  7.8  ms  offset  peak  (right).  The 
number  of  trigger  spikes  for  each  correlation:  43,891; 
39,208;  3,928;  and  the  detectability  index  of  each 
correlation:  45.2;  4.3;  5.3,  respectively.  Black  lines 
plotted  at  0.5  and  at  2.9  ms  to  indicate  how  data  was 
grouped  by  latency  (see  panel  B).  B:  A  histogram  of  the 
latency  relative  to  the  trigger  for  all  correlations  with  a 
significant  peak.  The  first  bar  in  the  histogram  is  a  count 
of  the  central  peaks  (0  ms  latency)  and  each  successive 
bar  in  the  plot  represents  counts  obtained  in  an 
increment  of  0.2.  Yellow  bars  plotted  at  0.5  and  at  2.9  ms 
indicate  how  data  was  grouped  based  on  latency.  Black 
lines  are  2nd  order  polynomial  fit  of  the  data.  C:  The  sum 
of  significant  cross-correlations  with  latencies  <  0.4  ms, 
between  0.6-2. 8  ms,  and  >  3.0  ms.  Percentages  reflect 
the  proportion  of  positive  correlations  out  of  the  total 
possible  (n=704).  D:  The  average  number  of  positive 
correlations  per  animal  when  both  recordings  were  on 
the  same  side  of  the  cord  (20  ±  3.9  positive 
correlations/animal),  or  on  opposite  sides  of  the  cord  (25 
±  7.5  positive  correlations/animal).  E:  The  number  of 
significant  correlations  expressed  as  a  percent  of  total 
possible  connections  obtained  when  both  recordings 
were  on  the  same  side  of  the  cord  (81/332  total 
possible),  or  on  opposite  sides  of  the  cord  (89/372  total 
possible). 


Mapping  anatomical  locations  and  quantitation  of 
interneuronal  discharge  across  mid-cervical  spinal  laminae 
are  shown  in  Figure  12.  A  representative  example  from 
one  recording  in  which  silver  labeling  was  used  to  “map” 
the  anatomical  locations  of  the  array  electrodes  and 
corresponding  electrophysiological  data  is  shown  in  Figure 
12A-C.  A  Camera  Lucida-style  drawing  of  the  anatomical 
locations  of  eight  silver  labeled  electrodes  was  constructed 
(Figure  12B).  Numerical  identification  of  the  anatomical 
positions  of  each  electrode  corresponds  to  the  schematic 
presented  in  Figure  12D  and  defines  the  rostral-caudal  and 
medial-lateral  positions  of  the  array  electrodes. 


Corresponding  mid-cervical  spinal  discharge 
and  integrated  phrenic  nerve  activity  during 
baseline  and  hypoxia  on  the  left  and  right  hemi- 
cord  is  shown  in  Figure  12C.  Many  these 
recordings  were  from  tonic  firing  interneurons 
(electrodes  3,  8,  9,  11,  12,  and  15)  and  one 
represented  a  phase  switching  (e.g.  from  tonic 
at  baseline  to  expiratory  during  hypoxia) 
interneuron  (electrode  6;  Figure  12C).  Using 
this  technique,  a  summary  of  the  anatomical 


Figure  14.  Anatomical  locations  of  functionally  connected 
interneurons.  Summary  maps  of  silver  labeled  interneurons  with 
significant  peaks  in  cross-correlograms.  The  number  of  positive 
correlations  in  each  lamina  was  normalized  to  the  number  of  neurons 
present  in  the  lamina  and  shaded  accordingly.  A:  When  both  trigger  and 
target  neurons  were  on  the  same  side  of  the  spinal  cord  (e.g.,  unilateral 
recordings),  interneurons  in  dorsal  lamina  (i.e.,  IV,  V  and  VI)  made  and 
received  the  greatest  number  of  excitatory  connections.  B:  When  the 
trigger  and  target  neurons  were  on  opposite  sides  of  the  spinal  cord 
(e.g.,  bilateral  recordings),  interneurons  in  lamina  V,  VIII  and  IX  made 
and  received  the  most  connections. 


locations  of  all  silver  labeled  interneurons  was 
constructed  (Figure  12D).  In  addition,  the  average  discharge  frequency  (Hz)  of  all  silver  labeled  interneurons 
was  presented  in  Figure  12E.  The  results  indicate  a  dorsal-ventral  discharge  gradient,  with  higher  discharge 
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rates  in  lamina  IV  and  X  and  lower  values  in  lamina  VIII  and  IX.  This  discharge  map  of  the  mid-cervical  spinal 
network  demonstrates  the  utility  of  using  the  silver  labeling  technique  to  identify  the  anatomical  locations  of 
ensembles  of  interneurons  recorded  with  a  MEA.  These  results  show  the  utility  of  using  this  technique  to 
standardize  the  sampling  distribution  of  recording  locations  across  experimental  groups. 

Temporal  Relationships  between  Cervical  Interneuron  Discharge.  Anatomical  data  indicate  that  cervical 
interneurons  are  part  of  a  diffuse  and  synaptically  coupled  propriospinal  network,  but  relatively  little  is  known 
about  functional  connectivity  between  cervical  interneurons.  Therefore,  initial  analyses  focused  on  all 
recorded  neurons  in  spinal  intact  rats  (i.e.,  silver  and  non-silver  labeled  neurons)  to  screen  for  short  time 
scale  (i.e.,  0-1  Oms)  features.  Significant  features  were  identified  as  departures  in  the  cross-correlation 
histogram  >  3  standard  deviations  from  the  background  noise.  Using  this  approach,  significant  central  (i.e., 
no  lag  time  from  zero)  and  offset  correlogram  peaks  were  detected  and  shown  in  Figure  13  along  with  a 
summary  of  the  latency  of  significant  features  relative  to  the  trigger  (i.e.,  time  zero). 

Neuronal  pairs  with  significant  correlogram  peaks  that  were  also  silver  labeled  were  used  to  construct 
correlation  summary  maps  to  illustrate  the  anatomical  locations  of  each  neuron  (Figure  14).  The  number  of 
positive  correlations  within  each  lamina  was  normalized  to  the  total  number  of  neurons  present  in  the 
corresponding  lamina.  Unilateral  neuronal  pairs  are  shown  in  Figure  14A  and  bilateral  pairs  are  shown  in 
Figure  14B.  Of  the  positive  unilateral  correlations,  interneurons  in  dorsal  lamina  (i.e.,  IV,  V  and  VI),  made  and 
received  the  greatest  number  of  excitatory  connections  (Figure  10A).  In  contrast,  when  the  trigger  and  target 
neuron  were  on  opposite  sides  of  the  spinal  cord,  interneurons  in  lamina  V,  VIII  and  IX  made  and  received 
the  most  connections  (Figure  14B). 

Conclusions 


To  our  knowledge,  no  prior  study  has  comprehensively  evaluated  the  temporal  characteristics  across 
multiple  mid-cervical  neurons  using  MEA  technology.  The  results  presented  herein  demonstrate  a  previously 
unappreciated  degree  of  connectivity  and  indicate  a  high  prevalence  of  temporally  related  discharge  patterns 
between  mid-cervical  interneurons  with  characteristics  consistent  with  mono-  and  polysynaptic  connections. 
Our  findings  also  provide  a  comprehensive  description  of  a  cost-effective  histological  approach  for  validating 
anatomical  locations  of  MEA  recording  sites  with  a  cost-effective  method  to  identify  the  anatomical  locations 
of  neuronal  ensembles  recorded  with  a  multi-electrode  array  during  acute  preparations  without  the 
requirement  of  specialized  array  electrodes. 

One  potential  application  of  MEA  technology  is  to  examine  how  SCI  alters  the  cervical  spinal  networks. 
Spinal  networks  undergo  substantial  remodeling  following  SCI  and  interneurons  have  been  implicated  in  SCI- 
induced  plasticity  and  motor  recovery.  Our  results  from  a  small  sample  of  spinal-contused  rats  demonstrate 
fundamental  feasibility  of  mapping  changes  in  circuitry;  however,  we  also  identified  potential  caveats  that 
need  further  investigation  to  refine  methods  for  spinal  injury  recordings  of  this  type. 

Collective  Summary  of  Papers  Submitted  or  Published  To-Date 

•  We  have  developed  technical  approaches  to  study  changes  in  phrenic  circuitry  following  cervical  SCI 

which  may  lend  to  future  studies  of  ISMS. 

•  Evaluation  of  activity  recorded  from  silver  labeled  sites  has  also  revealed  a  previously  unappreciated 

degree  of  connectivity  between  mid-cervical  interneurons. 

•  We  have  obtained  critical  proof-of-concept  data  supporting  efficacy  of  ISMS  in  reactivating  phrenic  motor 

circuitry  silenced  by  SCI. 

•  This  finding  has  been  independently  replicated  by  collaborators  at  the  University  of  Washington 

(subcontract). 

•  Our  ISMS  data  to-date  demonstrate  the  potential  for  a  lasting  effect  beyond  the  period  of  stimulation. 

•  Our  ISMS  data  also  very  importantly  demonstrate  the  potential  for  triggering  effective  stimulation  via 

spared  respiratory  activity. 

•  We  performed  the  first  study  to  examine  the  impact  of  high  frequency  epidural  stimulation  (FIF-ES)  on 

phrenic  motor  output  following  incomplete  cervical  spinal  cord  injury. 

•  Short  term  potentiation  of  phrenic  bursting  following  FIF-ES  illustrates  the  potential  for  spinal  stimulation 

to  induce  respiratory  neuroplasticity. 
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•  The  increases  in  tonic  phrenic  output  after  HF-ES  indicate  that  the  continuous  stimulation  paradigm  used 
here  is  unlikely  to  be  useful  for  respiratory  muscle  activation  after  spinal  cord  injury. 

Status  of  Primary  SOW  Tasks: 


Major  Task 

Status 

Major  Task  1 :  To  determine  chronic  changes  in 
neuronal  discharge  patterns  within  the  phrenic  circuit 
region  following  a  C3/4  lateralized  contusion  injury. 

Analysis  continues  to  be  in  progress.  One  publication 
thus  far. 

Major  Task  2:  Terminal  electrophysiological 
comparison  of  epidural  and  ISMS  stimulation  on 
activation  of  the  phrenic  motor  circuit  after  chronic  C2 
hemisection  (C2Hx;  one  month  post-SCI). 

Completed.  On  paper  published,  another  in  revision 

Major  Task  5:  Conduct  Phase  1  of  closed-loop 
intraspinal  stimulation  of  phrenic  motor  neurons 
(PhMNs)  after  C2Hx.  Contralateral  diaphragm  EMG 
triggers  intraspinal  microstimulation  (ISMS)  of  PhMNs 

In  progress  as  part  of  Major  Task  1 1 . 

Major  Task  6:  University  of  Florida  replication  of 
closed-loop/C2Hx  (contralateral  diaphragm  triggered 
stimulus)  results  at  U.W.  and  comparison  with  epidur; 
stimulation. 

Will  be  part  of  Major  Task  1 1 .  Partially  done  via  U.  W 
collaboration.  One  paper  in  preparation. 

Major  Task  7:  Conduct  Phase  II  of  closed-loop 
intraspinal  stimulation  of  phrenic  motor  neurons 
(PhMNs).  Contralateral  intercostal  EMG  triggers 
intraspinal  microstimulation  (ISMS)  of  PhMNs. 

Goal  is  part  of  Major  Task  1 1 . 

Major  Task  8:  Conduct  Open-loop  intraspinal 
stimulation  of  phrenic  motor  neurons  (PhMNs). 
Average  respiratory  activity  recorded  in  Phase  1  &  II 
triggers  ISMS 

Epidural  stimulation  results  have  reduced  the 
immediate  need  for  this  study. 

Major  Task  9:  Conduct  Phase  III  of  closed-loop 
intraspinal  stimulation  of  phrenic  motor  neurons 
(PhMNs).  Repeat  most  promising  closed-loop  methoc 
(Diaphragm  or  intercostal  triggered  ISMS)  using  mid- 
cervical  contusion  model. 

Contusion  injury  will  not  be  done  for  reasons 
discussed.  Otherwise,  goal  is  linked  to  Major  Task  1 1 

Major  Task  1 1 :  Neuroanatomical  studies  on  closed- 
loop  and  open-loop  ISMS  treated  rats. 

In  progress 

3.c.  Plans  for  Accomplishing  Goals: 

As  far  as  Major  Tasks,  Table  1  outlines  the  current  status  of  the  principal  Major  Tasks  proposed  (not 
including  basic  logistical  tasks).  We  have  now  been  approved  for  an  Extension  Without  Funds  which  will  be 
used  to  complete  a  major  objective  of  this  project  -  namely,  extend  our  proof  of  concept  study  of  ISMS  to 
living  rats.  Specifically,  our  goal  is  to  determine  whether:  (a)  chronic  ISMS  delivered  to  unanesthetized  rats 
with  C2Hx  injuries  will  restore  diaphragm  activity  on  the  side  of  injury,  (b)  the  effect  persists  beyond  the 
period  of  stimulation,  and  (c)  we  can  demonstrate  anatomical  neuroplasticity  as  a  result  (i.e.,  Major  Task  11). 
For  that  purpose,  IACUC  and  ACURO  modifications  are  being  prepared  at  this  time,  which  we  hope  will  not 
delay  us  to  any  significant  extent.  Given  costs  of  set-ups  for  this  challenging  experiment,  we  are  hoping  to 
obtain  sufficient  data  on  which  to  build  a  future  research  application  to  support  a  larger  scale  study,  while  at 
the  same  time  essentially  meeting  the  original  goals  of  this  project.  Collaborators  at  the  Univ.  Washington  are 
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also  ramping  up  to  do  chronic  ISMS  studies  as  part  of  the  independent  replication  we  proposed.  Chronic 
diaphragm  EMG  recordings  are  now  being  worked  out  and  those  signals  to  trigger  ISMS  will  be  tested  after 
C2Hx. 


4.  Impact: 

4. a.  Impact  on  Principal  Disciplines: 

The  results  we  have  obtained  to-date  provide  another  dimension  to  the  potential  benefits  of  ISMS  on 
motor  function.  Up  to  now,  much  of  what  we  know  about  the  potential  for  evoking  function  from  the  injured 
spinal  cord  with  electrical  stimulation  has  derived  from  studies  of  epidural  stimulation  or  ISMS  targeting  other 
somatic  motor  systems.  Our  studies  are  now  moving  incrementally,  but  significantly,  towards  a  better 
understanding  of  how  ISMS,  despite  its  more  invasive  nature,  may  not  only  restore  function  but  also  how  it 
may  promote  plasticity  such  that  functional  recovery  can  be  achieved  beyond  any  need  for  stimulation  via 
other  sources.  Electrical  modulation  of  spinal  circuits  is  becoming  an  area  of  expanding  interest  (Samaddar 
et  al.  2016;  Song  and  Martin  2016;  Yao  and  Li  2016). 

4.b.  Impact  on  Other  Disciplines 

Nothing  to  report 

4.c.  Impact  on  Technology  Transfer 

Nothing  to  report 

4.d.  Societal  Impact 
Nothing  to  report 

5.  Changes/Problems 

5a.  Changes  in  Approach 

We  originally  proposed  to  carry  out  our  studies  in  two  different  injury  models  -  after  C2  hemisection 
(C2Hx)  or  a  C3-4  lateralized  contusion.  Given  the  technical  challenges  encountered,  the  decision  was  made 
to  forgo  the  contusion  injury  model  for  two  reasons.  First,  there  is  still  controversy  in  the  literature  related  to 
the  type  of  deficit  that  occurs  following  a  mid-cervical  contusion  (Alvarez-Argote  et  al.  2015;  Rana  et  al. 
2016),  and  it  seems  that  significant  respiratory  deficits  may  involve  more  severe  and  debilitating  injuries  than 
logistically  feasible  or  ethical  in  experimental  animals  at  this  time.  At  the  very  least,  additional  experiments 
that  were  not  budgeted  would  be  required  to  support  any  effect  of  spinal  electrical  stimulation.  Second,  the 
C2Hx  model  is  often  the  lesion  of  choice  for  studies  of  cSCIs  and  respiratory  outcomes.  Therefore,  we  feel 
that  all  objectives  of  the  original  proposal  can  be  adequately  addressed  with  that  injury  model  as  it  results  in 
a  more  complete  inactivation  of  the  phrenic  circuit  on  the  side  of  injury  than  even  a  lateralized  contusion  can 
achieve.  In  sum,  our  revised  plan  is  not  a  departure  from  the  original  SOW  but  instead  a  reorganization  of 
experimental  priorities  based  upon  experiences  to  date.  These  changes  have  been  indicated  in  previous 
Quarterly  Reports. 

5b.  Anticipated  Problems  and  Proposed  Corrective  Action 

None 

5c.  Changes  Affecting  Expenditures 

None 

5d.  Changes  Affecting  Use  and  Care  of  Animals 
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None  for  this  reporting  period,  but  IACUC  and  ACURO  modifications  will  be  submitted  for  approval  for 
remaining  work  under  the  EWOF. 
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Project  Role 
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Month  Workec 
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Contribution  to  Project 

University  of  Florida 

Paul  J.  Reier,  Ph.D. 

Project  Director 

No  change 

Oversees  entire  project;  conducts 
experiments  related  to  Major  Task  3 

David  D.  Fuller,  Ph.D. 

Co-1 

No  change 

Functions  as  Associate  Project  Director 
responsible  for  overseeing 
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with  Seattle  researchers 

David  Baekey,  Ph.D. 

Co-1 

No  change 
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Task  1 
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Post-Doctoral 

Fellow 

4.0 

Conducted  recordings  for  Major  Task  1 

Elisa  Gonzalez-Rothi,  Ph.D. 

Post-Doctoral 

Fellow 

2.0 

Epidural  stimulation  studies  for  Major 
Task 

2. 

Lynne  Mercier2 

Graduate  Student 

No  change 

Conducted  ISMS  proof-of-concept  stud) 
for 

Major  Task  2. 

Lucy  Denholtz 

Technician 

No  change 

Assists  with  surgeries  and  histological 
procedures. 

Savannah  Posgai3 

Technician 

Left  program 
to  go  to 
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school. 

Assists  with  surgeries  and  histological 
procedures. 

Joey  Tringali 

Technician 

0.5 

Assists  with  surgeries  and  histological 
procedures.  Replaced  Ms.  Posgai. 

University  of  Washington 

Chet  Moritz,  Ph.D. 

Sub-contract 
project  director 

1.73 

Oversees  project  Major  Tasks 
to  be  done  at  U.  Wash. 

Michael  Sunshine 

Graduate  Student 

12.0 

Primary  person  responsible  for 
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carrying  out  U.  Washington  studies.  He 
has  now  come  to  the  Univ.  Florida  to 
pursue  graduate  studies  and  will  assist 
with  the  completion  of  this  project. 
Financial  support  derived  through  other 
sources. 

7b.  Changes  in  Active  Support/Personnel 
No  changes  during  this  reporting  period. 

8.  Special  Reporting  Requirements 
8a.  Quad  Chart  (Appendix  Attachment) 


9.  Appendices 

See  attached  PDF  files. 
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Mercier  LM,  Gonzalez-Rothi  EJ,  Streeter  KA,  Posgai  SS, 
Poirier  AS,  Fuller  DD,  Reier  PJ,  Baekey  DM.  Intraspinal  micro¬ 
stimulation  and  diaphragm  activation  after  cervical  spinal  cord  injury. 
J  Neurophysiol  117:  000-000,  2017.  First  published  November  23, 
2016;  doi:10.1 152/jn.00721.2016. — Intraspinal  microstimulation  (ISMS) 
using  implanted  electrodes  can  evoke  locomotor  movements  after 
spinal  cord  injury  (SCI)  but  has  not  been  explored  in  the  context  of 
respiratory  motor  output.  An  advantage  over  epidural  and  direct 
muscle  stimulation  is  the  potential  of  ISMS  to  selectively  stimulate 
components  of  the  spinal  respiratory  network.  The  present  study 
tested  the  hypothesis  that  medullary  respiratory  activity  could  be  used 
to  trigger  midcervical  ISMS  and  diaphragm  motor  unit  activation  in 
rats  with  cervical  SCI.  Studies  were  conducted  after  acute  (hours)  and 
subacute  (5-21  days)  C2  hemisection  (C2Hx)  injury  in  adult  rats. 
Inspiratory  bursting  in  the  genioglossus  (tongue)  muscle  was  used  to 
trigger  a  250-ms  train  stimulus  (100  Hz,  100-200  pA)  to  the  ventral 
C4  spinal  cord,  targeting  the  phrenic  motor  nucleus.  After  both  acute 
and  subacute  injury,  genioglossus  EMG  activity  effectively  triggered 
ISMS  and  activated  diaphragm  motor  units  during  the  inspiratory 
phase.  The  ISMS  paradigm  also  evoked  short-term  potentiation  of 
spontaneous  inspiratory  activity  in  the  previously  paralyzed  hemidi- 
aphragm  (i.e.,  bursting  persisting  beyond  the  stimulus  period)  in 
~70%  of  the  C2Hx  animals.  We  conclude  that  medullary  inspiratory 
output  can  be  used  to  trigger  cervical  ISMS  and  diaphragm  activity 
after  SCI.  Further  refinement  of  this  method  may  enable  “closed-loop- 
like”  ISMS  approaches  to  sustain  ventilation  after  severe  SCI. 

NEW  &  NOTEWORTHY  We  examined  the  feasibility  of  using 
intraspinal  microstimulation  (ISMS)  of  the  cervical  spinal  cord  to 
evoke  diaphragm  activity  ipsilateral  to  acute  and  subacute  hemisec¬ 
tion  of  the  upper  cervical  spinal  cord  of  the  rat.  This  proof-of-concept 
study  demonstrated  the  efficacy  of  diaphragm  activation,  using  an 
upper  airway  respiratory  EMG  signal  to  trigger  ISMS  at  the  level  of 
the  ipsilesional  phrenic  nucleus  during  acute  and  advanced  postinjury 
intervals. 

phrenic  motor  nucleus;  rat;  respiration;  hypoglossal  respiratory  activ¬ 
ity;  diaphragm  function 


severe  respiratory  compromise  often  occurs  after  spinal  cord 
injury  (SCI)  at  upper  cervical  to  midcervical  spinal  levels 
(Mansel  and  Norman  1990;  Winslow  and  Rozovsky  2003). 
When  independent  breathing  is  not  possible,  positive-pressure 
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of  Physiological  Sciences,  PO  Box  100144,  Gainesville,  FL  32610-0144 
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mechanical  ventilation  is  often  used  to  sustain  alveolar  venti¬ 
lation,  but  at  the  risk  of  potential  rapid  diaphragm  atrophy, 
atelectasis,  and  respiratory  tract  infections  (Bezzant  and 
Mortensen  1994;  Laghi  et  al.  2003;  Smuder  et  al.  2016). 
Diaphragm  and  phrenic  nerve  pacing  options  are  clinically 
available  (DiMarco  et  al.  2005;  Glenn  and  Phelps  1985; 
Onders  et  al.  2009),  and  recent  preclinical  work  has  begun 
exploring  high-frequency  epidural  stimulation  of  the  spinal 
cord  to  activate  the  diaphragm  (DiMarco  and  Kowalski  2013, 
2015;  Kowalski  et  al.  2013). 

Another  approach  that  may  effectively  activate  respiratory 
muscles  after  SCI  is  intraspinal  microstimulation  (ISMS).  Im¬ 
plementation  of  this  method  in  motor  systems  has  been  dem¬ 
onstrated  (Giszter  2015;  Kasten  et  al.  2013;  McPherson  et  al. 
2015;  Tator  et  al.  2012),  but  applications  to  respiratory  motor 
output  after  cervical  SCI  have  not  been  investigated.  Potential 
advantages  of  ISMS  include  the  ability  to  deliver  currents  over 
a  relatively  wide  range  of  intensities  and  with  high  selectivity 
for  specific  motor  systems  (Mondello  et  al.  2015;  Mushahwar 
et  al.  2000;  Sunshine  et  al.  2013). 

The  present  study  was  therefore  designed  to  obtain  proof- 
of-concept  evidence  in  support  of  cervical  ISMS  as  a  means  of 
activating  diaphragm  motor  units  after  cervical  SCI.  To  move 
toward  a  “closed-loop-like”  design  for  ISMS,  we  recorded  the 
inspiratory  output  of  an  upper  airway  muscle  (the  genioglos¬ 
sus,  innervated  by  cranial  nerve  XII)  and  used  this  to  trigger 
stimulation.  Thus  cervical  ISMS  was  activated  by  endogenous 
inspiratory  drive.  The  data  show  that  a  physiologically  relevant 
medullary  inspiratory  output  can  be  used  to  trigger  cervical 
ISMS  and  that  this  approach  can  effectively  activate  diaphragm 
motor  units  after  acute  and  subacute  SCI. 


MATERIALS  AND  METHODS 

Thirteen  female  Sprague-Dawley  rats  (255  ±  16  g;  Harlan,  India¬ 
napolis,  IN)  were  distributed  between  two  experimental  groups:  1) 
ISMS  after  acute  SCI,  tested  before  and  immediately  after  lateral 
hemisection  of  the  spinal  cord  at  C2  (C2Hx)  (n  =  8),  and  2)  ISMS 
after  subacute  SCI,  tested  5-21  days  after  C2Hx  (n  =  5).  All  data 
were  collected  in  terminal  procedures  using  protocols  approved  by  the 
Institutional  Animal  Care  and  Use  Committee  at  the  University  of 
Florida  and  the  US  Army  Medical  Research  and  Materiel  Command 
Animal  Care  and  Use  Review  Office. 
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C2Hx  surgery.  For  survival  surgery,  rats  were  initially  anesthetized 
with  isoflurane  (3-5%  in  02)  in  a  closed  chamber  and  then  anesthesia 
was  maintained  via  nosecone  (1.5-3%  in  02).  A  dorsal  incision  was 
made  over  the  cervical  spine,  followed  by  a  C2  laminectomy  and 
durotomy.  For  acute  C2Flx,  resection  lesions  were  performed  with 
a  needle-blade  microknife  (Fine  Science  Tools,  Foster  City,  CA). 
Subacute  C2Hx  resection  lesions  were  made  with  microscissors 
and  fine  forceps  followed  by  gentle  aspiration  of  tissue  to  complete 
the  hemisection  (Fuller  et  al.  2008).  The  dura  and  overlying 
muscles  were  closed  and  lactated  Ringer  solution  (5  ml  sc)  and 
buprenorphine  were  administered  (0.03  mg/kg  sc.;  Flospira,  Lake 
Forest,  IL).  Postoperative  care  consisted  of  daily  lactated  Ringer 
solution  (5  ml/day  sc)  and  oral  nutritional  supplement  (1-3  ml/day; 
Nutrical,  Webster  Veterinary).  Buprenorphine  was  given  at  ~12-h 
intervals  for  2  days  after  surgery. 

General  neurophysiology  protocols.  Rats  were  anesthetized  as 
above.  Rectal  temperature  was  maintained  at  ~37.5°C  by  a  heating 
pad  (CWE,  Ardmore,  PA).  A  femoral  artery  was  catheterized  (PE-50) 
for  blood  pressure  measurements  (Statham  P-10EZ  pressure  trans¬ 
ducer,  CP122  AC/DC  strain  gauge  amplifier;  Grass  Instruments,  West 
Warwick,  RI)  and  arterial  blood  samples.  The  femoral  vein  of  the 
same  hindlimb  was  also  catheterized  (PE-50)  for  supplemental  fluid 
administration  and  conversion  from  isoflurane  to  urethane  anesthesia 
(1.7  g/kg  iv;  Sigma,  St.  Louis,  MO).  Animals  received  a  tracheotomy 
and  were  mechanically  ventilated  (50-65%  02,  balance  N2;  6-7 
ml/kg  volume;  70-72  bpm  frequency;  Harvard  Apparatus,  Holliston, 
MA)  throughout  the  experimental  procedures.  Since  spontaneously 
breathing  anesthetized  rats  can  rapidly  become  hypercapnic,  we 
elected  to  employ  mechanical  ventilation  for  these  initial  experiments 
in  order  to  keep  blood  gases  stable  and  remove  Pa02  and/or  PaCOi 
fluctuations  as  confounding  variables. 

End-tidal  C02  (Capnogard  1265;  Respironics,  Wallingford,  CT) 
was  continuously  monitored,  and  arterial  blood  gases  (iSTATl;  Ab¬ 
bott,  Princeton,  NJ)  were  periodically  assessed  from  0.1 -ml  arterial 
samples.  On  the  basis  of  these  measures,  the  inspired  C02  content  was 
adjusted  to  maintain  Pacc,2  at  40  mmHg.  If  base  excess  was  greater 
than  —3  meq,  this  was  corrected  with  intravenous  administration  of 
sodium  bicarbonate  solution  [8.4%,  Hospira;  dose  (ml)  =  0.3-weight 
(kg)-standard  base  excess]. 

Recordings  of  EMG  activity  from  respiratory-related  muscles  (i.e., 
diaphragm,  genioglossus,  and  intercostal)  and  an  off-target,  nonrespi- 
ratory  muscle  [i.e.,  extensor  carpi  radialis  longus  (ECR)]  were  ob¬ 
tained  with  pairs  of  Teflon-coated  tungsten  hooked  wires  (A-M 
Systems,  Sequim,  WA).  Recordings  of  the  genioglossus  muscle  were 
obtained  at  the  base  of  the  tongue  (Fuller  et  al.  1999;  Fuller  and 
Fregosi  2000)  and  exhibited  a  centrally  driven  inspiratory  rhythm  that 
was  used  to  trigger  ISMS.  For  intercostal  EMG  recordings,  a  small 
incision  lateral  to  the  sternum  was  made  at  the  T2  level  and  wires  were 
placed  1-1.5  cm  lateral  to  the  midline. 

The  spinal  cord  was  exposed  via  a  cervical  middorsal  incision 
followed  by  laminectomy  and  durotomy  from  C2  to  C5.  Raw  EMG 
signals  were  amplified  at  100-10K,  band-pass  filtered  at  100  Hz-10 
kHz  (A-M  Systems,  Carlsborg,  WA),  and,  in  the  case  of  the  genio¬ 
glossus  EMG  recording,  passed  through  a  moving  time  averager  (50 
ms  time  constant;  CWE,  Ardmore,  PA).  The  moving  time  average 
signal  was  used  to  trigger  the  stimulator  for  ISMS.  All  EMG  signals 
were  digitized  at  25  kHz  (CED  Power  1401)  and  recorded  (Spike2  v8; 
CED,  Cambridge,  UK)  to  a  PC  and  then  analyzed  off-line. 

ISMS  at  midcervical  spinal  cord.  A  tungsten  microwire  electrode 
(FHC,  Bowdoin,  ME)  was  mounted  in  a  stereotaxic  micromanipulator 
(David  Kopf  Instruments,  Tujunga,  CA).  The  microwire  had  a 
100-p.m  segment  at  the  tip  stripped  of  all  insulation.  The  electrode 
was  placed  above  the  C4  segment  with  the  dorsal  root  entry  zone  as 
the  lateral  anatomical  landmark.  Electrical  activity  was  initially  re¬ 
corded  via  the  stimulating  electrode  and  assessed  visually  with  Spike2 
software  and  audibly  with  an  AM8  Audio  Monitor  (Grass  Technolo¬ 
gies,  Quincy,  MA).  The  software  program  simultaneously  displayed 


the  inspiratory  diaphragm  EMG  recording;  these  procedures  enabled 
placement  of  the  stimulating  electrode  tip  in  proximity  to  inspiratory 
neurons  (i.e.,  targeting  phrenic  motoneurons).  Inspiratory  bursting 
was  absent  in  rats  with  subacute  C2Hx  injury.  Therefore,  a  stimulator 
and  constant-current  stimulus  isolation  unit  (S88X  and  SIU-C;  Grass 
Technologies,  Warwick,  RI)  were  used  to  deliver  single  pulses 
(0.3  ms  duration)  during  electrode  descent  with  gradually  decreasing 
currents  (200,  100,  50  pA).  This  approach  allowed  us  to  determine  a 
location  for  ISMS  that  evoked  left  diaphragm  activation. 

Stimulation  of  the  spinal  cord  was  initiated  (triggered)  with  the 
EMG  signal  recorded  from  the  genioglossus  muscle.  Specifically,  a 
“threshold  crossing”  was  established  such  that  when  the  inspiratory 
integrated  EMG  burst  (50  ms  time  constant)  reached  a  preset  ampli¬ 
tude  the  spinal  cord  was  stimulated.  All  stimulations  during  the 
inspiratory  phase  were  made  at  the  onset  of  the  genioglossus  inspira¬ 
tory  burst,  while  stimulations  made  during  expiratory  periods  were 
accomplished  by  adding  a  400-ms  delay  to  the  trigger. 

Experimental  group  1:  ISMS  before  and  after  acute  C2Hx.  The 
genioglossus  EMG  signal  was  used  to  trigger  ISMS  during  two 
consecutive  inspiratory  cycles,  followed  by  ISMS  delivered  during 
two  consecutive  expiratory  cycles.  In  both  cases,  repeated  250-ms 
stimulus  trains  (200  pA)  were  delivered  with  0.3  ms  pulse  duration 
and  100  Hz  stimulus  frequency.  Once  these  initial  stimulations  were 
complete,  ISMS  was  delivered  during  the  inspiratory  cycle  for  1  min. 
All  of  the  aforementioned  ISMS  was  done  with  the  spinal  cord  intact. 
After  C2Hx,  the  ISMS  protocol  was  repeated,  as  described  above. 
Subsequent  to  protocols  performed  with  spinal  cord  intact  or  after 
C2Hx,  an  additional  bout  of  ISMS  was  administered  after  neuromus¬ 
cular  blockade  via  intravenous  pancuronium  bromide  (2.5  mg/kg; 
Hospira)  to  confirm  unequivocally  that  the  evoked  activity  was  not 
contaminated  by  a  stimulus  artifact. 

Selection  of  stimulus  parameters.  Pilot  experiments  were  done  with 
a  manually  triggered,  open-loop  approach  in  spinal-intact  animals. 
Repeated  250-ms  trains  of  ISMS  were  targeted  to  the  left  phrenic 
motor  nucleus,  and  EMG  activity  was  recorded  in  both  hemidi- 
aphragms,  the  genioglossus,  and  ipsilateral  intercostal  and  forelimb 
(ECR)  muscles  (see,  e.g..  Fig.  1).  A  range  of  stimulation  frequencies 
(50,  100,  200,  300  Hz)  and  currents  (50-200  pA)  were  tested  to 
optimize  parameters  for  eliciting  compound  motor  unit  action  poten¬ 
tials  (MUAPs)  in  the  left  diaphragm.  Activation  of  the  ipsilateral 
hemidiaphragm  was  more  pronounced  when  100-Hz  stimulation  was 
used,  with  stimulation  frequencies  above  200  Hz  thereby  producing 
large  contractior^of  nonrespiratory  muscles.  Preliminary  experiments 
also  indicated  that  progressive  attenuation  of  ISMS-induced  phrenic 
MUAP  amplitude  occurred  during  sustained  stimulation  (e.g.,  1  min) 
at  200  Hz  but  not  at  lower  stimulus  frequencies  (50  Hz,  100  Hz).  In 
the  first  series  of  experiments,  a  stimulus  current  of  200  pA  was  used 
to  activate  the  ipsilateral  diaphragm. 

Experimental  group  2:  ISMS  after  subacute  C2Hx.  The  intraspinal 
electrode  placement  was  determined  as  described  above.  Triggered 
ISMS  was  then  delivered  during  two  consecutive  expiratory  periods, 
followed  by  one  continuous  minute  of  stimulation  during  the  inspira¬ 
tory  phase.  In  both  cases,  a  continuous  250-ms  train  was  delivered 
with  0.3-ms  pulses  at  100  Hz  and  100  pA.  Stimulations  also  were 
repeated  after  neuromuscular  blockade  as  described  above. 

Consistent  with  other  reports,  we  used  cessation  of  the  ipsilateral 
diaphragm  inspiratory  EMG  burst  and  histology  as  functional  and 
anatomical  verification,  respectively,  of  the  subacute  C2Hx  lesions 
(Goshgarian  1981).  After  the  electrophysiology  protocols,  the  animals 
were  perfused  with  saline  followed  by  4%  paraformaldehyde  (Sigma). 
The  cervical  spinal  cords  from  the  subacute  C2Hx  animals  were 
subsequently  harvested,  and  a  tissue  block  including  the  C2  region 
was  paraffin  embedded,  sectioned  (8  pm),  and  counterstained  with 
cresyl  violet. 

Biomechanical  impact  of  ISMS.  Changes  in  tracheal  pressure  were 
assessed  as  a  crude  indicator  of  the  biomechanical  impact  of  ISMS.  It 
should  be  noted,  however,  that  this  was  not  considered  to  be  a  primary 
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outcome  variable  but  rather  an  accessory  measurement  intended  to 
provide  some  insight  regarding  the  functional  impact  of  spinal  stim¬ 
ulation.  The  pressures  associated  with  lung  inflation  and  deflation 
were  measured  at  the  tracheal  cannula.  Measurements  were  taken 
continuously  and  used  to  compare  baseline  conditions  to  stimulation 
periods.  For  each  period  of  interest,  tracheal  pressure  deflections  for 
20  consecutive  breaths  were  averaged  and  presented  in  a  box  and 
whisker  plot  (SPSS;  IBM,  Armonk,  NY).  Passive  tracheal  pressure 
values  during  exhalation  were  obtained  in  each  preparation  after 
neuromuscular  blockade  and  thus  provided  a  baseline  measurement 
that  was  devoid  of  any  respiratory  contribution  from  the  animal  (for 
additional  details,  see  results). 

Data  acquisition  and  analyses.  For  the  first  group  of  experiments, 
results  from  the  spinal-intact  condition  were  compared  to  data 
collected  10  min  after  acute  C2Flx.  In  addition  to  the  evoked 
responses  (i.e.,  response  during  ISMS),  we  also  evaluated  the 
impact  of  the  ISMS  on  spontaneous  EMG  activity  (i.e.,  following 
the  period  of  stimulation).  To  assess  ISMS-entrained  activation  of 
motor  units,  the  EMG  signals  were  averaged  with  respect  to 
stimulus  pulses  within  each  trial.  The  resulting  stimulus-triggered 
average  (McPherson  et  al.  2015;  Moritz  et  al.  2007)  served  to 
minimize  activity  unrelated  to  ISMS.  The  resultant  waveform 
averages  represented  ISMS-entrained  MUAPs.  All  data  were  col¬ 


lected  with  a  CED  Spike2  data  acquisition  system  and  subse¬ 
quently  analyzed  with  Spike2  v8  software  on  a  standard  PC. 
Values  are  reported  as  means  ±  SD. 

Statistics.  A  one-way  repeated-measures  ANOVA  followed  by  the 
Tukey  post  hoc  test  (SigmaPlot;  Systat  Software,  San  lose,  CA)  was 
used  to  compare  data  collected  during  baseline  vs.  ISMS.  Dependent 
variables  included  the  amplitude  of  MUAPs  and  spontaneous  EMG 
burst  amplitude,  blood  pressure,  tracheal  pressure,  heart  rate,  and 
blood  gas  data.  A  “detectability  index”  statistical  test  (Aertsen  and 
Gerstein  1985)  was  used  to  compare  ISMS-evoked  MUAP  amplitude 
to  background  activity  in  the  EMG  recordings.  We  used  the  modified 
version  of  this  test  with  a  more  stringent  value  of  D  a  3  to  prevent 
“false  positives”  (Meissen  and  Epping  1987). 

RESULTS 

ISMS  with  spinal  cord  intact.  A  representative  example  of 
triggered  cervical  ISMS  during  the  inspiratory  phase  is  pro¬ 
vided  in  Fig.  1.  Note  that  compound  MUAPs  in  the  left  Ft 
diaphragm  are  clearly  discernible  (Fig.  1C)  and  are  entrained  to 
each  stimulation  pulse  (e.g.,  Fig.  ID).  In  this  example,  the 
average  latency  between  stimulus  artifact  and  MUAP  peak  is 
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Fig.  1.  Cervical  ISMS  activates  the  ipsilateral  diaphragm  A:  schematic  diagram  illustrates  the  EMG  recording  sites  relative  to  the  placement  of  the  stimulating 
electrode  (S).  The  stimulating  electrode  was  placed  in  the  immediate  vicinity  of  the  phrenic  motor  nucleus,  and  ISMS  was  initiated  via  a  trigger  signal  based 
on  inspiratory  tongue  (genioglossus)  EMG  activity.  LI,  left  intercostal.  B:  pre-C2Hx  representative  EMG  recordings  from  the  genioglossus,  extensor  carpi  radialis 
(ECR),  and  both  sides  of  the  diaphragm  showing  baseline  (prestimulation)  activity  and  a  period  of  genioglossus-triggered  cervical  ISMS  (100  Hz,  200  pA;  0.3-ms 
pulse  duration,  250-ms  train  duration)  represented  by  a  gray  bar  in  the  / tongue  trace.  C:  expanded  trace  of  the  EMG  recording  from  the  left  diaphragm  during 
the  period  of  ISMS  shows  each  stimulus  artifact,  indicated  by  asterisks,  and  each  subsequent  MUAP.  D:  overlay  of  each  elicited  MUAP,  aligned  by  the  stimulus 
artifact,  demonstrates  constant  latency  and  amplitude.  I 
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2.14  ±  0.05  ms,  which  is  consistent  with  a  synaptic  delay 
before  motor  unit  activation.  An  example  of  ISMS  delivered 
during  the  expiratory  phase  is  provided  in  Fig.  2.  In  the  absence 
of  spontaneous  inspiratory  EMG  activity,  it  can  be  appreciated 
that  ISMS  targeting  the  left  phrenic  motor  nucleus  evoked  a 
marked  response  in  both  the  ECR  and  diaphragm  but  had 
minimal  impact  on  the  left  tongue  or  right  diaphragm.  To 
confirm  that  the  very  small-amplitude  EMG  signals  in  the 
tongue,  diaphragm,  and  intercostals  represented  activation  of 
motor  units,  we  employed  a  previously  published  “detectability 
index”  for  evoked  potentials  (Aertsen  and  Gerstein  1985; 
Meissen  and  Epping  1987).  This  assessment  indicated  that 
ISMS  evoked  a  small  but  statistically  significant  increase  in 
entrained  and  averaged  EMG  activity  in  the  tongue  (5  of  8 
animals),  left  intercostal  (8  of  8  animals),  and  right  diaphragm 
(8  of  8  animals).  It  should  be  noted,  however,  that  the  ampli¬ 
tudes  of  these  off-target  evoked  potentials  were  modest  com¬ 
pared  with  those  of  the  left  diaphragm  and  ECR  muscles. 
Neuromuscular  blockade  eliminated  evoked  potentials,  as  ex¬ 
pected  (Fig.  2C). 

ISMS  after  high  cervical  SCI.  Experiments  were  performed 
immediately  after  C2Hx  (acute  injury)  and  also  in  rats  that 
were  5-21  days  post-C2Hx  (subacute  injury).  These  tempo¬ 


rally  advanced  lesion  experiments  were  performed  to  deter¬ 
mine  whether  degenerative  processes  triggered  by  C2Hx  (e.g., 
axonal  retraction)  prevented  or  mitigated  the  impact  of  ISMS 
on  diaphragm  activation.  The  general  experimental  paradigm  is 
illustrated  in  Fig.  3A. 

In  the  acute  C2Hx  experiments,  diaphragm  MUAPs  were  first 
evaluated  with  the  spinal  cord  intact  (Fig.  3 B).  The  ISMS-evoked 
potentials  in  the  spinal-intact  condition  were  similar  to  those 
described  in  the  preceding  section.  Acute  C2Hx  abolished 
spontaneous  inspiratory  EMG  activity  in  the  ipsilateral  hemi- 
diaphragm  in  all  animals  (Fig.  3 C).  However,  cervical  ISMS 
still  produced  clear  MUAPs  in  the  ipsilateral  diaphragm  after 
the  acute  lesion  (Fig.  3D).  Before  C2Hx,  the  average  latency 
from  stimulus  to  the  peak  of  the  MUAP  was  1.99  ±  0.25  ms. 
After  C2Hx  the  value  tended  to  be  reduced  (1.69  ±  0.57  ms), 
but  this  did  not  approach  statistical  significance  (P  = 
0.190).  In  seven  of  eight  experiments,  the  amplitude  of  the 
ISMS-evoked  MUAP  remained  relatively  consistent  after 
the  acute  C2Hx,  ranging  from  80%  to  124%  of  the  preinjury 
value,  with  a  mean  of  94  ±  15%  (P  =  0.065  vs.  pre-C2Hx). 
There  was  one  significant  outlier  rat  (confirmed  with 
Grubbs’  outlier  test,  P  <  0.05)  that  demonstrated  a  profound 
reduction  in  MUAP  amplitude  after  acute  C2Hx,  with 
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B 


Fig.  2.  Examples  of  raw  EMG  activity  recorded 
in  the  tongue  (genioglossus),  forelimb  (ECR), 
and  the  left  and  right  hemidiaphragm  during 
ISMS.  A:  representative  examples  before 
C2Hx.  B:  stimulus-triggered  averages  of  EMG 
activity  show  prominent  MUAPs  in  the  fore¬ 
limb  and  ipsilateral  diaphragm.  C:  stimulus- 
triggered  averages  of  EMG  activity  after  neu¬ 
romuscular  blockade.  In  B  and  C,  the  stimulus 
artifact  is  highlighted  by  the  gray  boxes.  Stim¬ 
ulus-triggered  averages  were  scaled  to  the  same 
values  and  represent  25  stimulus  triggers.  In 
these  examples,  ISMS  was  delivered  at  100  Hz 
and  200  pA  during  expiration. 
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Fig.  3.  ISMS  after  acute  C2Hx.  A:  illustration  of  the  electrode  placement  relative  to  the  C2Hx  lesion.  The  solid  dark  bar  represents  the  injury  site;  the  dashed 
lines  represent  severed  pathways.  B\  examples  of  EMG  activity  in  the  right  and  left  diaphragm  before  C2Hx.  C:  examples  of  left  and  right  diaphragm  EMG 
activity  after  acute  C2Hx.  Note  that  activity  is  abolished  in  the  left  diaphragm.  In  B  and  C,  3  spontaneous  breaths  are  shown  followed  by  3  breaths  during  triggered 
cervical  ISMS  (solid  bars;  100  Hz,  200  /j,A).  D:  stimulus-triggered  averages  of  right  and  left  diaphragm  EMG  activity  obtained  from  the  period  of  ISMS 
highlighted  by  dashed  boxes  in  B  and  C.  Stimulus-triggered  averages  were  scaled  to  the  same  values  and  represent  25  stimulus  triggers.  In  D,  the  stimulus  artifact 
AQ:  6  is  highlighted  by  the  gray  boxes. 


evoked  responses  reaching  only  8%  of  the  preinjury  ampli¬ 
tude.  If  that  particular  data  point  is  included  in  the  calcula¬ 
tion  of  the  overall  mean  the  values  after  acute  C2Hx  are 
83  ±  35%  of  the  preinjury  condition,  but  the  added  variance 
actually  produces  a  higher  P  value  ( P  =  0.139  vs.  prein¬ 
jury).  Collectively,  the  data  are  consistent  with  a  slight 
reduction  in  cervical  ISMS-evoked  MUAP  amplitude  in 
most  animals  after  acute  C2Hx.  The  off-target  impact  of 
ISMS  after  C2Hx  (i.e.,  evoked  activity  in  ECR,  tongue,  and 
right  diaphragm)  was  indistinguishable  from  that  reported 
above  for  the  spinal-intact  condition  (data  not  shown). 

Tracheal  pressures  were  evaluated  to  explore  the  potential  of 
a  biomechanical  impact  of  ISMS-induced  muscle  contraction 
F4  (Fig.  4,  A  and  B).  Negative  pressure  swings  in  tracheal  pressure 
would  be  expected  in  this  preparation  if  diaphragm  contraction 
is  altering  the  dimensions  of  the  thoracic  cavity.  In  the  spinal- 
intact  condition,  genioglossus-triggered  cervical  ISMS  caused 
a  significant  negative  deflection  in  tracheal  pressure  ( P  =  0.02 
vs.  baseline;  Fig.  4Cj,  as  expected.  The  acute  C2Hx  injury 
resulted  in  a  small,  but  statistically  significant,  change  in  the 
tracheal  pressures  that  were  recorded  during  ventilator-induced 
lung  inflation  (i.e.,  independent  of  ISMS;  P  =  0.028;  Fig.  4 C). 
After  acute  C2Hx,  ISMS  also  induced  a  change  in  tracheal 
pressure  (P  =  0.005  vs.  baseline),  thus  suggesting  a  biome¬ 
chanical  impact  of  the  stimulation.  However,  the  relative  mag¬ 
nitude  of  ISMS-induced  changes  in  tracheal  pressure  were 
attenuated  after  C2Hx  (P  =  0.038;  Fig.  4 Cj. 


Spontaneous  ipsilateral  diaphragm  EMG  activity  was  ab¬ 
sent  in  rats  studied  5-21  days  after  C2Hx  (Fig.  5A),  and  fs 
subsequent  histological  evaluation  of  the  spinal  cord  indi¬ 
cated  anatomically  complete  hemilesion  in  all  animals  (Fig. 

5 B).  The  ISMS  procedure  evoked  clearly  discernible 
MUAPs  in  the  ipsilateral  (paralyzed)  hemidiaphragm  (Fig. 
5C).  Compared  with  the  acute  injury  group,  the  only  appar¬ 
ent  difference  in  the  ISMS  responses  in  the  animals  with  the 
subacute  lesions  was  a  trend  for  more  variable  and  longer 
latencies.  On  average,  however,  there  were  no  statistical 
differences  between  the  groups  for  latency  (acute  C2Hx: 
1.69  ±  0.57  ms,  subacute  C2Hx:  2.50  ±  0.67  ms),  but  the 
increase  in  latency  in  the  subacute  group  was  close  to 
threshold  for  significance  (P  =  0.065).  The  MUAP  ampli¬ 
tude  was  variable,  but  with  no  evidence  for  a  difference 
between  the  two  groups  (acute  C2Hx:  1.62  ±  1.45  mV, 
subacute  C2Hx:  1.76  ±  2.65  mV,  P  =  0.833).  ISMS  caused 
the  expected  negative  deflection  in  tracheal  pressure  in  the 
subacute  C2Hx  animals,  thus  confirming  that  there  was  a 
biomechanical  impact  of  the  stimulation  (Fig.  5 D).  The 
off-target  muscle  activation  was  similar  to  what  was  ob¬ 
tained  in  the  spinal-intact  and  acute  C2Hx  animals.  Thus 
very  small  but  statistically  significant  changes  in  EMG 
activity  during  ISMS  were  detected  in  the  tongue  (4/5 
animals),  intercostal  (3/4  animals),  and  right  diaphragm  (4/5 
animals).  The  left  ECR  showed  greater  EMG  responses  in 
five  of  five  animals. 
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Fig.  4.  Impact  of  ISMS  on  tracheal  pressure.  A:  representative  examples  of  tracheal  pressure  and  left  diaphragm  EMG  activity  before  and  during  cervical  ISMS 
(100  Hz,  200  pA).  The  period  of  ISMS  is  indicated  by  the  solid  bars.  Note  the  negative  deflection  in  tracheal  pressure  during  ISMS.  B:  overlay  plot  further 
illustrates  the  tracheal  pressure  before  stimulation  (solid  line)  and  during  ISMS  (dashed  line).  C:  average  change  in  tracheal  pressure  during  lung  inflation  at 
prestimulation  baseline  and  during  ISMS.  Data  are  shown  for  the  spinal-intact  condition  (left)  and  after  acute  C2Hx  (right).  *P  <  0.05,  **P  <  0.005. 


Inspiratory-related  diaphragm  EMG  bursting  is  enhanced 
after  ISMS.  While  our  primary  intent  was  to  determine  whether 
ISMS  evoked  diaphragm  activity  during  the  period  of  stimu¬ 
lation,  we  noted  that  spontaneous  diaphragm  EMG  activity 
often  was  present  after  the  stimulus  was  turned  off.  Thus,  in  six 
of  eight  animals  after  acute  C2Hx,  the  hemidiaphragm  that  was 
electrically  silent  before  ISMS  showed  both  tonic  and  inspira¬ 


tory-related  activity  during  and  after  the  l-min  period  of 
stimulation  (Fig.  6 A).  In  these  experiments,  “activity”  was  F6 
defined  as  one  or  more  clearly  discernible  motor  unit  potentials 
that  were  discharging  phasically  during  the  inspiratory  period. 
The  duration  of  the  effect  was  variable  and  did  not  persist 
beyond  2  min  in  any  experiment.  On  average,  ISMS-induced 
spontaneous  diaphragm  motor  unit  activity  lasted  for  43  ±  38 


Fig.  5.  ISMS  after  subacute  C2Hx.  A:  repre¬ 
sentative  diaphragm  EMG  activity  after  sub¬ 
acute  C2Hx.  In  the  example  traces,  genio- 
glossus-triggered  ISMS  was  delivered  during 
the  breaths  marked  by  the  black  bars.  B : 
histological  section  of  the  C2  spinal  cord 
stained  with  cresyl  violet.  The  example  dem¬ 
onstrates  an  anatomically  complete  hemile- 
sion  extending  to  the  midline  of  cervical 
cord.  C:  stimulus-triggered  averages  from  the 
ipsilesional  (solid  line)  and  contralesional  di¬ 
aphragm  (dashed  line);  data  were  obtained 
from  the  period  indicated  by  the  dashed 
boxes  in  A.  These  traces  illustrate  activation 
of  the  diaphragm  ipsilateral  to  the  C2Hx 
lesion.  D :  average  change  in  tracheal  pres¬ 
sure  during  lung  inflation  at  prestimulation 
baseline  and  during  ISMS.  *P  <  0.05. 
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Fig.  6.  Example  recordings  illustrate  short¬ 
term  potentiation  of  ipsilateral  diaphragm 
EMG  activity  following  ISMS.  A:  tongue  and 
ipsilateral  (left)  diaphragm  EMG  activity  af¬ 
ter  acute  C2Hx,  during  ISMS,  and  immedi¬ 
ately  after  cessation  of  ISMS.  Note  that  clear 
phasic  (inspiratory)  activity  can  be  seen  after 
ISMS,  whereas  the  baseline  showed  no  such 
activity.  B\  a  similar  response  can  be  ob¬ 
served  in  a  subacute  C2Hx  animal.  In  these 
examples,  the  ISMS  was  triggered  by  the 
inspiratory  genioglossus  EMG  signal  and 
was  delivered  at  100  Hz  for  1  min  of  respi¬ 
ratory  efforts. 


respiratory  cycles  after  ISMS  was  terminated.  The  first  breath 
immediately  after  cessation  of  ISMS  typically  demonstrated 
the  greatest  amount  of  motor  unit  recruitment,  and  this  was 
followed  by  a  gradual  decrease  in  motor  unit  activity  over 
subsequent  respiratory  cycles.  In  addition,  we  noted  an  in¬ 
crease  in  tonic  activity  (i.e.,  activity  persisting  across  the  entire 
respiratory  cycle)  after  ISMS  that  followed  a  time  course 
similar  to  the  phasically  active  motor  units.  A  similar  response 
to  ISMS  was  observed  in  three  of  five  rats  after  subacute  C2Hx 
injury  (Fig.  6 B). 

DISCUSSION 

This  study  has  demonstrated  the  feasibility  of  using  endog¬ 
enous  medullary  output  to  trigger  ISMS  to  the  midcervical 
spinal  cord.  The  ISMS  paradigm  was  successful  in  recruiting 
diaphragm  motor  units  in  spinal-intact  animals  and  also  after 
both  acute  and  subacute  interruption  of  bulbospinal  inspiratory 
drive  to  the  phrenic  motoneuron  pool  following  high  cervical 
SCI.  Further  refinement  of  this  approach  may  enable  develop¬ 
ment  of  a  “respiratory  neuroprosthesis”  that,  by  virtue  of  the 
endogenous  trigger,  could  adapt  to  temporally  changing  met¬ 
abolic  demands. 

Endogenously  triggered  ISMS  induces  diaphragm  activation. 
The  approach  used  in  this  study  was  a  “closed-loop-like” 
strategy  for  activating  phrenic  motoneurons  in  the  cervical 
spinal  cord.  As  an  initial  proof  of  concept,  cervical  ISMS  was 
triggered  from  EMG  signals  recorded  from  the  base  of  the 
tongue.  Such  recordings  consist  largely  of  activity  from  the 


primary  tongue  protrudor  muscle,  the  genioglossus  (Fuller  et 
al.  1998).  Respiratory-related  activation  of  the  genioglossus 
occurs  in  humans  (Mateika  et  al.  1999;  Saboisky  et  al.  2007) 
and  animal  models  (Bailey  et  al.  2005;  Fregosi  and  Fuller 
1997)  and  acts  to  stiffen  and/or  dilate  the  oropharyngeal 
airway  and  thus  preserve  upper  airway  patency  during 
breathing.  The  respiratory-related  discharge  of  genioglossus 
motor  units  increases  during  chemoreceptor  stimulation  in  a 
manner  similar  to  respiratory  pump  muscles  such  as  the 
diaphragm  (Bailey  2011;  Fuller  et  al.  1998).  Genioglossus 
and  tongue  muscle  output  is  also  subject  to  vagal  modula¬ 
tion,  with  progressive  inhibition  occurring  as  the  lung  inflates 
(Bailey  et  al.  2001;  Fukuda  and  Honda  1982).  Collectively,  the 
literature  establishes  that  a  complex^but  coordinated  interaction 
between  the  genioglossus  and  the  respiratory  pump  musclesjs 
important  for  optimizing  airway  resistance  across  the  respi¬ 
ratory  cycle  and  maintaining  upper  airway  patency.  For  the 
present  study,  the  most  salient  point  is  that  an  inspiratory 
signal  was  recorded  from  the  genioglossus  that  provided  a 
marker  for  the  endogenous  respiratory  rhythm.  Thus  the 
animals’  own  “decision  to  breathe”  triggered  the  cervical 
ISMS. 

The  genioglossus  EMG-triggered  cervical  ISMS  approach 
resulted  in  diaphragm  motor  unit  activation  after  both  acute 
(i.e.,  minutes)  and  subacute  (i.e.,  days-weeks)  C2Hx  injury.  It 
was  important  to  repeat  the  studies  under  these  advanced  lesion 
settings  because  severed  descending  pathways  have  not  under¬ 
gone  degeneration  after  acute  injury  and  therefore  could  still  be 
capable  of  conducting  action  potentials  during  stimulation 
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(Gandevia  and  Kirkwood  2011;  Kowalski  et  al.  2013)  as  seen 
in  other  experimental  conditions  (Moldovan  et  al.  2009).  In¬ 
deed,  the  latencies  of  the  diaphragm  MUAPs  in  the  acute 
condition  (i.e.,  >1.6  ms)  were  not  consistent  with  direct 
phrenic  motoneuron  activation,  and  therefore  it  is  likely  that 
synaptic  inputs  to  the  phrenic  motoneuron  pool  were  being 
activated  via  ISMS.  Accordingly,  it  was  uncertain  if  the  ISMS 
approach  would  be  effective  at  later  C2Hx  postinjury  intervals 
when  degeneration  of  severed  bulbospinal  respiratory  axons 
would  be  well  underway.  The  results  showed  that  diaphragm 
MUAPs  could  still  be  effectively  evoked  after  subacute  injury 
but  with  a  trend  toward  longer  and  more  variable  diaphragm 
MUAP  latencies  compared  with  acute  C2Hx  animals.  In  addi¬ 
tion,  the  biomechanical  impact  of  the  ISMS  (indirectly  inferred 
from  changes  in  tracheal  pressure)  also  appeared  to  be  reduced 
after  subacute  vs.  acute  injury.  Refinement  of  the  ISMS  method 
to  enable  direct  phrenic  motoneuron  stimulation  will  likely  be 
necessary  to  optimize  diaphragm  activation  after  long-term 
cervical  SCI. 

An  interesting  and  potentially  important  observation  that 
was  not  part  of  our  original  hypothesis  was  that  ISMS-induced 
inspiratory  bursting  in  the  hemiparetic  diaphragm  persisted  for 
several  respiratory  cycles  after  cessation  of  the  stimulation. 
This  occurred  in  the  majority  of  acute  and  subacute  C2Hx 
lesions  and  is  similar  to  the  previously  described  phenomenon 
of  respiratory  short-term  potentiation  (STP)  (Lee  et  al.  2015). 
Exposure  to  brief  periods  of  reduced  oxygen  (hypoxia)  triggers 
phrenic  STP,  which  is  manifested  as  a  progressive  enhance¬ 
ment  of  activity  followed  by  a  slow  decline  to  baseline  levels 
after  removal  of  the  hypoxic  stimulus  (Lee  et  al.  2009;  Powell 
et  al.  1998).  We  previously  demonstrated  that  hypoxia  can 
induce  STP  of  phrenic  motor  output  in  rats  with  subacute 
C2Hx,  and  the  response  is  greater  in  the  nerve  ipsilateral  vs. 
contralateral  to  the  lesion  (Lee  et  al.  2015).  The  STP  in  C2Hx 
rats  appeared  to  reflect  recruitment  of  a  population  of  phrenic 
motoneurons  that  had  been  silenced  by  injury,  which  continued 
to  burst  beyond  the  period  of  hypoxic  stimulation.  It  may  be 
that  ISMS  can  trigger  similar  STP-like  mechanisms  of  short¬ 
term  plasticity  and  ultimately  have  value  in  the  context  of 
neurorehabilitation.  This  increased  phrenic  bursting  after  ces¬ 
sation  of  ISMS  also  may  be  analogous  to  longer-duration 
functional  recoveries  others  have  observed  with  chronic  ISMS 
in  other  motor  systems  (Kasten  et  al.  2013;  McPherson  et  al. 
2015). 

Methodological  caveats  including  off-target  effects  of  cervi¬ 
cal  ISMS.  The  method  described  here  is  not  truly  a  “closed- 
loop”  system  because  the  ISMS  intensity  did  not  vary  with  the 
relative  strength  of  endogenous  “respiratory  drive.”  Varying 
the  stimulus  intensity  or  duration  in  proportion  to  the  relative 
magnitude  of  the  EMG  burst  used  to  trigger  ISMS  (vs.  a  simple 
“onset  trigger”)  could  potentially  address  that  issue,  but  this 
could  also  prove  problematic  since  ISMS  needs  to  occur  at  the 
onset  (vs.  peak)  of  the  inspiratory  effort.  Another  caveat  is  that 
rats  were  mechanically  ventilated  because  in  our  experience 
spontaneously  breathing  and  anesthetized  rats  rapidly  develop 
arterial  hypercapnia  and  respiratory  acidosis.  However,  me¬ 
chanical  ventilation  will  impact  endogenous  respiratory  drive 
and  thereby  alter  the  interactions  between  spontaneous  and 
ISMS  input  to  the  phrenic  region.  It  should  also  be  noted  that 
tongue  EMG  activity  is  unlikely  to  provide  an  effective  ISMS 
trigger  signal  in  the  vagally  intact  and  spontaneously  breathing 


rodent  or  human.  It  is  possible  to  record  inspiratory-related 
discharge  from  the  tongue  muscles  in  the  awake  human  and 
rodent,  but  available  data  indicate  that  tongue  EMG  may  not 
provide  a  consistent  inspiratory  trigger  signal  (Bailey  2011; 
Sood  et  al.  2005).  In  the  present  experiments  on  anesthetized 
rats,  however,  the  tongue  EMG  output  enabled  a  rigorous  test 
of  the  fundamental  hypothesis,  and  we  suggest  that  future  work 
should  focus  on  additional  potential  sources  of  respiratory 
output  to  serve  as  endogenous  trigger  signals  or  alternative 
detection  algorithms  (Dow  et  al.  2006). 

Neural  substrate  considerations.  One  of  the  more  challeng¬ 
ing  issues  related  to  electrical  stimulation  of  the  spinal  cord  is 
identification  of  the  neural  substrate  being  affected  directly  or 
indirectly.  Here  our  immediate  goal  was  to  target  ISMS  of  gray 
matter  at  the  level  of  the  phrenic  motor  nucleus.  Whether 
phrenic  motoneurons  are  being  directly  activated  by  ISMS 
cannot  be  determined  from  the  present  findings,  although 
latencies  in  the  range  of  1.5-2. 5  ms  suggest  indirect  stimula¬ 
tion  via  polysynaptic  pathways.  In  the  subacute  C2Hx  group, 
stimulation  of  ipsilateral  bulbospinal  inputs  to  phrenic  mo¬ 
toneurons  can  be  effectively  ruled  out  because  these  pathways 
will  be  undergoing  degeneration.  It  is  possible  that  ISMS 
activated  commissural  bulbospinal  projections  from  the  con¬ 
tralateral  cord  [i.e.,  pathways  associated  with  the  “crossed 
phrenic  phenomenon”  (Goshgarian  2003;  Goshgarian  et  al. 
1991)].  However,  past  studies  in  which  intraspinal  stimulation 
of  the  cervical  cord  was  used  to  activate  crossed  phrenic 
pathways  reported  latencies  of  <1  ms  (Fuller  et  al.  2003),  and 
these  are  considerably  shorter  than  the  values  reported  here. 
Collectively,  the  evidence  is  most  consistent  with  indirect 
activation  of  phrenic  motoneurons  via  ISMS,  possibly  via 
spinal  prephrenic  interneurons  (Lane  2011)  or  even  ascending 
afferent  projections  originating  below  the  injury  (Decima  and 
von  Euler  1969). 

ISMS  is  capable  of  discrete  activation  of  circuitry  with 
minimal  undesired  physiological  effects  (Pikov  et  al.  2007). 
Here  we  observed  clear  responses  in  the  ipsilateral  diaphragm, 
as  intended,  but  also  off-target  motor  unit  recruitment.  Very 
small  and  inconsistent  activation  of  the  contralateral  dia¬ 
phragm,  intercostal  muscles,  and  genioglossus  were  present 
during  ISMS.  This  indicates  current  spread  beyond  the  target 
region  or  activation  of  other  circuits  via  interneurons  (Perlmut- 
ter  et  al.  1998).  For  example,  in  the  case  of  ISMS-evoked 
responses  in  the  contralateral  diaphragm  and  intercostal  mus¬ 
cles,  our  previous  transneuronal  tracing  studies  provide  an 
interneuronal  basis  for  those  off-target  responses  (Lane  et  al. 
2008).  Intercostal  and/or  contralateral  diaphragm  recruitment 
could  in  fact  be  advantageous  when  cervical  SCIs  involve 
much  larger-scale  compromise  of  diaphragm  and  inspiratory 
intercostal  function  compared  with  C2Hx.  Activation  of  the 
genioglossus  was  barely  above  the  threshold  for  detection  and 
probably  reflected  activation  of  ascending  projections  to  med¬ 
ullary  respiratory  control  neurons.  For  example,  we  previously 
reported  that  some  cervical  interneurons  associated  with  the 
phrenic  motor  circuit  have  ascending  projections  to  the  me¬ 
dulla  (Lane  et  al.  2008). 

The  relatively  stronger  coactivation  of  forelimb  muscles 
induced  by  ISMS  presents  a  more  challenging  technical  issue. 
As  with  the  ipsilateral  diaphragm,  constant-latency  entrain¬ 
ment  occurred  in  the  forelimb,  and  this  may  reflect  overlap  of 
phrenic  and  forelimb  circuits  along  the  rostro-caudal  axis 
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(Gonzalez-Rothi  et  al.  2015).  ECR  motoneurons,  however, 
appear  to  occupy  a  more  dorso-lateral  position  than  the  phrenic 
motoneuron  pool  (Tosolini  and  Morris  2012).  This  raises  the 
possibility  that  even  if  ISMS  was  precisely  restricted  to  phrenic 
motoneurons,  concomitant  excitation  of  ECR  axons  as  they 
course  toward  their  ventral  roots  could  still  occur.  Further 
investigations  are  needed  to  sort  out  the  neuroanatomical  basis 
for  direct  vs.  indirect  effects  of  ISMS.  The  latter  may  be 
especially  beneficial  by  activating  desired  motoneuron  pools  in 
a  more  natural  order  than  achieved  with  direct  stimulation. 

Conclusions.  Previous  studies  directed  at  spinal  cord  stim¬ 
ulation  and  the  respiratory  system  after  SCI  have  focused  on 
epidural  stimulation  as  a  means  of  enhancing  abdominal  or 
phrenic  motor  output.  Epidural  stimulation  has  been  successful 
at  recruiting  diaphragm  motor  units  in  animal  studies  (DiMarco 
and  Kowalski  2009,  2011;  Kowalski  et  al.  2013)  and  has  been 
used  to  enhance  cough  in  spinally  injured  humans  (DiMarco  et 
al.  2014).  Epidural  stimulation  also  shows  promise  in  the 
context  of  locomotor  function  after  SCI  in  both  animals  and 
humans  (Edgerton  and  Harkema  2011;  Rejc  et  al.  2015). 
Despite  these  successes,  there  is  unquestionably  a  need  for 
improvement  and  refinement  of  spinal  cord  stimulation  meth¬ 
ods.  In  that  regard,  the  initial  findings  from  our  study  raise  the 
possibility  of  utilizing  endogenous,  physiologically  relevant 
respiratory  signals,  which  entail  peripheral  and  chemosensory 
feedback,  for  triggering  neuromodulation  of  the  phrenic  motor 
circuit  and  alleviation  of  ventilatory  insufficiency  following 
mid-  to  high-cervical  SCI.  Many  ISMS  parameters  (Bamford 
and  Mushahwar  2011;  Giszter  2015;  Tator  et  al.  2012)  must  be 
refined,  however,  before  the  benefits  of  invasive  ISMS  proto¬ 
cols  to  promote  respiratory  improvements  after  cervical  SCI 
can  be  fully  evaluated  or  compared  with  other  approaches, 
high-frequency  epidural  stimulation  (Kowalski  et  al.  2013)  in 
particular.  It  also  is  possible  that  enq  spinal  stimulation  ap¬ 
proach  will  be  better  suited  for  some  SCI  cases  than  others 
Future  investigations  involving  awake,  spontaneously  breath¬ 
ing  animals  will  be  crucial  for  determining  whether  chronic 
ISMS  delivery  will  function  solely  as  a  neuroprosthetic  or  can 
be  used  short  term  to  promote  functional  and  anatomical 
neuroplasticity  (Moritz  et  al.  2007)  alone  or  in  combination 
with  other  therapeutic  approaches  leading  to  long-lasting  im¬ 
provement  in  respiratory  function. 
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52  Abstract 

53  Mid-cervical  spinal  interneurons  form  a  complex  and  diffuse  network  and  may  be  involved  in 

54  modulating  phrenic  motor  output.  The  intent  of  the  current  work  was  to  enable  a  better 

55  understanding  of  mid-cervical  “network  level”  connectivity  by  pairing  the  neurophysiological 

56  multi-electrode  array  (MEA)  data  with  histological  verification  of  the  recording  locations.  We  first 

57  developed  a  method  to  deliver  100  nA  currents  to  electroplate  silver  onto  and  subsequently 

58  deposit  silver  from  electrode  tips  after  obtaining  mid-cervical  (C3-5)  recordings  using  an  MEA  in 

59  anesthetized  and  ventilated  adult  rats.  Spinal  tissue  was  then  fixed,  harvested,  and 

60  histologically  processed  to  “develop”  the  deposited  silver.  Histological  studies  verified  that  the 

61  silver  deposition  method  discretely  labeled  (50  pm  resolution)  spinal  recording  locations 

62  between  laminae  IV-X  in  cervical  segments  C3-C5.  Using  correlative  techniques,  we  next  tested 

63  the  hypothesis  that  mid-cervical  neuronal  discharge  patterns  are  temporally  linked.  Cross 

64  correlation  histograms  produced  few  positive  peaks  (5.3%)  in  the  range  of  0  -  0.4  ms,  but  21.4% 

65  of  neuronal  pairs  had  correlogram  peaks  with  a  lag  of  >  0.6  ms.  These  results  are  consistent 

66  with  synchronous  discharge  involving  mono-  and  polysynaptic  connections  among  mid-cervical 

67  neurons.  We  conclude  that  there  is  a  high  degree  of  synaptic  connectivity  in  the  mid-cervical 

68  spinal  cord,  and  that  the  silver  labeling  method  can  reliably  mark  metal  electrode  recording  sites 

69  and  “map”  interneuron  populations,  thereby  providing  a  low  cost  and  effective  tool  for  use  in 

70  MEA  experiments.  We  suggest  that  this  method  will  be  useful  for  further  exploration  of  mid- 

71  cervical  network  connectivity. 


72  New  and  Noteworthy 


73  We  describe  a  method  that  reliably  identifies  the  locations  of  multi-electrode  array  (MEA) 

74  recording  sites  while  preserving  the  surrounding  tissue  for  immunohistochemistry.  To  our 

75  knowledge,  this  is  the  first  cost-effective  method  to  identify  the  anatomical  locations  of  neuronal 

76  ensembles  recorded  with  a  MEA  during  acute  preparations  without  the  requirement  of 

77  specialized  array  electrodes.  In  addition,  evaluation  of  activity  recorded  from  silver  labeled  sites 

78  revealed  a  previously  unappreciated  degree  of  connectivity  between  mid-cervical  interneurons. 


79  Introduction 


80  Mid-cervical  spinal  interneurons  form  a  complex  and  diffuse  network  that  is  synaptically 

81  coupled  to  both  respiratory  and  non-respiratory  motor  pools  (Gonzalez-Rothi  et  al.  2015;  Lane 

82  2011;  Lane  et  al.  2008b).  Several  groups  have  advanced  the  hypothesis  that  mid-cervical  spinal 

83  interneurons  can  modulate  phrenic  motoneuron  excitability  and  thereby  influence  the  neural 

84  control  of  the  diaphragm  (Bellingham  and  Lipski  1990;  Douse  and  Duffin  1993;  Lane  et  al. 

85  2008a;  Lane  et  al.  2008b;  Palisses  et  al.  1989).  While  there  is  some  direct  evidence  to  support 

86  this  hypothesis  (Marchenko  et  al.  2015;  Sandhu  et  al.  2015),  other  studies  have  concluded  the 

87  opposite  (Duffin  and  Iscoe  1996).  A  significant  hurdle  in  testing  that  specific  hypothesis,  or 

88  related  hypotheses  regarding  cervical  interneuronal  circuits,  is  the  difficulty  of  studying  the 

89  “functional  connectivity”  in  diffuse  spinal  cord  networks.  One  of  the  foremost  challenges  is 

90  simultaneously  recording  numerous  cells,  and  this  can  be  addressed  through  the  use  of  multi- 

91  electrode  arrays  (MEA).  The  MEA  approach  enables  simultaneous  recordings  of  multiple  sites, 

92  but  histologically  identifying  each  recording  location  (versus  the  electrode  track)  while  also 

93  preserving  tissue  integrity,  poses  a  further  challenge  (Borg  et  al.  2015;  Li  et  al.  2015;  Nuding  et 

94  al.  2015).  Thus,  the  initial  thrust  of  the  current  work  was  modification  and  validation  of  a  silver 

95  labeling  technique  (Spinelli  1975)  to  enable  post-recording  deposition  of  a  small  amount  of  silver 

96  (i.e. ,  for  histological  marking)  from  the  tip  of  each  electrode  in  an  MEA.  Additionally,  we 

97  developed  an  electrical  circuit  to  enable  the  use  of  small  currents  (100  nA)  for  silver 

98  electroplating  and  deposition  to  prevent  tissue  and  electrode  damage  associated  with  high 

99  levels  of  current  (Fung  et  al.  1998).  Using  the  electroplated  MEA,  we  recorded  discharge  from 

100  ensembles  of  neurons  in  the  mid-cervical  (C3-C5)  spinal  cord  in  adult  rats  and  demonstrated  a 

101  practical  application  of  this  technique  by  “matching”  the  array  electrodes  to  the  corresponding 

102  anatomical  locations  marked  by  silver. 


103  The  overall  intent  was  to  enable  a  better  understanding  of  mid-cervical  spinal  discharge 

104  and  “network  level”  connectivity  by  pairing  the  neurophysiological  MEA  data  with  histological 

105  verification  of  the  recording  locations.  Thus,  using  correlative  techniques  (Moore  et  al.  1970),  we 

106  tested  the  hypothesis  that  the  discharge  patterns  of  mid-cervical  (C3-4)  spinal  neurons  are 

107  temporally  linked  in  time  domains  consistent  with  mono-  and  poly-synaptic  connections.  In 

108  addition,  mid-cervical  neuronal  discharge  patterns  were  assessed  relative  to  the  endogenous 

109  inspiratory  pattern  -  measured  via  bilateral  phrenic  nerve  recordings  -  to  determine  if  bursting 

110  was  temporally  linked  to  phrenic  motoneuron  activity.  To  our  knowledge,  no  prior  study  has 

111  comprehensively  evaluated  the  temporal  characteristics  across  multiple  mid-cervical  neurons 

112  using  MEA  technology.  The  results  presented  herein  demonstrate  a  previously  unappreciated 

113  degree  of  connectivity  and  indicate  a  high  prevalence  of  temporally  related  discharge  patterns 

114  between  mid-cervical  interneurons  with  characteristics  consistent  with  mono-  and  polysynaptic 

115  connections;  and  also  provides  a  comprehensive  description  of  a  cost-effective  histological 

116  approach  for  validating  anatomical  locations  of  MEA  recording  sites. 

117  Materials  and  Methods 

1 1 8  Animals 

119  All  experiments  were  conducted  with  adult  Sprague-Dawley  rats  obtained  from  ENVIGO 

120  (formally  Harlan  Laboratories).  Most  experiments  (n=12)  were  performed  with  untreated,  spinal 

121  intact  rats.  A  subset  of  rats  (n=2)  received  a  cervical  spinal  injury  (C3/C4  lateralized  contusion, 

122  force:  205  KD,  displacement:  1225  pm;  Infinite  Horizon  pneumatic  impactor,  Precision  Systems 

123  &  Instrumentation)  using  published  methods  from  our  group  (Lane  et  al.  2012).  Spinal  injured 

124  rats  were  allowed  to  recover  for  12  weeks  prior  to  electrophysiology.  Histological  results  from 

125  spinal  injured  animals  were  used  to  determine  effectiveness  of  silver  labeling  and  to  compare 

126  micro-motor  depths  relative  to  silver  labeling  depths.  Since  there  is  evidence  suggesting 


127  neurophysiological  properties  such  as  connectivity  may  be  altered  following  SCI  (Lane  et  al. 

128  2009;  Lane  et  al.  2008b),  analysis  of  neurophysiology  was  limited  to  spinal  intact  animals  (n=4). 

129  All  rats  were  housed  in  pairs  in  a  controlled  environment  (1 2h  light/dark  cycles)  with  food  and 

130  water  ad  libitum.  All  experimental  protocols  were  approved  by  the  Institutional  Animal  Care  and 

131  Use  Committee  at  the  University  of  Florida. 

1 32  In  Vivo  Electrophysiology 

133  All  rats  were  anesthetized  with  3%  isoflurane  (in  100%  02)  and  transferred  to  a  heated 

134  station  where  established  surgical  methods  (Lee  and  Fuller  2010;  Lee  et  al.  2009;  Mahamed  et 

135  al.  2011;  Sandhu  et  al.  2015;  Streeter  and  Baker-Herman  2014a;  b;  Strey  et  al.  2012)  were 

136  used  to  set  up  electrophysiological  experiments.  Core  body  temperature  was  maintained  at  37  ± 

137  0.5°C  with  a  servo  controlled  heating  device  (CWE,  model  700  TC-1000).  A  tail  vein  catheter 

138  was  placed  for  /'.v  delivery  of  urethane  anesthesia  and  fluids.  The  trachea  was  cannulated  (PE- 

139  240  tubing)  and  rats  were  pump-ventilated  (Harvard  Apparatus,  Rodent  Ventilator  683;  volume: 

140  ~3-2.5  mL;  frequency:  70/minute).  Once  ventilated,  C02  was  added  (FiC02:  <  3%)  to  maintain 

141  end-tidal  C02  (EtC02)  between  -40-50  mmHg  throughout  the  protocol  (Capnogard; 

142  Respironics,  Inc.).  Tracheal  pressure  was  monitored  and  lungs  were  periodically  hyperinflated 

143  (2-3  breaths)  to  prevent  atelectasis  (~1/hour).  A  bilateral  vagotomy  was  performed.  Rats  were 

144  slowly  converted  (6ml/hr;  Harvard  Apparatus  syringe  pump)  to  urethane  anesthesia  (1.7  g/kg, 

145  i.v.  in  distilled  water)  and  isoflurane  was  withdrawn.  A  femoral  arterial  catheter  (PE-50  tubing) 

146  was  placed  to  monitor  blood  pressure  (CWE,  TA-100)  and  sample  blood  gases  (i-STATI 

147  Analyzer  Abbot)  throughout  the  protocol.  Using  a  dorsal  approach,  the  left  and  right  phrenic 

148  nerves  were  isolated,  cut  distally,  and  partially  desheathed  (-1/2  the  length  of  the  exposed 

149  nerve).  A  midline  incision  extending  from  the  base  of  the  skull  to  mid-thoracic  region  was  made 

150  and  spinal  vertebrae  C3-T2  were  exposed.  Using  a  nose  clamp  and  T2  spinous  process,  the  rat 

151  was  slightly  elevated  off  the  table  to  reduce  ventilator  induced  motion  artifact  and  level  the 


152  spinal  cord.  A  laminectomy  was  performed  from  C3-C6,  the  dura  was  cut  and  reflected  back, 

153  and  the  pia  was  gently  removed  at  the  MEA  insertion  point.  A  unilateral  pneumothorax  was 

154  performed  to  decrease  chest  wall  movement  and  positive  end-expiratory  pressure  (PEEP)  of 

155  ~1-2  cm  of  water  was  applied  to  prevent  atelectasis.  Animals  received  the  neuromuscular 

156  paralytic  pancuronium  bromide  (2.5  mg/kg,  i.v.,  Hospira,  Inc.)  to  eliminate  spontaneous 

157  breathing  efforts.  Adequate  depth  of  anesthesia  was  monitored  by  assessing  blood  pressure 

158  responses  to  toe-pinch  and  urethane  supplements  were  given  as  necessary.  Blood  pressure 

159  and  fluid  homeostasis  were  maintained  by  a  slow  infusion  of  a  1:3  solution  (8.4%  sodium 

160  bicarbonate/Lactated  Ringers,  i.v.). 

161  Bilateral  phrenic  nerve  output  was  recorded  using  custom  made  bipolar  suction 

162  electrodes  filled  with  0.9%  saline.  Compound  action  potentials  were  amplified  (X20  k,  Grass 

163  Instruments,  P511),  band-pass  filtered  (3  Hz  to  3  kHz),  digitized  (16-bit,  25k 

164  samples/second/channel;  Power1401,  Cambridge  Electronic  Design,  Cambridge,  UK  (CED)), 

165  and  integrated  (time  constant:  20ms)  in  Spike2.v8  software  (CED).  A  custom  made  MEA  similar 

166  to  that  originally  designed  at  the  University  of  South  Florida  by  Lindsey  BG  and  colleagues 

167  (Morris  1996)  and  used  in  our  previous  report  (Sandhu  et  al.  2015)  was  used  to  record  bilateral 

168  mid-cervical  (C3-C5)  spinal  activity  (Figure  1A).  The  array  contained  16  tungsten  electrodes 

169  coated  with  Epoxylite  Insulation  (impedance:  10  ±1  MO;  shank  diameter  125  pm;  tip  diameter 

170  <1  pm;  FHC,  #UEWLEGSE0N1  E)  (Figure  IB).  A  key  feature  of  this  array  is  the  ability  to 

171  independently  control  the  depth  of  each  of  the  electrodes  using  micro-motors.  This  greatly 

172  improves  the  recording  yield  by  allowing  the  experimenter  to  “hold”  a  recording  on  one  electrode 

173  while  continuing  to  search  for  neurons  with  other  electrodes.  The  array  was  mounted  on  a 

174  stereotaxic  frame  and  eight  electrodes  arranged  in  two  staggered  rows  of  four  were  placed  into 

175  each  hemi-cord  at  the  dorsal  root  entry  zone  (indicated  by  black  arrows  in  Figure  1C).  The  inner 

176  distance  between  the  two  sets  of  eight  electrodes  was  approximately  1  mm,  while  the  distance 

177  between  electrodes  within  each  row  was  -300  pm  (Figure  ID).  Electrodes  tips  were  maintained 


178  in  this  “fixed  matrix”  by  the  array  guide.  One  by  one,  electrodes  were  advanced  into  the  spinal 

179  cord  while  audio  was  monitored  until  single  units  with  ~3:1  signal  to  noise  ratio  were 

180  discriminated.  Neural  signals  from  single  electrodes  were  amplified  (X5k),  band-pass  filtered  (3- 

181  3kHz),  digitized  (16-bit,  25k  samples/second/channel;  Power  1401,  CED)  and  recorded  with 

182  Spike2.v8  software  (CED).  Once  phrenic  nerve  activity  and  spinal  discharge  was  stable, 

183  “baseline  activity”  was  recorded  with  Fi02  set  to  0.50.  Using  experimental  procedures  similar  to 

184  our  previous  studies  (Sandhu  et  al.  2015;  Strey  et  al.  2012),  the  inspiratory  gas  mixture  was 

185  altered  using  adjustable  flow  meters  to  expose  rats  to  a  5  min  period  of  hypoxia  (Fi02:  0.11), 

186  followed  by  15  min  of  baseline  oxygen  levels  (Fi02:  0.50).  Silver  deposition  was  performed 

187  immediately  after  the  neurophysiology  protocol.  An  arterial  blood  gas  was  sampled  during 

188  baseline  and  hypoxia  to  assure  blood  gases  were  within  physiological  limits. 

1 89  Analysis  of  Electrophysiological  Signals 

190  All  electrophysiological  data  was  collected  and  sorted  using  Spike2.v8  software  (CED). 

191  Extracellular  action  potentials  from  individual  neurons  were  extracted  from  continuous 

192  recordings  and  converted  to  waveforms  using  spike-sorting  tools  of  the  acquisition  software. 

193  Briefly,  spikes  were  clustered  using  at  least  80%  template  matching  and  principle  component 

194  analyses.  Sorted  spikes  were  exported  and  analyzed  using  custom  MATLAB  software 

195  (MathWorks  R2015a).  A  standard  set  of  analyses  (spike  interval  histograms,  cycle  triggered 

196  histograms,  and  spike-triggered  averages)  were  performed  to  electrophysiologically  phenotype 

197  each  neuron.  Interval  histograms  were  created  for  each  spike  train  to  assure  only  single  unit 

198  activity  was  represented  in  each  waveform.  Similar  to  previous  studies  (Galan  et  al.  2010; 

199  Sandhu  et  al.  2015),  cycle-triggered  histograms  were  used  to  assess  the  preferred  respiratory 

200  modulation  of  each  recorded  neuron  relative  to  the  phase  of  the  respiratory  cycle  defined  by 

201  phrenic  nerve  activity  (Cohen  1968).  Using  the  integrated  phrenic  nerve  output,  the  beginning 

202  and  end  of  the  inspiratory,  and  therefore  expiratory  phases,  were  calculated  as  a  departure  of 


203  phrenic  nerve  activity  >  15  standard  deviations  above  the  average  activity  during  the  expiratory 

204  phase.  Cycle-triggered  histograms  were  constructed  for  each  neuron  by  dividing  the  respiratory 

205  period  into  20  bins  of  equal  size  and  spikes  were  counted  within  each  bin  and  summated  over 

206  50  consecutive  breaths  during  baseline,  hypoxia,  and  post-hypoxia.  To  determine  if  neurons 

207  were  respiratory  modulated,  the  spikes  occurring  during  bins  of  the  inspiratory  phase  were 

208  separated  from  those  occurring  during  the  expiratory  phase  and  the  Wilcoxon  signed-rank  test 

209  was  used  to  test  the  null  hypothesis  (i.e.,  no  difference  between  inspiration  and  expiration). 

210  Cervical  spinal  neurons  were  classified  into  four  categories:  1)  neurons  that  discharged  primarily 

211  during  the  inspiratory  period  (i.e.,  inspiratory  modulated;  waveform  1  Figure  7C,  D);  2)  neurons 

212  that  discharged  primarily  during  the  expiratory  period  (i.e.,  expiratory  modulated;  waveform  1 

213  Figure  7H,  I);  3)  neurons  that  discharged  without  respiratory  modulation  (i.e.,  tonic;  waveform  2 

214  Figure  7C,  D,  H,  and  I);  and  4)  neurons  which  ceased  firing  at  time  points  after  baseline  were 

215  labeled  as  inactive  at  that  time  point. 

216  Spike-triggered  averaging  of  raw  and  rectified  phrenic  nerve  activity,  was  used  to 

217  examine  the  temporal  relationship  of  neuronal  waveforms  and  phrenic  motor  output  (Lipski  et  al. 

218  1983).  Short-latency,  offset  peaks  in  the  raw  and  rectified  phrenic  nerve  average  provided 

219  evidence  that  the  recorded  neuron  was  a  phrenic  motoneuron  (Christakos  et  al.  1994).  If 

220  features  were  not  detected  using  spike-triggered  averaging,  the  recorded  cell  was  classified  as 

221  a  spinal  interneuron.  Cross-correlation  histograms  were  constructed  for  all  possible  pairs  of 

222  simultaneously  recorded  neurons  using  a  bin  width  of  0.2ms  to  evaluate  functional  connectivity 

223  (Moore  et  al.  1970).  Similar  to  published  methods  (Aertsen  and  Gerstein  1985),  the  detectability 

224  index  (Dl)  was  calculated  for  each  cycle-triggered  histogram  as  the  peak  relative  to  average 

225  background  activity  (calculated  during  the  interval  -15  to  -3  ms  prior  to  zero),  divided  by  the 

226  standard  deviation.  Only  significant  features  occurring  with  a  positive  lag  (e.g.,  >  0)  were 

227  counted.  Features  were  considered  significant  if  the  Dl  was  >  3  (Meissen  and  Epping  1987). 

228  Central  peaks  (0-0.4  ms)  between  a  pair  of  neurons  supported  synchronous  firing  between 


229  neurons;  whereas  offset  peaks  (>  0.6  ms)  indicated  functional  excitation  between  the  trigger  and 

230  target  neurons. 

231  Plating  Electrodes  and  Depositing  Silver 

232  Methods  for  electroplating  and  depositing  silver  described  previously  (Spinelli  1975) 

233  were  adapted  for  use  with  MEA  recordings.  For  each  electrode,  silver  was  electroplated  (Figure 

234  2A)  and  deposited  (Figure  2B)  using  tightly  controlled  and  monitored  DC  circuits.  The  circuits 

235  consisted  of  an  electrode  selection  interface  (i.e.,  break  out  box),  two  9V  batteries,  a  custom 

236  built  p-current  control  device  (Bare  Electronics,  Gainesville,  FL),  a  p-current  precision  nA 

237  current  measurement  assistant  v3  (EEVblog  Adafruit  Industries,  #882),  and  a  digital  multimeter 

238  (Ideal  Industries,  #61-340)  all  connected  with  copper  patch  cables.  Specifications  of  the 

239  pcurrent  control  circuit  are  provided  in  Figure  2C.  Prior  to  each  experiment,  15  of  the  16  fine 

240  wire  tungsten  electrodes  were  electroplated  with  dissolved  silver  cations  in  an  electrolyte 

241  solution  (AgN03/KCN)  and  an  un-insulated  silver  wire  (100  mm,  0.025  in  diameter,  A-M 

242  Systems,  #783500)  was  used  as  the  anode  in  the  solution  (Figure  2A).  The  electrolyte  solution 

243  was  prepared  by  mixing  equal  parts  of  1%  potassium  cyanide  and  1%  silver  nitrate.  The  circuit 

244  was  tested  before  each  use  by  electroplating  silver  onto  copper  wire  to  assure  proper 

245  connections  and  current  settings.  Immediately  following  the  testing  procedure,  the  tip  of  each 

246  individual  electrode  was  introduced  to  the  electrolyte  solution  (as  the  cathode)  for  100  sec  at 

247  100  nA.  Electrodes  used  to  deposit  silver  were  re-plated  prior  to  each  experiment.  In  a  subset  of 

248  experiments  (n=3),  2  of  the  16  electrodes  were  coated  with  the  lipophilic  fluorescent  dye,  Dil 

249  (1,1'-Dioctadecyl-3,3,3',3'-Tetramethylindocarbocyanine  Perchlorate  ('Dil';  DilCI  8(3))) 

250  (Thermofisher;  #D-282)  after  electrodes  were  electroplated  with  silver  and  prior  to 

251  electrophysiology.  Each  electrode  was  repeatedly  dipped  (~10  times)  in  Dil  (50  mg/ml  in 

252  ethanol)  and  allowed  to  dry  for  5  sec  between  each  dip  (DiCarlo  et  al.  1996). 


253  Subsequent  to  ensemble  recordings  of  individual  neurons,  electrodes  were  left  in  place 

254  while  the  MEA  was  disconnected  from  the  acquisition  system  and  attached  to  the  electrode 

255  selection  interface  (i.e.,  breakout  box)  (Figure  2B).  Sites  to  be  labeled  were  chosen  based  on 

256  the  quality  of  the  signals  during  the  electrophysiologic  recordings.  The  circuitry  was  connected 

257  as  in  the  plating  procedure  with  the  exception  that  electrode  16  was  used  as  the  cathode  and 

258  the  individual  electrode  chosen  for  labeling  was  treated  as  the  anode.  Electroplated  silver  was 

259  released  from  the  electrode  tip  using  a  100  nA  current  for  10-100  sec.  Tissue  was  fixed 

260  approximately  15-30  minutes  after  silver  was  deposited  (described  below). 

261  Tissue  Preparation 

262  Animals  were  transcardially  perfused  with  ice  cold  saline  followed  by  4% 

263  paraformaldehyde  (PFA;  Electron  Microscopy  Sciences,  #19210)  in  IX  Dulbecco’s  phosphate 

264  buffered  saline  (DPBS,  Mediatech,  Inc.,  #21-030-CV).  Spinal  tissue  was  harvested  and  post 

265  fixed  in  4%  PFA  overnight  at  4°C.  The  cervical  spinal  cord  was  subsequently  blocked  (C3-C6; 

266  ~8-9  mm  of  tissue)  and  cryoprotected  in  30%  sucrose  in  IX  DPBS  overnight  at  4°C.  The  spinal 

267  cord  was  embedded  in  optimal  cutting  temperature  compound  (OCT,  Fisher  Scientific,  #23-730- 

268  571),  flash-frozen  using  2-methylbutane  and  dry  ice,  and  transversely  cut  at  40  pm  on  a  cryostat 

269  (Microm  HM  500,  GMI).  Sections  for  silver  labeling  were  placed  into  12-well  trays  (Corning®, 

270  #3737)  containing  Corning®  Netwell™  inserts  (2  sections  per  well;  15  mm  dia.,  mesh  size  of  74 

271  pm,  Corning®,  #3477)  filled  with  glass-distilled  water  (~4  ml.  per  well).  Sections  for 

272  immunohistochemistry  were  placed  in  96  well  plates  containing  IX  DPBS. 

273  Silver  Staining  Protocol 

274  Silver  staining  was  performed  immediately  after  sectioning.  Due  to  the  sensitivity  of  the 

275  reaction,  all  trays,  glassware,  and  instruments  used  for  staining  protocol  were  thoroughly  acid 

276  washed  (1%  HCI)  prior  to  staining  to  eliminate  contamination  and  remove  non-specific  catalysts. 


277  Plastic  12-well  trays  containing  the  developer  solution  were  discarded  after  used  for  a  single 

278  tray  of  tissue,  while  all  other  trays  were  discarded  after  each  animal.  Netwell  mesh  inserts  were 

279  soaked  in  bleach  between  each  animal.  The  deposited  silver  was  “developed”  as  previously 

280  described  (Spinelli  1975).  Briefly,  stock  solutions  were  prepared  as  outlined  in  Table  1.  The 

281  hydrogen  peroxide  solution  was  covered  in  foil  throughout  entire  protocol  due  to  its  light 

282  sensitivity  and  the  hydroquinone  and  citric  acid  solution  was  made  in  the  same  bottle.  The 

283  developer  was  made  fresh  for  each  12-well  tray  of  tissue  by  mixing  the  three  stock  solutions 

284  (acacia,  silver  nitrate,  and  hydroquinone/citric  acid)  five  minutes  before  use  (at  the  start  of  step 

285  9).  All  steps  were  carried  out  in  a  fume  hood,  and  light  exposure  was  limited  while  mixing 

286  developer  solution  and  during  “developing”  by  turning  off  the  fume  hood  light. 

287  The  solutions  were  poured  into  either  200  ml.  glass  petri  dishes  (Corning®, 

288  #3160152BO)  or  12-well  plastic  trays  according  to  the  step  (see  Table  1)  and  arranged  in  order 

289  of  the  staining  sequence.  To  maintain  precision  during  the  staining  protocol,  staining  was 

290  performed  in  “rounds”,  where  each  round  contained  three  trays  of  tissue  (~1  mm  tissue/tray). 

291  Trays  were  staggered  by  5  min,  and  all  three  trays  of  tissue  were  stained  from  start  to  finish. 

292  Rounds  were  repeated  until  all  tissue  was  stained.  The  Netwell  mesh  inserts  were  transported 

293  through  the  staining  sequence  using  a  Netwell  carrier  kit  with  handles  (Corning,  #3520)  to  allow 

294  simultaneous  processing  and  assure  all  wells  were  in  each  solution  for  the  appropriate  time. 

295  Following  staining,  tissue  was  stored  at  4°C  in  Netwells  filled  with  dH20.  Sections  were  mounted 

296  onto  charged  microscope  slides  (Fisher  Scientific™,  #12-550-15)  using  IX  PBS  and  allowed  to 

297  dry  for  48-72  hrs.  Sections  were  counterstained  with  cresyl  violet  covered  with  mounting 

298  medium  (Thermo  Scientific™,  #4111)  and  coverslipped.  Images  were  captured  with  bright-field 

299  microscopy  using  a  Microscope  Axio  lmager.A2  (Carl  Zeiss  Microscopy).  Images  were  stitched 

300  and  white  balanced  using  Adobe  Photoshop.  Anatomical  locations  of  the  silver  labeled  sites 

301  within  the  spinal  laminae  were  identified  using  Sengul,  Watson,  Tanaka,  and  Paxinos  Atlas  of 

302  the  spinal  cord  (Sengul  et  al.  2012).  All  representative  spinal  cord  images  reflect  Camera 


303  Lucida-style  drawings  of  our  histological  images  or  C4  spinal  section  in  the  Atlas  of  the  spinal 

304  cord  (Sengul  et  al.  2012).  The  depth  of  silver  labeled  sites  were  determined  using  10X 

305  photomicrographs  and  adjusted  by  10%  to  account  for  tissue  contraction  during  the  fixation  and 

306  histological  processing  based  on  estimates  of  prior  investigations  (Deutsch  and  Hillman  1977; 

307  Quester  and  Schroder  1997). 

308  Immunohistochemistry 

309  Initial  immunohistochemistry  (IHC)  experiments  determined  the  order  in  which  IHC 

310  should  be  performed  to  achieve  optimal  fluorescent  staining  using  tissue  which  did  not  contain 

31 1  deposited  silver.  IHC  was  performed  either  prior  to,  or  immediately  after  silver  staining  using  two 

312  common  markers  (e.g.  GFAP  and  NeuN).  Using  the  optimal  staining  order,  a  subsequent 

313  experiment  was  performed  with  tissue  containing  deposited  silver.  All  IHC  was  performed  with 

314  free-floating  sections  in  96  well  plates.  All  tissue  was  washed  with  IX  PBS  with  0.2%  Triton  (3  * 

315  5  min).  For  NeuN  staining,  tissue  was  incubated  for  40  min  at  room  temp  in  10%  normal  goat 

316  serum  (NGS)  in  IX  PBS  with  0.2%  Triton.  Spinal  sections  were  then  incubated  for  1  hr  at  room 

317  temp  followed  by  72hrs  at  4°C  in  the  primary  antibody  solution:  10%  NGS  in  IX  PBS  with  0.2% 

318  Triton  and  NeuN  (mouse  anti-NeuN,  1:1400;  Millipore,  #MAB377).  Tissue  was  washed  with  IX 

319  PBS  (3><5  min)  and  incubated  in  the  secondary  antibody  solution  for  1  hr  at  room  temp:  10% 

320  NGS  in  IX  PBS  and  Alexa  Fluor®  488  (goat  anti-mouse,  1:500;  Invitrogen,  #A1 1029).  Following 

321  incubation,  tissue  was  washed  with  IX  PBS  (3x5  min).  For  GFAP  staining,  tissue  was 

322  incubated  for  60  min  at  room  temp  in  5%  normal  goat  serum  (NGS)  in  IX  PBS  with  0.2%  Triton. 

323  Spinal  sections  were  then  incubated  for  1  hr  at  room  temp,  and  overnight  at  4°C  in  the  primary 

324  antibody  solution:  5%  NGS  in  IX  PBS  with  0.2%  Triton  and  GFAP  (mouse  anti-GFAP,  1:200; 

325  Sigma-Aldrich,  G8393).  Tissue  was  washed  with  IX  PBS  with  0.2%  Triton  (3x5  min)  then 

326  incubated  in  the  secondary  antibody  solution  for  1  hr  at  room  temp:  5%  NGS  in  IX  PBS  with 

327  0.2%  Triton  and  Alexa  Fluor®  594  (goat  anti-mouse  1:500;  Life  Technologies,  #A11005). 


328  Following  incubation,  tissue  was  washed  with  cold  IX  PBS  (3><5  min).  Sections  were  mounted 

329  on  positively  charged  glass  slides  (Fisher  Scientific™,  #12-550-15),  covered  with 

330  VECTASHIELD  Antifade  mounting  medium  (Vector  Laboratories,  #H1-1200)  and  cover-slipped. 

331  Slides  were  air-dried  and  stored  at  4°C.  All  fluorescent  images  were  captured  at  10X,  20X  and 

332  40X  magnification  with  a  Microscope  Axio  lmager.A2  (Carl  Zeiss  Microscopy). 

333  Results 

334  Optimization  of  Silver  Deposition  from  MEAs  and  Subsequent  Histological  Development 

335  to  Enable  Visualization  of  Recording  Locations 

336  Initial  experiments  determined  the  parameters  needed  to  deposit  and  identify  silver  in 

337  the  spinal  cord  without  emphasis  on  recording  neuronal  bursting  from  the  electrodes.  As  the  DC 

338  resistance  of  the  commercially  purchased  recording  electrodes  was  constant  and  the  current 

339  used  for  electroplating  silver  onto  the  electrodes  was  standardized  (as  described  in  the 

340  methods),  we  assessed  the  impact  of  two  critical  variables  on  the  histological  appearance  of 

341  silver  in  the  spinal  cord:  1)  the  duration  of  deposition  current,  and  2)  the  duration  of  histological 

342  development.  Figure  3  shows  representative  images  of  silver  labeling  obtained  at  different 

343  durations  of  deposition  and  histological  development  which  demonstrate  the  variability  of  silver 

344  staining  achieved  by  altering  these  parameters.  Consistent  and  discrete  labeling  at  a  resolution 

345  between  -25-50  pm  was  obtained  when  current  used  to  deposit  silver  was  applied  for  50  sec 

346  and  tissue  was  developed  for  7-10  minutes  (Figure  3A  and  B).  When  deposition  current  of  less 

347  than  50  sec  was  applied  few  apparent  labels  were  detected,  presumably  due  to  an  inadequate 

348  amount  of  silver  ions  necessary  for  visualization.  In  contrast,  deposition  of  80-100  sec  resulted 

349  in  a  “halo  effect”  which  increased  the  silver  labeled  area  to  -100-150  pm  (Figure  3C  and  D), 

350  most  likely  due  to  a  broad  distribution  of  deposited  silver.  Similarly,  histological  development  of 

351  less  than  7  minutes  was  not  long  enough  to  “develop”  the  deposited  silver  required  to  visualize 


352  the  label;  whereas  prolonged  exposure  to  the  developing  solution  tended  to  darken  large  areas 

353  of  the  section  thus  confounding  differentiation  due  to  insufficient  contrast  (data  not  shown). 

354  We  next  verified  that  the  sites  of  silver  deposition  represented  the  location  of  the  tip  (i.e., 

355  recording  location)  of  a  given  electrode.  In  a  subset  of  animals  (n=2  with  silver  processing,  n=1 

356  unprocessed  positive  control),  two  of  the  sixteen  silver  plated  electrodes  were  coated  with  the 

357  fluorescent  lipophilic  dye,  Dil  to  visualize  the  tracks  associated  with  the  electrode  insertion 

358  (DiCarlo  et  al.  1996;  Marton  et  al.  2016;  Naselaris  et  al.  2005).  While  Dil  was  readily  detected  in 

359  unprocessed  tissue  (data  not  shown),  Dil  could  not  be  detected  after  the  histological  processing 

360  required  for  silver  labeling.  Therefore,  we  conclude  these  two  techniques  are  not  compatible. 

361  However,  in  several  histological  sections  an  “electrode  track”  could  be  identified  terminating  at 

362  the  silver  labeled  site  (Figure  3E).  This  is  consistent  with  the  interpretation  that  the  silver  labeled 

363  site  is  located  at  the  tip  of  the  electrode. 

364  Another  objective  was  to  determine  if  the  silver  labeling  method  could  be  coupled  with 

365  fluorescent  immunohistochemistry  (IHC)  procedures.  Accordingly,  we  performed  IHC  with  two 

366  markers  commonly  used  in  the  central  nervous  system  (GFAP  and  NeuN)  either  prior  to,  or 

367  immediately  after  histological  processing  for  silver  development.  These  markers  were  chosen 

368  due  to  their  ability  to  label  astrocytes  (Bignami  et  al.  1972)  and  neurons  (Wolf  et  al.  1996). 

369  Fluorescent  staining  was  detected  in  both  cases,  but  was  more  robust  when  IHC  was  performed 

370  after  the  silver  staining  procedures  (data  not  shown),  indicating  IHC  (at  least  for  these 

371  commonly  used  markers)  is  compatible  with  the  histological  processing  steps  necessary  to 

372  “develop”  the  silver.  Using  this  staining  order,  we  determined  the  location  of  deposited  silver 

373  relative  to  NeuN  positive  cells.  Representative  photomicrographs  from  a  C4  spinal  section 

374  depicting  positive  silver  labeling  in  the  intermediate  gray  is  shown  in  Figure  4A.  Fluorescent 

375  labeling  of  the  same  section  stained  with  NeuN  and  DAPI  are  shown  in  Figure  4B.  High 

376  magnification  images  of  the  silver  labeled  site  suggest  the  recording  electrode  was  in  close 


377  proximity  to  a  NeuN  positive  neuron  (Figure  4B).  These  qualitative  results  demonstrate  that  the 

378  silver  labeling  technique  can  be  coupled  with  fluorescent  IHC. 

379  Using  Silver  Labeling  to  Identify  Anatomical  Locations  of  Metal  Array  Electrodes 

380  In  six  animals  (n=4  spinal  intact  and  n=2  SCI),  we  used  the  optimized  deposition  and 

381  development  parameters  discussed  previously  to  label  select  electrodes.  Overall,  using  this 

382  approach  to  identify  the  anatomical  locations  of  the  recording  electrodes,  we  obtained  positive 

383  labeling  in  35  of  39  attempts.  This  indicates  that  the  silver  labeling  method  can  be  used  to 

384  identify  the  anatomical  locations  of  array  electrodes  with  a  success  rate  of  approximately  90%. 

385  Positive  silver  labeling  was  identified  in  the  spinal  gray  matter  between  C3-C5.  The  locations  of 

386  each  silver  labeled  site  in  spinal  intact  and  SCI  animals  was  plotted  (unilaterally  for  simplicity) 

387  according  to  their  anatomical  locations  (Figure  5A).  Silver  labeled  sites  were  identified  between 

388  lamina  IV-X,  with  the  greatest  number  found  in  lamina  VII  (Figure  5A). 

389  The  success  of  the  discrete  silver  labeling  afforded  the  chance  to  compare  the  intended 

390  recording  location  (i.e. ,  the  micro-motor  coordinates  used  during  the  neurophysiology 

391  experiment)  with  actual  the  location  of  the  labeling  (Figure  5B,  C).  Each  silver  labeled  site  was 

392  measured  from  the  dorsal  surface  of  the  cord.  Based  on  estimates  of  prior  investigations 

393  regarding  tissue  shrinkage  during  paraformaldehyde  fixation  and  subsequent  tissue  processing 

394  (Deutsch  and  Hillman  1977;  Quester  and  Schroder  1997),  the  measured  histological  depth  of 

395  each  silver  labeled  site  was  adjusted  by  10%.  This  analyses  indicated  that  the  micro-motor 

396  depth  is  likely  to  overestimate  the  actual  depth  of  the  electrode,  and  perhaps  even  more 

397  importantly  that  this  relationship  is  altered  by  SCI.  Linear  regression  analysis  indicated  a 

398  significant  relationship  between  recording  coordinates  and  histological  staining  in  spinal-intact 

399  animals  (motor:  1.6  ±  0.04  mm;  histological:  1.2  ±  0.05  mm;  P  =  0.0003),  but  not  after  chronic 

400  SCI  (motor:  1.4  ±  0.07  mm;  histological:  0.6  ±  0.05  mm;  P  =  0.3014;  Figure  5B).  The  depth  of 

401  the  silver  labeled  site  differed  from  the  micro-motor  coordinate  in  both  groups  (Figure  5C),  but 


402  was  significantly  increased  following  SCI  (0.8  ±  0.07)  compared  to  spinal  intact  animals  (0.4  ± 

403  0.04).  These  data  highlight  the  need  to  histologically  identify  multi-electrode  recording  sites 

404  especially  following  experimental  conditions  such  as  SCI  when  tissue  fibrosis  and  scarring  can 

405  be  expected  to  alter  electrode  movement  within  the  CNS.  Future,  studies  can  apply  the  linear  fit 

406  as  a  proxy  to  more  accurately  calculate  the  “actual  depth”  of  the  electrode  tips. 

407  Coupling  Spinal  Discharge  with  Anatomical  Locations 

408  All  silver  electroplated  electrodes  were  capable  of  discriminating  single  units,  indicating 

409  that  the  pretreatment  had  a  minimal,  if  any,  impact  on  the  ability  to  record  and  discriminate 

410  extracellular  signals.  Three  silver  labeled  sites  were  identified  within  lamina  IX  of  the  ventral 

41 1  horn  (Figure  5A  and  6A).  Spike-triggered  averaging  (STA)  of  the  raw  and  rectified  phrenic  nerve 

412  activity  in  relation  to  neuronal  discharge  produced  a  distinct  peak  with  an  average  lag  time  of 

413  0.45  ±  0.16  ms  (Figure  6E),  and  therefore  indicated  that  the  recorded  discharge  was  from 

414  phrenic  motoneurons  (Christakos  et  al.  1994;  Mitchell  et  al.  1992).  All  three  of  these  cells  were 

415  active  primarily  during  the  inspiratory  phase  (Figure  6B  and  C),  which  is  the  typical  firing  pattern 

416  of  phrenic  motoneurons  in  this  preparation  (Sandhu  et  al.  2015).  The  silver  labeling  data  verified 

417  that  the  recordings  were  from  the  region  of  the  phrenic  motor  nucleus  (Furicchia  and 

418  Goshgarian  1987;  Goshgarian  and  Rafols  1981;  Kinkead  et  al.  1998;  Mantilla  et  al.  2009; 

419  Prakash  et  al.  2000).  This  close  matching  between  neurophysiological  data  and  anatomical  data 

420  illustrates  the  effectiveness  of  the  silver  labeling  method. 

421  The  remaining  silver  labeled  sites  were  located  between  laminae  IV-X  (Figure  7).  Spike- 

422  triggered  averaging  provided  no  evidence  of  discharge  synchrony  in  relation  to  phrenic  motor 

423  output;  thus  these  28  recordings  were  considered  to  represent  interneurons  (Figure  7E  and  J). 

424  At  baseline,  most  recorded  interneurons  (18/28  or  66%;  Figure  7K)  fired  tonically  throughout  the 

425  respiratory  cycle  (e.g.,  waveform  2  in  Figure  7C,  D,  H  and  I).  These  tonically  discharging  cells 

426  were  not  restricted  to  a  particular  lamina,  but  rather  were  recorded  throughout  the  cervical  grey 


427  matter  (Figure  7L).  A  smaller  proportion  of  interneurons  (5/28  or  18%;  Figure  7K)  primarily  fired 

428  during  the  inspiratory  phase  (waveform  1  in  Figure  7C  and  D),  and  these  cells  were  found  in 

429  lamina  VI  (n=1),  VII  (n=3),  and  X  (n=1).  During  hypoxia,  two  of  these  neurons  (in  lamina  VII  and 

430  X)  switched  to  a  tonic  firing  pattern,  but  then  resumed  an  inspiratory  pattern  post-hypoxia.  In 

431  addition,  one  inspiratory  neuron  in  lamina  VII  was  inspiratory  modulated  during  baseline  and 

432  hypoxia,  but  then  assumed  a  tonic  firing  pattern  post-hypoxia.  Similarly,  during  baseline 

433  conditions,  five  cells  (5/28  or  18%;  Figure  7K)  discharged  primarily  during  the  expiratory  period 

434  (e.g.,  waveform  1  in  Figure  7H  and  I)  and  were  located  in  lamina  VIII  (n=3),  IX  (n=1),  and  X 

435  (n=1).  During  hypoxia,  three  of  these  neurons  (in  lamina  VIII  (n=2)  and  IX  (n=1))  switched  to  a 

436  tonic  firing  pattern  which  was  then  maintained  post-hypoxia.  In  addition,  one  expiratory- 

437  modulated  neuron  in  lamina  VIII  ceased  bursting  immediately  post-hypoxia.  Therefore,  due  to 

438  phase  switching,  the  proportion  of  tonic  firing  interneurons  increased  to  82%  (23/28)  during 

439  hypoxia  and  remained  elevated  post-hypoxia  (20/28  or  71%).  Interneuron  bursting  patterns 

440  during  baseline,  hypoxia  and  post-hypoxia  are  summarized  in  Figure  7K,  and  the  corresponding 

441  anatomical  locations  are  provided  in  Figure  7L.  In  this  sample  of  neurons,  respiratory-related 

442  discharge  was  observed  only  in  the  ventral  grey  matter  (i.e.,  lamina  VI,  VII,  VII,  IX  and  X),  and 

443  tonic  discharge  was  recorded  throughout  the  mid-cervical  grey  matter. 

444  Mapping  anatomical  locations  and  quantitation  of  interneuronal  discharge  across  mid- 

445  cervical  spinal  laminae  are  shown  in  Figure  8.  A  representative  example  from  one  recording  in 

446  which  silver  labeling  was  used  to  “map”  the  anatomical  locations  of  the  array  electrodes  and 

447  corresponding  electrophysiological  data  is  shown  in  Figure  8A-C.  A  Camera  Lucida-style 

448  drawing  of  the  anatomical  locations  of  eight  silver  labeled  electrodes  was  constructed  (Figure 

449  8B).  Numerical  identification  of  the  anatomical  positions  of  each  electrode  corresponds  to  the 

450  schematic  presented  in  Figure  ID  and  defines  the  rostral-caudal  and  medial-lateral  positions  of 

451  the  array  electrodes.  Corresponding  mid-cervical  spinal  discharge  and  integrated  phrenic  nerve 

452  activity  during  baseline  and  hypoxia  on  the  left  and  right  hemi-cord  is  shown  in  Figure  8C.  The 


453  majority  of  these  recordings  were  from  tonic  firing  interneurons  (electrodes  3,  8,  9,  11,  12,  and 

454  15)  and  one  represented  a  phase  switching  (e.g.  from  tonic  at  baseline  to  expiratory  during 

455  hypoxia)  interneuron  (electrode  6;  Figure  8C).  Using  this  technique,  a  summary  of  the 

456  anatomical  locations  of  all  silver  labeled  interneurons  was  constructed  (Figure  8D).  In  addition, 

457  the  average  discharge  frequency  (Hz)  of  all  silver  labeled  interneurons  was  presented  in  Figure 

458  8E.  The  results  indicate  a  dorsal-ventral  discharge  gradient,  with  higher  discharge  rates  in 

459  lamina  IV  and  X  and  lower  values  in  lamina  VIII  and  IX.  This  discharge  map  of  the  mid-cervical 

460  spinal  network  demonstrates  the  utility  of  using  the  silver  labeling  technique  to  identify  the 

461  anatomical  locations  of  ensembles  of  interneurons  recorded  with  a  MEA.  These  results  show 

462  the  utility  of  using  this  technique  to  standardize  the  sampling  distribution  of  recording  locations 

463  across  experimental  groups. 

464  Temporal  Relationships  between  Cervical  Interneuron  Discharge 

465  Anatomical  data  indicate  that  cervical  interneurons  are  part  of  a  diffuse  and  synaptically 

466  coupled  propriospinal  network  (Lane  et  al.  2008b),  but  relatively  little  is  known  about  functional 

467  connectivity  between  cervical  interneurons.  Therefore,  initial  analyses  focused  on  all  recorded 

468  neurons  in  spinal  intact  rats  (i.e.,  silver  and  non-silver  labeled  neurons)  to  screen  for  short  time 

469  scale  (i.e.,  0-1 0ms)  features  (Aertsen  and  Gerstein  1985;  Aertsen  et  al.  1989;  Meissen  and 

470  Epping  1987).  Significant  features  were  identified  as  departures  in  the  cross-correlation 

471  histogram  >  3  standard  deviations  from  the  background  noise  (Meissen  and  Epping  1987). 

472  Using  this  approach,  significant  central  (i.e.,  no  lag  time  from  zero)  and  offset  correlogram 

473  peaks  were  detected  (Figure  9A).  A  summary  of  the  latency  of  significant  features  relative  to  the 

474  trigger  (i.e.,  time  zero)  are  shown  in  Figure  9B. 

475  Inspection  of  the  data  suggested  that  latencies  for  correlogram  peaks  were  distributed  in 

476  three  ranges:  1)  <  0.4  ms;  2)  0.6-2. 8ms,  and  3)  >  3.0  ms  (discrete  integers  reflect  0.2ms  bins 

477  used  to  construct  histograms;  Figure  9B  and  C).  Few  correlations  (37/704  or  5.3%  of  possible 


478  correlations)  had  a  peak  with  latency  <  0.4  ms  which  is  consistent  with  excitation  from  a 

479  common  synaptic  input  (Aertsen  and  Gerstein  1985;  Kirkwood  1979;  Kirkwood  et  al.  1991; 

480  Meissen  and  Epping  1987).  A  greater  number  of  positive  correlations  had  latencies  between 

481  0.6-2. 8  ms  and  >  3.0  ms  (54/704  or  7.7%  and  97/704  or  13.8%,  respectively).  Correlogram 

482  peaks  with  latencies  above  0.6  ms  are  consistent  with  the  interpretation  of  functional  excitation 

483  involving  mono-,  or  polysynaptic  connections  (Aertsen  and  Gerstein  1985;  Kirkwood  1979; 

484  Moore  et  al.  1970).  We  next  evaluated  the  number  of  positive  correlations  relative  to  the 

485  recording  locations.  For  this,  neuronal  pairs  were  classified  as  “unilateral”  (i.e. ,  recordings  on 

486  the  same  side  of  the  cord)  or  “bilateral”  (i.e.,  recordings  on  opposite  sides  of  the  cord). 

487  Evaluation  of  unilateral  pairs  showed  the  average  number  of  significant  correlogram  peaks  was 

488  20  ±  3.9  correlations  per  animal.  A  similar  value  was  obtained  for  bilateral  neuronal  pairs  (25  ± 

489  7.5  correlations  per  animal;  Figure  9D).  When  normalized  to  the  total  number  of  possible 

490  unilateral  (332)  or  bilateral  (372)  correlations,  similar  values  were  observed  (24  ±  0.5%  and  26  ± 

491  5.5%  respectively;  Figure  9E).  No  significant  differences  in  latencies  were  observed  when  data 

492  was  separated  into  unilateral  and  bilateral  neuronal  pairs  (data  not  shown). 

493  Neuronal  pairs  with  significant  correlogram  peaks  that  were  also  silver  labeled  were 

494  used  to  construct  correlation  summary  maps  to  illustrate  the  anatomical  locations  of  each 

495  neuron  (Figure  10).  The  number  of  positive  correlations  within  each  lamina  was  normalized  to 

496  the  total  number  of  neurons  present  in  the  corresponding  lamina.  Unilateral  neuronal  pairs  are 

497  shown  in  Figure  10A  and  bilateral  pairs  are  shown  in  Figure  10B.  Of  the  positive  unilateral 

498  correlations,  interneurons  in  dorsal  lamina  (i.e.,  IV,  V  and  VI),  made  and  received  the  greatest 

499  number  of  excitatory  connections  (Figure  10A).  In  contrast,  when  the  trigger  and  target  neuron 

500  were  on  opposite  sides  of  the  spinal  cord,  interneurons  in  lamina  V,  VIII  and  IX  made  and 

501  received  the  most  connections  (Figure  10B). 

502  We  also  found  33  troughs  in  cross-correlograms  which  were  significantly  different  than 

503  background  activity  using  the  same  criteria  used  to  detect  peaks  (data  not  shown).  For 


504  comparison  purposes,  the  number  of  troughs  occurring  <  0.4  was  20/704  or  2.8%;  whereas 

505  those  between  0.6-2. 8  ms  was  1/704  or  0.1%;  and  those  >  3.0  ms  was  12/704  or  1.7%. 

506  Correlogram  troughs  with  these  latencies  are  generally  thought  to  reflect  functional  inhibition 

507  (Aertsen  and  Gerstein  1985;  Kirkwood  1979). 

508  Discussion 

509  Here  we  describe  a  silver  labeling  technique  that  reliably  identifies  the  position  of  MEA 

510  tips  following  in  vivo  neurophysiologic  recordings.  The  results  indicate  that  silver  labeling  can  be 

51 1  histologically  identified  with  a  resolution  of  ~50  pm  while  also  preserving  the  structural  integrity 

512  of  the  tissue  with  a  90%  success  rate.  A  practical  application  of  this  technique  was 

513  demonstrated  by  matching  the  location  of  mid-cervical  MEA  recording  sites  with  neuronal 

514  discharge  to  create  maps  of  bursting  patterns  (e.g.,  Figures  7  and  8).  Moreover,  using 

515  correlative  techniques,  we  also  provide  to  our  knowledge  the  first  neurophysiological  evidence 

516  for  the  presence  of  extensive  connections  between  interneurons  in  mid-cervical  grey  matter 

517  (e.g.,  Figures  9  and  10).  Coupling  silver  labeling  with  MEA  neurophysiology  provides  a  powerful 

518  tool  to  investigate  network  level  properties  of  the  spinal  cord. 

51 9  Commentary  Regarding  the  Silver  Labeling  Method 

520  The  first  challenge  associated  with  adapting  the  silver  labeling  method  (Spinelli  1975)  for 

521  MEAs  was  accurately  measuring  the  low  currents  (nA)  needed  for  electroplating  and  depositing 

522  silver  from  tungsten  array  electrodes.  Although  devices  which  measure  low  currents  are 

523  commercially  available  (e.g.,  for  patchclamp/iontophoresis),  they  are  expensive  and  not  easily 

524  adapted  to  electroplating  or  electrodeposition.  Additionally,  multimeters  cannot  resolve  currents 

525  in  the  nanoamp  range  and  have  limited  accuracy  due  to  voltage  drop  caused  by  current  flow 

526  through  the  current  measuring  device  (i.e. ,  “burden  voltage”;  (Jones  2010).  To  overcome  these 

527  limitations,  we  utilized  an  EEV  precision  p-current  assistant  which  permits  nanoamp  currents  to 


528  be  measured  on  a  standard  multimeter  and  also  improves  the  accuracy  to  0.2%  (Jones  2010). 

529  We  also  designed  a  p-current  control  device  capable  of  generating  100  nA  while  also  allowing 

530  current  adjustments  of  less  than  1  nA  (circuit  details  provided  in  Figure  2C).  The  ability  to 

531  precisely  produce  and  control  this  current  was  critical  since  passing  higher  currents  can  alter 

532  electrode  impedance  and  recording  characteristics  (Fung  et  al.  1998).  The  current  generating 

533  system  was  specifically  designed  to  be  compatible  with  both  single  electrodes  and  MEAs. 

534  Published  methods  for  marking  the  location  of  MEAs  include  iontophoresis  of  neural 

535  tracers  (Fekete  et  al.  2015;  Haidarliu  et  al.  1999;  Kovacs  et  al.  2005),  topical  application  of 

536  fluorescent  dyes  (DiCarlo  et  al.  1996;  Naselaris  et  al.  2005),  lesioning  (Brozoski  et  al.  2006; 

537  Townsend  et  al.  2002) ,  electrical  imaging  (Li  et  al.  2015),  and  imaging-based  approaches  (Borg 

538  et  al.  2015;  Fung  et  al.  1998;  Koyano  et  al.  201 1 ;  Matsui  et  al.  2007).  The  silver  labeling  method 

539  offers  several  advantages  over  these  previous  techniques.  First,  this  technique  can  be  used 

540  with  the  standard  metal  electrodes  typically  used  with  MEAs  (Borg  et  al.  2015;  Cogan  2008). 

541  This  is  in  contrast  to  the  specialized  electrodes  necessary  for  iontophoresis  (Fekete  et  al.  2015; 

542  Kovacs  et  al.  2005;  Naselaris  et  al.  2005).  Second,  the  silver  labeling  method  does  not  require 

543  implanted  electrodes  or  expensive  and  often  unavailable  equipment  such  as  MRI  to  visualize 

544  electrode  placement  (Borg  et  al.  2015;  Fung  et  al.  1998;  Koyano  et  al.  2011;  Matsui  et  al.  2007). 

545  Third,  the  histologically  identified  silver  deposition  spanned  an  area  of  ~50  pm,  and  this 

546  provides  an  approximately  3X  greater  resolution  as  compared  to  previously  published 

547  techniques  compatible  with  metal  MEAs.  For  example,  MRI  based  approaches  enables 

548  visualization  of  microelectrodes  sites  at  an  in-plane  resolution  of  150-200  pm  (Fung  et  al.  1998; 

549  Koyano  et  al.  2011;  Matsui  et  al.  2007).  Fluorescent  dyes  achieve  resolutions  between  50-400 

550  pm,  but  identify  electrode  track  trajectories  rather  than  discrete  locations  of  electrode  tips 

551  (DiCarlo  et  al.  1996;  Naselaris  et  al.  2005).  Finally,  whereas  lesioning  techniques  can  damage 

552  surrounding  tissues  (Townsend  et  al.  2002),  low  current  levels  required  for  silver  deposition 


553  preserves  the  integrity  of  the  surrounding  tissue.  This  is  a  particularly  important  point  since  it 

554  permits  neurochemical  phenotyping  of  neurons  at  or  near  the  electrode  tip. 

555  One  of  the  features  enabling  “matching”  of  electrophysiological  signals  with  anatomical 

556  locations  is  that  the  relative  position  of  each  electrode  was  known  since  the  electrodes  exist  in  a 

557  fixed  matrix  (Figure  ID).  Nevertheless,  with  multiple  silver  labeled  sites  in  a  relatively  discrete 

558  area,  matching  each  silver  labeled  site  with  its  corresponding  recording  electrode  can  pose  a 

559  challenge.  To  overcome  this  hurdle,  several  safeguards  were  implemented.  First,  since 

560  recordings  were  performed  with  a  bilateral  recording  array,  a  small  longitudinal  cut  in  the  dorsal 

561  spinal  cord  was  made  prior  to  sectioning  to  demarcate  the  left  vs.  right  hemi-cord.  Second,  the 

562  rostral-caudal  orientation  of  spinal  sections  was  maintained  throughout  the  histological  process 

563  by  sectioning  sequentially  into  12-well  plates.  Most  importantly,  the  histologically  identified  silver 

564  deposition  was  linked  to  the  corresponding  recording  data  by  using  the  lateral  -  medial  aspect 

565  of  the  electrode  matrix.  The  phrenic  motoneuron  data  provided  further  verification  of  the  ability 

566  to  match  silver  labeling  with  the  neurophysiology  data.  The  anatomical  location  of  phrenic 

567  motoneurons  is  well  defined  in  rat  (Goshgarian  and  Rafols  1981;  Mantilla  et  al.  2009;  Prakash  et 

568  al.  2000;  Zhan  et  al.  1989),  and  the  electrophysiological  analyses  (e.g.  STA)  confirmed  an 

569  appropriate  “match”  between  the  anatomical  location  of  silver  labeling  (i.e,  ventral  grey  matter, 

570  lamina  IX)  and  the  recording  electrode.  Moreover,  cells  with  histological  features  consistent 

571  phrenic  motoneurons  (Furicchia  and  Goshgarian  1987;  Prakash  et  al.  2000)  were  clearly 

572  identified  near  the  silver  label  (-50-100  pm). 

573  To  assure  that  our  recordings  (and  silver  labeled  sites)  were  from  neurons  (vs.  fibers  of 

574  passage)  in  the  immediate  vicinity  of  electrode,  we  used  high  impedance  tungsten  wires  coated 

575  with  Epoxylite  Insulation  with  an  exposed  tip  of  <  1  pm.  We  suggest  that  using  these  insulated 

576  high  impedance  electrodes,  the  following  assumptions  are  reasonable:  1)  the  recorded  neural 

577  discharge  is  of  somatic  origin  ( versus  axonal );  and  2)  the  amplitude  of  the  recorded  action 

578  potentials  decrease  with  distance  between  the  soma  and  the  recording  electrode  tip.  The 


579  reasons  for  these  assumptions  are  as  follows.  Extracellular  action  potentials  recorded  from 

580  axons  are  very  brief  (i.e. ,  0.1 -0.3  ms),  small,  and  extremely  sensitive  to  electrode  movement 

581  (Bellingham  and  Lipski  1990;  Kirkwood  et  al.  1988).  In  contrast,  action  potentials  recorded  in  the 

582  current  study  were  broad  (i.e.,  >  1  ms),  biphasic,  and  the  recordings  were  stable  as  electrodes 

583  were  moved  up  to  100  pm.  These  features  are  consistent  with  extracellularly  recorded  somatic 

584  action  potentials  (Bellingham  and  Lipski  1990;  Nelson  1959).  The  literature  indicates  that  the 

585  amplitude  of  extracellularly  recorded  action  potentials  is  inversely  proportion  to  the  distance 

586  between  the  soma  and  the  recording  electrode  (Kirkwood  et  al.  1988;  Nelson  1959).  Therefore, 

587  by  achieving  a  3:1  signal-to-noise  ratio  when  placing  each  electrode,  we  can  be  confident  that 

588  the  electrode  tip  is  close  to  the  cell  soma  generating  the  largest  action  potential.  Lastly,  the 

589  electrode  design  described  above  reduces  the  surface  area  available  for  silver  plating,  and  thus 

590  silver  is  deposited  from  only  the  ~1  pm  plated  length  at  each  electrode  tip. 

591  The  silver  labeling  method  was  not  100%  successful,  and  we  suspect  the  small  number 

592  of  failures  (N=4  of  39  possible  sites)  occurred  due  to  a  lack  of  deposited  silver  or  tissue  damage 

593  during  histological  processing  (e.g.,  tearing  of  sections  or  error  during  the  placement  of 

594  longitudinal  notch  in  the  spinal  cord).  A  final  methodological  commentary  relates  to  the 

595  difference  between  the  histologically  verified  “depth”  of  the  electrode  and  the  coordinates  used 

596  in  the  micro-drive  while  placing  each  electrode  (e.g.,  Figure  5).  The  micro-drive  coordinates 

597  overestimated  the  actual  depth  of  the  recording,  and  probably  reflects  some  degree  of 

598  “pillowing”  of  the  tissue  as  the  electrodes  were  advanced.  The  most  salient  point,  however,  is 

599  that  the  difference  between  histological  and  micro-motor  coordinates  was  increased  in  rats  with 

600  chronic  SCI  (e.g.,  Figure  5C).  These  data  highlight  the  need  to  histologically  identify  multi- 

601  electrode  recording  sites  following  SCI  when  tissue  fibrosis  and  scarring  can  be  expected  to 

602  impede  electrode  movement  within  the  spinal  cord. 


603  Mid-Cervical  Spinal  Interneurons  and  Respiratory  Motor  Output 

604  Several  laboratories  have  suggested  that  phrenic  motoneuron  discharge  can  be 

605  modulated  by  synaptic  input  from  spinal  interneurons  (Bellingham  and  Lipski  1990;  Davies  et  al. 

606  1985;  Sandhu  et  al.  2015).  Anatomical  studies  show  spinal  interneurons  are  uniquely  situated  to 

607  modulate  phrenic  motor  output  sine  they  are  in  close  opposition  with  medullary  projections 

608  (Davies  et  al.  1985;  Fedorko  et  al.  1983;  Hayashi  et  al.  2003;  Lane  et  al.  2008b)  and  have 

609  synaptic  connections  to  phrenic  motoneurons  (Dobbins  and  Feldman  1994;  Lane  et  al.  2009; 

610  Lane  et  al.  2008b;  Lois  et  al.  2009;  Yates  et  al.  1999).  Electrophysiological  recordings  in 

611  multiple  species  have  demonstrated  that  cervical  interneurons  have  respiratory-related 

612  discharge  patterns  (Bellingham  and  Lipski  1990;  Duffin  and  Iscoe  1996;  Marchenko  et  al.  2015; 

613  Palisses  et  al.  1989;  Sandhu  et  al.  2015)  and  respond  to  respiratory  stimuli  such  as  phrenic 

614  afferent  stimulation  (Iscoe  and  Duffin  1996;  Speck  and  Revelette  1987)  and  hypoxia  (Sandhu  et 

615  al.  2015).  Our  results  add  to  this  literature  in  two  primary  ways.  First,  to  our  knowledge  only  one 

616  prior  study  (Marchenko  et  al.  2015),  has  combined  neurophysiology  with  the  histological 

617  identification  of  the  anatomical  location  of  respiratory-related  spinal  interneurons  associated  with 

618  the  phrenic  motor  pool.  Our  data  show  that  interneurons  with  respiratory  related  discharge  are 

619  located  throughout  laminae  VI,  VII,  VII,  IX,  and  X.  Second,  the  current  results  provide  new 

620  information  regarding  the  impact  of  respiratory  stimulation  with  hypoxia  on  mid-cervical 

621  interneuronal  discharge  patterns.  Indeed,  during/following  hypoxia  approximately  29%  of  spinal 

622  interneurons  alter  their  discharge  in  relation  to  the  respiratory  cycle.  We  noted  that  both 

623  inspiratory  and  expiratory  modulated  neurons  had  a  tendency  to  adopt  a  tonic  firing  pattern 

624  during  hypoxia.  However,  post-hypoxia  interneurons  with  baseline  inspiratory  discharge  tended 

625  to  resume  an  inspiratory  firing  pattern,  whereas  the  expiratory  interneurons  maintained  the  tonic 

626  pattern.  In  addition,  we  found  one  example  of  a  tonic  firing  interneuron  during  baseline  and 

627  hypoxia  which  become  inspiratory  modulated  post-hypoxia.  Neurons  exhibiting  “phase 


628  switching”  were  located  in  intermediate  and  ventral  laminae  (i.e.,  VII,  VIII  and  X)  rather  than 

629  dorsal  laminae.  To  our  knowledge,  this  is  the  first  report  to  describe  hypoxia  induced  phase 

630  switching  of  mid-cervical  interneuron  bursting  patterns.  When  neuronal  discharge  frequency 

631  was  evaluated  relative  to  anatomical  locations,  higher  frequencies  were  noted  in  lamina  IV  and 

632  X  (e.g.,  Figure  8E).  These  results  again  highlight  the  importance  of  matching  anatomical 

633  location  with  discharge  properties,  and  show  that  mid-cervical  interneuron  discharge  varies 

634  across  cervical  lamina  which  is  consistent  with  previous  descriptions  of  lumbosacral  neural 

635  activity  (Borowska  et  al.  2013;  Ruscheweyh  and  Sandkuhler  2002). 

636  Discharge  Synchrony  between  Mid-Cervical  Spinal  Interneurons 

637  Despite  both  historical  (Duffin  and  Iscoe  1996;  Palisses  et  al.  1989)  and  recent 

638  publications  (Lane  2011;  Marchenko  et  al.  2015;  Sandhu  et  al.  2015)  related  to  cervical 

639  respiratory  interneurons,  the  functional  connectivity  of  mid-cervical  spinal  circuit  remains  largely 

640  unexplored.  This  can  be  attributed  to  the  fact  that  studies  investigating  respiratory-related 

641  interneurons  in  the  spinal  cord  have  primarily  utilized  a  single  unit  recording  approach 

642  (Bellingham  and  Lipski  1990;  Lipski  and  Duffin  1986;  Lipski  et  al.  1993)  and  have  instead 

643  focused  on  connections  between  motor-  and  interneurons  (Davies  et  al.  1985;  Duffin  and  Iscoe 

644  1996;  Lipski  et  al.  1993).  In  the  current  study,  we  utilized  cross-correlation  analysis  to 

645  characterize  the  functional  connectivity  of  pairs  of  spinal  interneurons.  Similar  to  previous 

646  studies  investigating  spinal  interneurons  in  non-respiratory  related  networks  (Brown  et  al.  1979; 

647  Prut  and  Perlmutter  2003),  we  found  relatively  few  examples  of  synchronous  discharge 

648  consistent  with  a  shared  excitatory  pre-synaptic  input  (Aertsen  et  al.  1989;  Kirkwood  1979).  That 

649  is,  when  interneuronal  pairs  were  examined,  correlogram  peaks  with  latencies  between  0  -  0.4 

650  ms  were  only  observed  in  5.3%  of  recordings.  In  contrast,  21.4%  of  the  total  possible  positive 

651  correlations  had  latencies  >  0.6  ms.  These  longer  latency  peaks  are  consistent  with  mono-  and 

652  polysynaptic  connections  between  cervical  interneuron  pairs.  Of  the  positive  correlations 


653  involving  neurons  on  the  same  side  of  the  spinal  cord,  interneurons  in  dorsal  lamina  (i.e.,  IV,  V 

654  and  VI),  made  and  received  the  greatest  number  of  excitatory  connections.  In  contrast,  when 

655  the  trigger  and  target  neuron  were  on  opposite  sides  of  the  spinal  cord,  interneurons  in  lamina 

656  V,  VII  and  IX  made  and  received  the  most  connections.  Taken  together,  our  results  provide 

657  support  for  the  hypothesis  that  cervical  interneurons  form  a  dynamic  network  which  is  capable 

658  of  rapid  reconfiguration  and  modulation  of  cervical  motor  outputs.  We  suggest  that  the  MEA 

659  silver  labeling  method  will  help  advance  our  understanding  of  the  functional  and  anatomical 

660  correlates  of  ensembles  of  spinal  neurons. 

661  Application  to  Spinal  Cord  Injury 

662  One  potential  application  of  MEA  technology  is  to  examine  how  SCI  alters  the  cervical 

663  spinal  networks.  Spinal  networks  undergo  substantial  remodeling  following  SCI  (Bareyre  et  al. 

664  2004;  Sperry  and  Goshgarian  1993)  and  interneurons  have  been  implicated  in  SCI-induced 

665  plasticity  and  motor  recovery  (Alilain  et  al.  2008;  Bareyre  et  al.  2004;  Harkema  2008).  However, 

666  relatively  little  is  known  regarding  their  contribution  to  the  spontaneous  recovery  of  phrenic 

667  output  following  SCI  (Lane  et  al.  2008b).  In  the  current  study,  we  performed  recordings  on  n=2 

668  rats  with  chronic  cervical  SCI.  The  intent  of  these  experiments  was  not  to  “map”  changes  in  the 

669  spinal  network  post-injury,  but  rather  to  confirm  that  the  method  developed  herein  could  be 

670  effectively  utilized  in  rats  with  SCI.  We  successfully  deposited  and  detected  silver  labeling 

671  following  SCI  using  the  same  parameters  optimized  in  spinal  intact  animals.  However,  we  noted 

672  a  greater  discrepancy  between  micro-motor  and  silver  labeling  depth  in  spinal  injured  animals 

673  (e.g.,  Figure  5C).  While  multiple  factors  may  contribute  to  this  disparity,  structural  changes 

674  resulting  from  fibrosis  and  scarring  induced  following  injury  (Cregg  et  al.  2014)  may  increase 

675  electrode  drag,  thereby  reducing  the  effectiveness  of  utilizing  the  electrode  micro-motor  to 

676  predict  the  actual  depth  of  the  electrode  tips.  Regardless  of  the  specific  cause,  our  results 

677  highlight  the  importance  of  histologically  verifying  the  locations  of  recorded  neurons,  especially 


678  following  experimental  conditions  such  as  SCI  which  may  alter  electrode  movement  within  the 

679  CNS.  We  suggest  that  future  studies  utilizing  the  methods  described  here  will  shed  light  on  the 

680  contribution  of  the  propriospinal  network  to  the  recovery  of  function  following  SCI. 

681  Conclusions 

682  The  methodology  described  here  provides  a  solution  to  the  fundamental  limitation  of 

683  using  MEAs  containing  standard  metal  electrodes  by  identifying  the  discrete  anatomical 

684  locations  of  the  recording  sites.  Although  prior  studies  have  described  methods  to  determine  the 

685  approximate  locations  of  MEA  electrodes  (Borg  et  al.  2015;  Brozoski  et  al.  2006;  DiCarlo  et  al. 

686  1996;  Fekete  et  al.  2015;  Haidarliu  et  al.  1999;  Kovacs  et  al.  2005;  Koyano  et  al.  2011; 

687  Naselaris  et  al.  2005;  Townsend  et  al.  2002),  to  date  labeling  the  discrete  recording  locations  of 

688  non-specialized  metal  array  electrodes  in  an  acute  preparation  while  preserving  the  surrounding 

689  tissue,  has  only  been  achieved  with  single-cell  recordings.  With  MEA  recordings,  higher  order 

690  analyses  can  be  used  to  investigate  the  functional  connectivity  of  the  neural  network  during 

691  physiological  stimuli  and  following  neurologic  injury  (e.g.,  Figures  9  and  10  also  see:  (Aertsen 

692  and  Gerstein  1985;  Kirkwood  1979;  Meissen  and  Epping  1987).  In  this  regard,  the  current 

693  results  are  consistent  with  a  high  degree  of  synaptic  connections  between  mid-cervical  neurons. 

694  Collectively,  our  experiments  show  that  1)  MEA  silver  labeling  method  enables  the 

695  electrophysiological  output  of  neural  networks  to  be  coupled  with  histological  verification  of 

696  electrode  locations;  2)  mid-cervical  interneurons  are  capable  of  rapid  “phase  switching”  of  burst 

697  patterns  relative  to  the  respiratory  cycle  during  and  after  hypoxia;  and  3)  a  high  percentage  of 

698  mid-cervical  interneuronal  pairs  have  temporally  related  discharge  patterns. 


699  Figure  Legends 

700  Figure  1:  Multi-electrode  array.  A:  Multi-electrode  recording  array  containing  16  tungsten  fine 

701  wire  electrodes  each  controlled  by  micro-motors  and  held  in  place  by  the  array  guide.  B:  High 

702  resolution  image  of  the  electrode  tips  maintained  in  a  “fixed  matrix”  by  the  array  guide.  C:  Image 

703  of  the  electrodes  positioned  in  the  dorsal  C4  spinal  cord  at  the  dorsal  root  entry  zone  (e.g.,  black 

704  arrows).  D:  Schematic  of  the  recording  matrix  consisting  of  eight  electrodes  arranged  in  two 

705  staggered  rows  of  four.  The  inner  distance  between  the  two  sets  of  eight  electrodes  was 

706  approximately  1  mm  and  the  distance  between  electrodes  within  each  row  was  -300  pm. 

707  Figure  2:  Silver  electroporation  and  deposition  circuits.  Both  circuits  consisted  of  an 

708  electrode  selection  interface  ( i.e .  breakout  box),  two  9V  batteries  connected  in  series,  a  custom 

709  built  p-current  control  device  (constant  current  board),  a  p-current  precision  nA  current 

710  measurement  assistant  (p-current  Box),  and  a  digital  multimeter  connected  via  copper  patch 

71 1  cables  (thick  black  lines).  A:  For  the  silver  electroplating  circuit,  the  array  was  connected  as  the 

712  cathode  to  the  electrode  selection  box  (breakout  box)  while  an  un-insulated  silver  wire  was 

713  placed  in  the  electrolyte  solution  (KCN/AgN03)  and  used  as  the  anode.  The  p-current  box  was 

714  connected  to  breakout  box  corresponding  to  the  electrode  to  be  plated  (Chi-15).  B:  For  the 

715  silver  deposition  circuit,  the  polarity  of  the  circuit  was  reversed  with  breakout  box  used  to  select 

716  the  reference  channel  (Cathode/Chi 6)  and  the  electrode  to  be  deposited  (Anode/Chl-15).  C: 

717  The  p-current  control  device  utilizes  an  NPN  transistor  (Q1  PN2222A)  to  regulate  current  from 

718  the  power  source  (two  9V  batteries;  B1  and  B2)  through  a  load  between  J1  and  J2.  Two  small- 

719  signal  diodes  D3  and  D4  (1N4148)  provide  a  fixed  voltage  (~1.25VDC)  used  in  parallel  with 

720  voltage  divider  (R4).  This  adjustable  voltage  divider  controls  the  base  current  delivered  to  the 

721  transistor  so  that  the  emitter  current  of  the  transistor  produces  a  stable  current  of  lOOnA.  As 

722  load  between  J1  and  J2  changes,  the  transistor  actively  alters  the  voltage  -  and  therefore 


723  current  -  drop  across  collector  (C)  to  emitter  (E),  maintaining  a  constant  voltage  (across  R3). 

724  Current  was  monitored  by  an  EEVblog  p-current  precision  nA  current  measurement  assistant  v3 

725  with  a  consumer  grade  multi-meter  and  adjusted  to  maintain  100  nA. 

726  Figure  3:  Representative  images  of  silver  labeling.  A-D:  Representative  histological  sections 

727  from  the  cervical  spinal  cord  containing  silver  labeled  sites  counterstained  with  cresyl  violet  and 

728  high  resolution  images  of  the  boxed  areas.  These  images  depict  the  variability  of  silver  labeling 

729  at  different  durations  of  silver  deposition  (50-100  sec)  and  histological  development  (7-10 

730  minutes).  E:  Images  of  a  cervical  spinal  section  depicting  an  electrode  track  (black  arrows) 

731  coursing  through  the  spinal  tissue  and  terminating  at  the  silver  labeled  site.  CC:  central  canal; 

732  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  A-D:  0.5  mm  and  50  pm  (callouts);  E:  0.5  mm, 

733  100  pm,  and  50  pm,  respectively. 

734  Figure  4:  Silver  labeling  coupled  with  fluorescent  immunohistochemistry.  A: 

735  Representative  cervical  spinal  section  of  positive  silver  labeling  ( brown/black )  and  high 

736  resolution  image  of  the  boxed  area.  B:  Fluorescent  labeling  of  neurons  stained  with  NeuN 

737  (green)  and  nuclei  stained  with  DAPI  (blue)  of  the  same  section  presented  in  panel  A,  and  high 

738  resolution  image  of  the  boxed  area  illustrating  the  silver  labeled  site  was  in  close  proximity  of  a 

739  NeuN  positive  cell.  CC:  central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  0.5  mm 

740  and  50  pm  (callouts). 

741  Figure  5:  Electrode  and  silver  labeled  depth.  A:  Representative  C4  section  of  silver  labeled 

742  sites  in  spinally  intact  (•)  and  spinally  injured  (o)  animals  projected  onto  one  side  of  the  cord  for 

743  simplicity.  Each  lamina  was  shaded  according  to  the  number  of  silver  labeled  sites  within  that 

744  lamina.  B:  Scatter  plot  of  the  micro-motor  depths  versus  the  depths  of  the  corresponding  silver 

745  labeled  sites  in  spinally  intact  (•)  and  spinally  injured  (o)  animals.  Linear  regression  analysis 


746  was  applied  to  determine  the  line  of  best-fit  and  linear  equation  for  each  group.  The  line  of 

747  identity  is  displayed  to  indicate  the  location  where  micro-motor  and  silver  labeling  depths  are 

748  equal.  C:  Normalized  depth  calculated  as  the  difference  between  the  micro-motor  depth  silver 

749  labeled  site  for  spinally  intact  (•)  and  spinally  injured  (o)  animals.  *  P<0.001  Unpaired  T-test. 

750  Figure  6:  Silver  labeling  of  a  phrenic  motoneuron.  A:  Representative  C4  section  containing 

751  a  silver  labeled  site  counterstained  with  cresyl  violet  and  high  resolution  image  containing  the 

752  silver  labeled  site  located  in  the  medial  aspect  of  the  ventral  horn  in  lamina  IX  (callout).  B: 

753  Corresponding  raw  spinal  discharge  and  integrated  phrenic  output  during  baseline,  hypoxia,  and 

754  post-hypoxia  depicting  discharge  during  the  inspiratory  phase.  C:  Integrated  phrenic  output,  raw 

755  neuronal  discharge,  and  sorted  spikes  (waveform)  during  hypoxia.  D:  Cycle-triggered  histogram 

756  during  50  consecutive  breaths  overlaid  with  the  average  integrated  phrenic  waveform  during 

757  hypoxia.  E:  Spike-triggered  average  of  the  raw  and  rectified  phrenic  nerve  depicting  a  delay  of 

758  0.76  ms.  CC:  central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  A  0.5mm  and  50pm 

759  (callout);  B,  C:  1.5  sec. 

760  Figure  7:  Silver  labeling  of  mid-cervical  spinal  interneurons.  A:  Representative  C4  section 

761  containing  a  silver  labeled  site  counterstained  with  cresyl  violet  and  high  resolution  image 

762  containing  the  silver  labeled  site  (callout)  B:  Corresponding  neuronal  output  and  integrated 

763  phrenic  output  during  baseline,  hypoxia,  and  post-hypoxia  depicting  a  single  tonic  firing  neuron 

764  at  baseline  and  recruitment  of  a  phasic  inspiratory  neuron  during  hypoxia  and  post-hypoxia.  C: 

765  Integrated  phrenic  output,  raw  neuronal  discharge,  and  sorted  spikes  (waveforms)  during 

766  hypoxia.  D:  Cycle-triggered  histograms  of  both  neurons  during  50  consecutive  breaths  overlaid 

767  with  the  average  integrated  phrenic  waveform  during  hypoxia.  E:  Spike-triggered  average  of  the 

768  raw  and  rectified  phrenic  nerve  depicting  a  lack  of  positive  features.  F:  Representative  C4 

769  section  containing  a  silver  labeled  site  counterstained  with  cresyl  violet  and  high  resolution 


770  image  containing  the  silver  labeled  site  (callout)  G:  Corresponding  neuronal  output  and 

771  integrated  phrenic  output  during  baseline,  hypoxia,  and  post-hypoxia  depicting  an  expiratory 

772  and  tonic  firing  neuron  at  baseline  and  hypoxia,  and  only  the  tonic  neuron  post-hypoxia.  H: 

773  Integrated  phrenic  output,  raw  neuronal  discharge,  and  sorted  spikes  (waveforms)  during 

774  hypoxia.  I:  Cycle  triggered  histograms  of  both  neurons  during  50  consecutive  breaths  overlaid 

775  with  the  average  integrated  phrenic  waveform.  J:  Spike  triggered  average  of  the  raw  and 

776  rectified  phrenic  nerve  depicting  a  lack  of  positive  features.  K:  Stacked  bar  graphs  depicting  the 

777  number  of  each  cell  type  at  baseline,  hypoxia,  and  post-hypoxia  each  normalized  to  the  total 

778  number  of  interneurons  (n=28).  L:  Anatomical  location  of  each  silver  labeled  interneuron 

779  projected  onto  one  side  of  the  cord  for  simplicity,  identified  by  the  bursting  pattern  during 

780  baseline.  CC:  central  canal;  DH:  dorsal  horn;  VH:  ventral  horn.  Scale  bars:  A,  F,  L:  0.5  mm  and 

781  50  pm  (callouts),  B,  C,  G,  H:  1.5  sec. 

782  Figure  8:  Mapping  anatomical  locations  of  mid-cervical  spinal  neurons.  A:  Representative 

783  photomicrograph  of  a  C4  section  containing  a  silver  labeled  site  (from  electrode  nine) 

784  counterstained  with  cresyl  violet  and  high  resolution  image  of  the  silver  labeled  site  (callout).  B: 

785  Camera  Lucida-style  drawing  of  the  cervical  spinal  cord  summarizing  the  anatomical  locations 

786  of  eight  silver  labeled  sites  obtained  in  one  animal.  C:  Integrated  phrenic  motor  output  and  mid- 

787  cervical  spinal  discharge  on  the  left  and  right  hemi-cord  corresponding  to  identified  silver 

788  labeled  sites  in  panel  B  during  baseline  and  hypoxia.  Electrode  one  recorded  phrenic 

789  motoneuron  discharge  and  the  remaining  electrodes  recorded  interneuron  discharge.  D: 

790  Representative  C4  section  summarizing  the  total  number  of  silver  labeled  interneurons  within 

791  each  lamina  in  spinal  intact  animals.  E:  Average  discharge  frequency  (Hz)  of  spinal 

792  interneurons  within  each  lamina  represented  in  panel  D.  CC:  central  canal;  DH:  dorsal  horn;  VH: 

793  ventral  horn.  Scale  bars:  A,  B:  0.5  mm  and  50  pm  (callout);  C:  1.5  sec. 


794  Figure  9:  Cross-correlation  analysis  of  mid-cervical  spinal  interneurons.  A: 

795  Representative  cross-correlations  obtained  from  three  pairs  of  neurons  depicting  a  significant 

796  central  peak  (left),  1.8  ms  offset  peak  (middle),  and  7.8  ms  offset  peak  (right).  The  number  of 

797  trigger  spikes  for  each  correlation:  43,891;  39,208;  3,928;  and  the  detectability  index  of  each 

798  correlation:  45.2;  4.3;  5.3,  respectively.  Black  lines  plotted  at  0.5  and  at  2.9  ms  to  indicate  how 

799  data  was  grouped  by  latency  (see  panel  B).  B:  A  histogram  of  the  latency  relative  to  the  trigger 

800  for  all  correlations  with  a  significant  peak.  The  first  bar  in  the  histogram  is  a  count  of  the  central 

801  peaks  (0  ms  latency)  and  each  successive  bar  in  the  plot  represents  counts  obtained  in  an 

802  increment  of  0.2.  Yellow  bars  plotted  at  0.5  and  at  2.9  ms  indicate  how  data  was  grouped  based 

803  on  latency.  Black  lines  are  2nd  order  polynomial  fit  of  the  data.  C:  The  sum  of  significant  cross- 

804  correlations  with  latencies  <  0.4  ms,  between  0.6-2. 8  ms,  and  >  3.0  ms.  Percentages  reflect  the 

805  proportion  of  positive  correlations  out  of  the  total  possible  (n=704).  D:  The  average  number  of 

806  positive  correlations  per  animal  when  both  recordings  were  on  the  same  side  of  the  cord  (20  ± 

807  3.9  positive  correlations/animal),  or  on  opposite  sides  of  the  cord  (25  ±  7.5  positive 

808  correlations/animal).  E:  The  number  of  significant  correlations  expressed  as  a  percent  of  total 

809  possible  connections  obtained  when  both  recordings  were  on  the  same  side  of  the  cord  (81/332 

810  total  possible),  or  on  opposite  sides  of  the  cord  (89/372  total  possible). 

811  Figure  10:  Anatomical  locations  of  functionally  connected  interneurons.  Summary  maps 

812  of  silver  labeled  interneurons  with  significant  peaks  in  cross-correlograms.  The  number  of 

813  positive  correlations  in  each  lamina  was  normalized  to  the  number  of  neurons  present  in  the 

814  lamina  and  shaded  accordingly.  A:  When  both  trigger  and  target  neurons  were  on  the  same 

815  side  of  the  spinal  cord  (e.g.,  unilateral  recordings),  interneurons  in  dorsal  lamina  (i.e.,  IV,  V  and 

816  VI)  made  and  received  the  greatest  number  of  excitatory  connections.  B:  When  the  trigger  and 

817  target  neurons  were  on  opposite  sides  of  the  spinal  cord  (e.g.,  bilateral  recordings), 

818  interneurons  in  lamina  V,  VIII  and  IX  made  and  received  the  most  connections. 


819 


Table  1:  Silver  staining  protocol.  Sequence  of  solutions  necessary  to  histologically  develop 

820  the  deposited  silver  adapted  from  (Spinelli  1975).  Solutions  were  prepared  and  arranged  in  the 

821  order  of  the  staining  sequence  in  either  glass  petri  dishes  or  plastic  trays.  Volumes  of  solution 

822  were  calculated  for  one  12-well  tray  of  tissue  (24  sections;  ~1  mm  of  tissue).  Shaded  text 

823  indicates  the  steps  in  which  solutions  were  changed  after  processing  one  tray  of  tissue,  while 

824  remaining  non-shaded  steps  indicate  solutions  were  chanced  after  processing  three  trays  of 

825  tissue. 
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Table  1. 


Dish 

no.  Contents 

Time 

Volume/1  tray 

1 

dH20 

5  min 

150mL/dish 

2 

dH20 

5  min 

150mL/dish 

3* 

1%  (NH4)2S 

5  min 

48mL/tray 

4*# 

0.125%  H202 

5  min 

48mL/tray 

5 

dH20 

5  min 

150mL/dish 

6 

dH20 

5  min 

150mL/dish 

7* 

1%  ascorbic  acid 

5  min 

48mL/tray 

8 

dH20 

5  min 

150mL/dish 

9 

dH20 

5  min 

150mL/dish 

10*# 

Developer 

8  min 

25%  acacia 

16mL 

2%  hydroquinone,  5%  citric  acid 

16mL 

“  48mL/tray 

1%  AgN02 

16mL 

11 

dH20 

Quick  rinse 

150mL/dish 

12 

dH20 

Quick  rinse 

150mL/dish 

13 

dH20 

Quick  rinse 

150mL/dish 

14 

1%  Na2S204.  2H20 

3-4  min 

150mL/dish 

15 

dH20 

Quick  rinse 

150mL/dish 

Store  in  dH20  at  4°C  until  ready  to  mount  on 

charged  slides,  then  counterstain,  and 

cover  with  mounting  medium. 

(*)  Place  solution  in  a  12-well  tray. 
(#)  Light  sensitive  solutions. 

□  Change  solution  every  3  times. 

□  Change  solutions  every  time. 


